CN114114050A - Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system - Google Patents
Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system Download PDFInfo
- Publication number
- CN114114050A CN114114050A CN202111311781.2A CN202111311781A CN114114050A CN 114114050 A CN114114050 A CN 114114050A CN 202111311781 A CN202111311781 A CN 202111311781A CN 114114050 A CN114114050 A CN 114114050A
- Authority
- CN
- China
- Prior art keywords
- value
- stress
- charge
- temperature
- parameter
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000012937 correction Methods 0.000 claims abstract description 113
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000002474 experimental method Methods 0.000 claims abstract description 11
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 9
- 238000010606 normalization Methods 0.000 claims description 10
- 238000012544 monitoring process Methods 0.000 claims description 7
- 230000015556 catabolic process Effects 0.000 claims description 6
- 238000006731 degradation reaction Methods 0.000 claims description 6
- 230000006399 behavior Effects 0.000 description 7
- 230000007613 environmental effect Effects 0.000 description 3
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 2
- 229910001416 lithium ion Inorganic materials 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/392—Determining battery ageing or deterioration, e.g. state of health
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
Abstract
The invention discloses a battery attenuation attribution method based on a rain flow counting method and a vehicle-mounted management systemPlossIs the percentage of the decay of the battery,is a first parameter, N is the number of cycles,as the second parameter, the parameter is,is a third parameter; setting a plurality of conditions, carrying out cyclic charge and discharge experiments on a plurality of batteries with the same model under each condition,obtaining a plurality of condition correction values; the third parameter is adjusted according to a plurality of condition correction valuesCorrecting to obtain a corrected attenuation reference curve; and processing the condition correction values and the corrected attenuation reference curve by using a rain flow counting method, and outputting the proportion of each stress in the total stress. The technical scheme provided by the invention can enable the user to adjust the behavior of the user according to each occupation situation, thereby enabling the service life of the battery to be longer.
Description
Technical Field
The invention relates to the field of power batteries, in particular to a battery attenuation attribution method based on a rain flow counting method and a vehicle-mounted management system.
Background
With the development of new energy industry, more and more users use electric vehicles. The attribution of the attenuation factor of the power battery is particularly important as the almost most expensive equipment on board the vehicle.
At present, the mainstream judgment methods for the attenuation reasons of the battery performance are all performance-based judgment methods, and little attention is paid to analysis of environmental factors and driving behaviors. For example, although patent application No. 201510332424.2 proposes a method for testing and diagnosing causes of performance degradation of lithium ion batteries, it is necessary to determine the causes of performance degradation of the batteries by comparing the absolute values of the ohmic impedance increase rate, the relaxation impedance increase rate, and the temperature entropy coefficient increase rate of the batteries to be tested with respect to the reference battery. Therefore, there is a need for a battery attenuation attribution method and an on-board management system based on the rain flow counting method.
Disclosure of Invention
The invention aims to provide a battery attenuation attribution method based on a rain flow counting method and an on-board management system, wherein the reason of the battery attenuation is analyzed according to environmental factors and driving behaviors.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a method for attributing battery attenuation based on rain flow counting, comprising:
carrying out cyclic charge-discharge experiment on a plurality of batteries with the same model to obtain the cycle number and the attenuation percentage of each battery, thereby obtaining the cycle number and the attenuation percentage of each batteryDetermining a decay reference curve for the battery of that type In the formula, PlossIs the percentage of the decay of the battery,is a first parameter, N is the number of cycles,as the second parameter, the parameter is,is a third parameter;
setting a plurality of conditions, wherein the conditions comprise temperature, depth of discharge, average state of charge and service life, and carrying out cyclic charge-discharge experiments on a plurality of batteries of the same type under each condition to obtain a plurality of condition correction values, wherein the condition correction values comprise a temperature correction value, a depth of discharge correction value, an average state of charge correction value and a calendar time correction value;
correcting the value according to a plurality of conditions and by the formula fd,1=[Sδ(δ)+St(tc)]Sσ(σ)ST(Tc) For the third parameterCorrecting to obtain a corrected attenuation reference curve, wherein fd,1Is a third parameterCorrection value of (S)δ(delta) is a depth of discharge correction value, St(tc) For calendar time correction values, Sσ(σ) is the mean state of charge correction value, ST(Tc) Is a temperature correction value;
and processing the condition correction values and the corrected attenuation reference curve by using a rain flow counting method, and outputting the proportion of each stress in the total stress, wherein the stress comprises temperature stress, depth of discharge stress, average state of charge stress and duration stress. The temperature stress is also the average temperature stress, which reflects the approximate temperature condition of the battery used by the user, and the average state of charge stress and the depth of discharge stress reflect the use habit of the user.
Further, before the determining the attenuation reference curve of the battery of the model, the method further comprises the following steps: fitting curves according to cycle number and percent decay for each cell Wherein, PlossIs the percent fade, alpha, of the cellseiIs a first parameter, N is the number of cycles, βseiIs a second parameter, fd,1Is a third parameter; solving alpha corresponding to each battery according to the curve obtained by fittingsei、βseiAnd fd,1And respectively carrying out mean value calculation to obtainAnd
further, the temperature correction value ST(T) is obtained by the following formula:
in the formula, ST(T) is a temperature correction value at temperature T, kTAs a parameter of temperature, TrefIs a reference temperature, ST(Tc) To be at a temperature TcCorrection value of temperature of fd,t[T=Tc]To be at a temperature TcThe value of the third parameter below is,to be at temperatureA third parameter value of; wherein f isd,t[T=Tc]Anddetermined by the decay reference curve.
Further, the average state of charge correction value Sσ(σ) is obtained by the following formula:
in the formula, Sσ(σ) is the mean state of charge correction value at mean state of charge σ, kσIs the average state of charge parameter, σrefTo reference the average state of charge, Sσ(σA) To be at an average state of charge σACorrected value of the average state of charge of fd,t[σ=σA]To be at an average state of charge σAValue of a third parameter ofd,t[σ=σref]To be at a reference average state of charge σrefA third parameter value of; wherein f isd,t[σ=σA]And fd,t[σ=σref]Determined by the decay reference curve.
Further, the calendar time correction value St(t) is obtained by the following formula:
St(t)=ktt
in the formula, St(t) is a calendar time correction value at calendar time t, ktAs a calendar time parameter, fd,t[T=TA,σ=σA]To be at a temperature TAAverage state of charge σAValue of the third parameter, ST(TA) To be at a temperature TACorrection value of temperature of Sσ(σA) To be at an average state of charge σAA lower average state of charge correction value; wherein f isd,t[T=TA,σ=σA]Determined by the decay reference curve.
Further, the depth of discharge correction value Sδ(δ) is obtained by the following formula:
in the formula, Sδ(delta) is a depth of discharge correction value at depth of discharge delta, kδ,e1Is a first depth of discharge parameter that is,as a second depth of discharge parameter, Sδ(δA) To a depth of discharge deltaAA correction value for the depth of discharge at the time of discharging,at depth of discharge deltaAAverage state of charge σATemperature TACalendar time tp,AValue of the third parameter, ST(TA) To be at a temperature TACorrection value of temperature of Sσ(σA) To be at an average state of charge σAA lower average state of charge correction value; wherein,determined by the decay reference curve.
Further, the processing the plurality of condition correction values and the corrected attenuation reference curve by using a rain flow counting method includes: setting an upper bound and a lower bound for each stress, and normalizing each stress value according to the following formula:
in the formula, SNormalizationIs normalized stress value, x is temperature or depth of discharge or calendar time or average state of charge, S (x) is stress value corresponding to x, and the stress value is corrected by the temperatureT(T) or depth of discharge correction value Sδ(delta) or calendar time correction values St(tc) Or mean state of charge correction value Sσ(sigma) determining that S (upper bound) is the upper bound stress value of the stress corresponding to x, and S (lower bound) is the lower bound stress value of the stress corresponding to x.
Further, the processing the plurality of condition correction values and the corrected attenuation reference curve by using a rain flow counting method further includes: calculating a total stress value according to the normalized stress values, and calculating the proportion of each stress value in the total stress value;
wherein the total stress value SGeneral assemblyCalculated by the following formula:
in the formula, SGeneral assemblyIs the total stress value, ST(T) is the temperature stress value, ST(upper bound) is the temperature upper bound stress value, ST(lower bound) is the value of the temperature lower bound stress, Sσ(σ) is the average state of charge stress value, Sσ(Upper)Bound) is the average state of charge upper bound stress value, Sσ(lower bound) is the value of the stress at the lower bound of the mean state of charge, St(t) is the value of the stress for the duration, St(upper bound) is the value of the stress at the upper bound of the duration, St(lower bound) is the value of the stress at the lower bound of the duration, Sδ(delta) is the depth of discharge stress value, Sδ(upper bound) is the stress value at the upper bound of the depth of discharge, Sδ(lower bound) is the stress value of the lower bound of the depth of discharge;
the proportion of each stress value to the total stress value is calculated by the following formula:
p=Snormalization/SGeneral assembly
Wherein p is the ratio of stress to total stress, SNormalizationIs a normalized stress value, SGeneral assemblyIs the total stress value.
Preferably, the upper bound of the temperature is 70 ℃, the lower bound of the temperature is-20 ℃, the upper bound of the average state of charge is 100%, and the lower bound of the average state of charge is 0%.
An on-vehicle management system based on the battery attenuation attribution method described above, for processing usage data of an on-vehicle battery and outputting an analysis result, the on-vehicle battery management system comprising:
a monitoring unit configured to acquire usage data of an on-vehicle battery, wherein the usage data includes temperature, depth of discharge, average state of charge, and usage duration data;
and the control unit is electrically connected with the monitoring unit, a data processing unit is arranged in the control unit, and the data processing unit is configured to process the acquired use data of the battery and output an analysis result.
The invention has the advantages that: the proportion of the stresses such as temperature stress, depth of discharge stress, average state of charge stress, duration stress and the like is obtained through analysis, so that a user can adjust own behavior according to each proportion condition, and the service life of the battery is longer.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart of a battery attenuation cause method according to an embodiment of the present invention.
Detailed Description
In order to make the technical solutions of the present invention better understood and more clearly understood by those skilled in the art, the technical solutions of the embodiments of the present invention will be described below in detail and completely with reference to the accompanying drawings. It should be noted that the implementations not shown or described in the drawings are in a form known to those of ordinary skill in the art. Additionally, while exemplifications of parameters including particular values may be provided herein, it is to be understood that the parameters need not be exactly equal to the respective values, but may be approximated to the respective values within acceptable error margins or design constraints. It is to be understood that the described embodiments are merely exemplary of a portion of the invention and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention. In addition, the terms "comprises" and "comprising," and any variations thereof, in the description and claims of this invention, are intended to cover a non-exclusive inclusion, such that a process, method, apparatus, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
In one embodiment of the invention, a battery attenuation attribution method based on a rain flow counting method is provided, and comprises the following steps:
firstly, carrying out cyclic charge and discharge (attenuation) experiments on a plurality of batteries of the same type to obtain the cycle number and the attenuation percentage of each battery, thereby determining the attenuation reference curve of the battery of the type.
Specifically, n batteries of the same type are taken and are heated at a standard temperature t0The cyclic charge and discharge experiment was performed, wherein n>3,t0At 25 ℃, the average state of charge (soc), i.e. soc cycle range, is 0% -100%, but the scope of the present invention is not limited thereto, and the percentage of decay P corresponding to each cycle of the number of cycles N is recordedlossAnd fitting a curve according to the cycle number and the attenuation percentage of each battery by the following formula:
in the formula, PlossIs the percent fade, alpha, of the cellseiIs a first parameter, which is related to the first effect of the battery and represents the proportion of lithium ions consumed by the battery at an early stage for forming a stable sei film, N is the number of cycles, βseiIs a second parameter, fd,1Is the third parameter. Then, α corresponding to each cell is determinedsei、βseiAnd fd,1And respectively carrying out mean value calculation to obtainAndthus determining the attenuation reference curve for this model of battery:
in the formula, PlossIs the percentage of the decay of the battery,is a first parameter, N is the number of cycles,as the second parameter, the parameter is,is the third parameter.
Secondly, setting a plurality of conditions, wherein the conditions comprise temperature, depth of discharge, average state of charge and service life, and carrying out cyclic charge-discharge experiments on a plurality of batteries of the same type under each condition to obtain a plurality of condition correction values, wherein the condition correction values comprise a temperature correction value, a depth of discharge correction value, an average state of charge correction value and a calendar time correction value; according to a plurality of condition correction values, and the third parameter is corrected by the following formulaCorrection is performed to obtain a corrected attenuation reference curve:
fd,1=[Sδ(δ)+St(tc)]Sσ(σ)ST(Tc)
in the formula (f)d,1Is a third parameterCorrection value of (S)δ(delta) is a depth of discharge correction value, St(tc) For calendar time correction values, Sσ(σ) is the mean state of charge correction value, ST(Tc) Is a temperature correction value.
Specifically, a plurality of batteries with the same model are taken, the batteries are divided into a plurality of groups, and a control variable method is adopted to perform cyclic charge and discharge experiments on the batteries at different temperatures, different average socs and different discharge depths, wherein at least 2 different temperatures, 2 different average socs and 4 different discharge depths are taken, but the protection scope of the invention is not limited by the above. Then, according to the cycle times N obtained by carrying out cycle charge and discharge experiments under different experimental conditions and the attenuation percentage P corresponding to each cyclelossTo determine the attenuation reference curve and then obtain fd,1OrThen, for the temperature correction value ST(T) has the following formula:
in the formula, ST(Tc) To be at a temperature TcCorrection value of temperature of fd,t[T=Tc]To be at a temperature TcThe value of the third parameter below is,to be at temperatureThe third parameter value of. Wherein f isd,t[T=Tc]Andf obtained from the attenuation reference curved,1OrIs determined so as to obtain a temperature TcCorrection value S of temperatureT(Tc) Then substituted into S as shown belowT(T) solving the equation to obtain kT:
In the formula, ST(T) is a temperature correction value at temperature T, kTAs a parameter of temperature, TrefAs the reference temperature, in the present embodiment, the reference temperature
Similarly to the temperature correction value STTemperature parameter k in (T)TFor mean state of charge correction values Sσ(σ) has the followingThe following formula:
in the formula, Sσ(σA) To be at an average state of charge σACorrected value of the average state of charge of fd,t[σ=σA]To be at an average state of charge σAValue of a third parameter ofd,t[σ=σref]To be at a reference average state of charge σrefThe third parameter value of. Wherein f isd,t[σ=σA]And fd,t[σ=σref]F obtained from the attenuation reference curved,1OrDetermination to obtain the average state of charge σALower average state of charge correction value Sσ(σA) Then substituted into S as shown belowσSolving formula of (sigma) to obtain kσ:
In the formula, Sσ(σ) is the mean state of charge correction value at mean state of charge σ, kσIs the average state of charge parameter, σrefTo refer to the average state of charge, in the present embodiment, the average state of charge σ is referred torefIs avg (soc)cycle)。
Similarly to the temperature correction value STTemperature parameter k in (T)TFor calendar time correction values St(t) has the following formula:
St(t)=ktt (5)
in the formula, ktIs a calendar time parameter, which is determined by:
in the formula (f)d,t[T=TA,σ=σA]To be at a temperature TAAverage state of charge σAA third parameter value of f obtained from the attenuation reference curved,1OrDetermining; t is the current calendar time; sT(TA) To be at a temperature TAThe temperature correction value of (c) is obtained from the equation (2), and it should be noted that, in the present embodiment, k in the equation (2)TFound and TrefIs t0Bringing in temperature TACan calculate ST(TA);Sσ(σA) To be at an average state of charge σAA lower average state of charge correction value derived from equation (4), and k in equation (4)σFound and σrefIs avg (soc)cycle) Bringing into an average state of charge σACan calculate Sσ(σA) Thereby f will bed,t[T=TA,σ=σA]Current calendar time t, ST(TA) And Sσ(σA) Substituting into formula (6) to obtain calendar time parameter ktThen k is addedtSubstituting into equation (5) allows determination of a calendar time correction value S at a calendar time tt(t)。
Similarly to the temperature correction value STTemperature parameter k in (T)TFor a correction value S of depth of dischargeδ(δ) has the following formula:
in the formula, Sδ(δA) To a depth of discharge deltaAA lower depth of discharge correction value; sT(TA) To be at a temperature TACorrection value of temperature of Sσ(σA) To be at an average state of charge σAAverage charge ofState correction value, St(tp,A) To be at calendar time tp,ACorrection of calendar time ST(TA)、Sσ(σA)、St(tp,A) Obtained from the above formulas (2), (4), (5) for which the respective parameters have been found;to a depth of discharge deltaAAverage state of charge σATemperature TACalendar time tp,AA third parameter value of f obtained from the attenuation reference curved,1OrDetermined and then substituted into S as shown belowδSolving the equation of (delta) to obtain kδ,e1、
In the formula, Sδ(delta) is a depth of discharge correction value at depth of discharge delta, kδ,e1Is a first depth of discharge parameter that is,is a second depth of discharge parameter.
After determining parameters in a formula solved by each correction value, obtaining a plurality of condition correction values and a corrected attenuation reference curve through experiments according to different set conditions, namely different temperatures, different average socs and different depths of discharge, then processing the condition correction values and the corrected attenuation reference curve by using a rain flow counting method, and outputting the proportion of a plurality of stresses in the total stress, wherein the stresses comprise temperature stress, depth of discharge stress, average state-of-charge stress and duration stress, the temperature stress reflects the approximate temperature condition of a user such as a driver using a battery, the average state-of-charge stress and the depth of discharge stress reflect the use habit of the driver, the duration stress reflects the duration of the driver using the battery, it needs to be noted that the temperature stress is analyzed based on the temperature correction value, and the depth of discharge stress is analyzed based on the depth of discharge correction value, the average state of charge stress is analyzed based on the average state of charge correction value, and the duration stress is analyzed based on the calendar time correction value.
Specifically, when a plurality of condition correction values and the modified attenuation reference curve are processed by a rain flow counting method, an upper bound and a lower bound are set for each stress, and each stress value is normalized according to the following formula:
in the formula, SNormalizationIs normalized stress value, x is temperature or depth of discharge or calendar time or average state of charge, S (x) is stress value corresponding to x, and the stress value is corrected by the temperatureT(T) or depth of discharge correction value Sδ(delta) or calendar time correction values St(tc) Or mean state of charge correction value Sσ(sigma) determining that S (upper bound) is the upper bound stress value of the stress corresponding to x, and S (lower bound) is the lower bound stress value of the stress corresponding to x.
Taking temperature stress as an example, in this embodiment, 70 ℃ is set as an upper limit of temperature, and-20 ℃ is set as a lower limit of temperature, there is normalization of temperature stress:
the depth of discharge, average state of charge, and duration stress are similar to the temperature stress and are not described herein.
Calculating a total stress value S according to the following formula and the normalized temperature stress value, the normalized depth of discharge stress value, the normalized average state of charge stress value and the normalized duration stress valueGeneral assembly:
In the formula, SGeneral assemblyIs the total stress value, ST(T) is the temperature stress value, i.e. the temperature correction value as described above, ST(upper bound) is the temperature upper bound stress value, ST(lower bound) is the value of the temperature lower bound stress, Sσ(σ) is the mean state of charge stress value, i.e. the mean state of charge correction value, Sσ(upper bound) is the average state of charge upper bound stress value, Sσ(lower bound) is the value of the stress at the lower bound of the mean state of charge, St(t) is the duration stress value, i.e. the calendar time correction value as described above, St(upper bound) is the value of the stress at the upper bound of the duration, St(lower bound) is the value of the stress at the lower bound of the duration, Sδ(delta) is the depth of discharge stress value, i.e. the depth of discharge correction value, Sδ(upper bound) is the stress value at the upper bound of the depth of discharge, Sδ(lower bound) is the depth of discharge lower bound stress value.
And then calculating the proportion of each stress value to the total stress value through the following formula and outputting the proportion:
p=Snormalization/SGeneral assembly
Wherein p is the ratio of stress to total stress, SNormalizationIs a normalized stress value, SGeneral assemblyIs the total stress value. Based on the analysis, the user can know the use behavior of the user according to each stress proportion, and the user can adjust the use behavior of the user, so that the service life of the battery is longer.
In this embodiment, the upper bound of the average state of charge is 100%, the lower bound of the average state of charge is 0%, and the upper bound and the lower bound of each stress are selected according to actual conditions, which does not limit the protection scope of the present invention. In addition, the temperature T iscAverage state of charge σADepth of discharge deltaACalendar time tp,AAre illustrative and do not limit the scope of the invention by their numerical values or indices.
In an embodiment of the present invention, an on-board management system based on the battery attenuation attribution method described above is provided, configured to process usage data of an on-board battery and output an analysis result, the on-board battery management system includes a monitoring unit and a control unit, wherein the monitoring unit is configured to acquire the usage data of the on-board battery, the usage data includes temperature, depth of discharge, average state of charge and usage duration data; the control unit is electrically connected with the monitoring unit, a data processing unit is arranged in the control unit, and the data processing unit is configured to process the acquired use data of the battery and output an analysis result.
The idea of the embodiment of the vehicle-mounted management system and the working process of the battery attenuation attribution method in the embodiment belong to the same idea, and the whole content of the embodiment of the battery attenuation attribution method is incorporated into the embodiment of the vehicle-mounted management system in a full-text reference manner, so that the detailed description is omitted.
The technical scheme provided by the invention is a judgment method based on environmental factors and driving behaviors, which is beneficial for users to adjust own use behaviors and enables the service life of a battery to be longer.
The above description is only for the preferred embodiment of the present invention and is not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes that can be directly or indirectly applied to other related technical fields using the contents of the present specification and the accompanying drawings are included in the scope of the present invention.
Claims (10)
1. A method for attributing battery attenuation based on a rain flow counting method, comprising:
carrying out cyclic charge and discharge experiments on a plurality of batteries with the same model to obtain the cycle number and the attenuation percentage of each battery, thereby determining the attenuation reference curve of the battery with the model, namely In the formula, PlossIs the percentage of the decay of the battery,is a first parameter, N is the number of cycles,as the second parameter, the parameter is,is a third parameter;
setting a plurality of conditions, wherein the conditions comprise temperature, depth of discharge, average state of charge and service life, and carrying out cyclic charge-discharge experiments on a plurality of batteries of the same type under each condition to obtain a plurality of condition correction values, wherein the condition correction values comprise a temperature correction value, a depth of discharge correction value, an average state of charge correction value and a calendar time correction value;
correcting the value according to a plurality of conditions and by the formula fd,1=[Sδ(δ)+St(tc)]Sσ(σ)ST(Tc) For the third parameterCorrecting to obtain a corrected attenuation reference curve, wherein fd,1Is a third parameterCorrection value of (S)δ(delta) is a depth of discharge correction value, St(tc) For calendar time correction values, Sσ(σ) is the mean state of charge correction value, ST(Tc) Is a temperature correction value;
and processing the condition correction values and the corrected attenuation reference curve by using a rain flow counting method, and outputting the proportion of each stress in the total stress, wherein the stress comprises temperature stress, depth of discharge stress, average state of charge stress and duration stress.
2. The method for attributing battery attenuation according to claim 1, further comprising, before the determining the attenuation reference curve for the battery of the model: fitting curves according to cycle number and percent decay for each cellWherein, PlossIs the percent fade, alpha, of the cellseiIs a first parameter, N is the number of cycles, βseiIs a second parameter, fd,1Is a third parameter; solving alpha corresponding to each battery according to the curve obtained by fittingsei、βseiAnd fd,1And respectively carrying out mean value calculation to obtain And
3. the battery degradation attributable method according to claim 1, wherein the temperature correction value ST(T) is obtained by the following formula:
in the formula, ST(T) is a temperature correction value at temperature T, kTAs a parameter of temperature, TrefIs a reference temperature, ST(Tc) To be at a temperature TcCorrection value of temperature of fd,t[T=Tc]To be at a temperature TcThe value of the third parameter below is,to be at temperatureA third parameter value of; wherein f isd,t[T=Tc]Anddetermined by the decay reference curve.
4. The battery degradation attributable method of claim 1 wherein the average state of charge correction value Sσ(σ) is obtained by the following formula:
in the formula, Sσ(σ) is the mean state of charge correction value at mean state of charge σ, kσIs the average state of charge parameter, σrefTo reference the average state of charge, Sσ(σA) To be at an average state of charge σACorrected value of the average state of charge of fd,t[σ=σA]To be at an average state of charge σAValue of a third parameter ofd,t[σ=σref]To be at a reference average state of charge σrefA third parameter value of; wherein f isd,t[σ=σA]And fd,t[σ=σref]Determined by the decay reference curve.
5. The battery degradation attribution method according to claim 1, wherein the calendar time correction value St(t) is obtained by the following formula:
St(t)=ktt
in the formula, St(t) is a calendar time correction value at calendar time t, ktAs a calendar time parameter, fd,t[T=TA,σ=σA]To be at a temperature TAAverage state of charge σAValue of the third parameter, ST(TA) To be at a temperature TACorrection value of temperature of Sσ(σA) To be at an average state of charge σAA lower average state of charge correction value; wherein f isd,t[T=TA,σ=σA]Determined by the decay reference curve.
6. The battery degradation attributable method of claim 1 wherein the depth of discharge correction value Sδ(δ) is obtained by the following formula:
in the formula, Sδ(delta) is a depth of discharge correction value at depth of discharge delta, kδ,e1Is a first depth of discharge parameter that is,as a second depth of discharge parameter, Sδ(δA) To a depth of discharge deltaAIs as followsA correction value for the depth of discharge,at depth of discharge deltaAAverage state of charge σATemperature TACalendar time tp,AValue of the third parameter, ST(TA) To be at a temperature TACorrection value of temperature of Sσ(σA) To be at an average state of charge σAA lower average state of charge correction value; wherein,determined by the decay reference curve.
7. The battery attenuation attribution method according to claim 1, wherein the processing of the plurality of condition correction values and the corrected attenuation reference curve using a rain flow counting method comprises: setting an upper bound and a lower bound for each stress, and normalizing each stress value according to the following formula:
in the formula, SNormalizationIs normalized stress value, x is temperature or depth of discharge or calendar time or average state of charge, S (x) is stress value corresponding to x, and the stress value is corrected by the temperatureT(T) or depth of discharge correction value Sδ(delta) or calendar time correction values St(tc) Or mean state of charge correction value Sσ(sigma) determining that S (upper bound) is the upper bound stress value of the stress corresponding to x, and S (lower bound) is the lower bound stress value of the stress corresponding to x.
8. The battery attenuation attribution method according to claim 7, wherein the processing the plurality of condition correction values and the corrected attenuation reference curve using a rain flow counting method further comprises: calculating a total stress value according to the normalized stress values, and calculating the proportion of each stress value in the total stress value;
wherein the total stress value SGeneral assemblyCalculated by the following formula:
in the formula, SGeneral assemblyIs the total stress value, ST(T) is the temperature stress value, ST(upper bound) is the temperature upper bound stress value, ST(lower bound) is the value of the temperature lower bound stress, Sσ(σ) is the average state of charge stress value, Sσ(upper bound) is the average state of charge upper bound stress value, Sσ(lower bound) is the value of the stress at the lower bound of the mean state of charge, St(t) is the value of the stress for the duration, St(upper bound) is the value of the stress at the upper bound of the duration, St(lower bound) is the value of the stress at the lower bound of the duration, Sδ(delta) is the depth of discharge stress value, Sδ(upper bound) is the stress value at the upper bound of the depth of discharge, Sδ(lower bound) is the stress value of the lower bound of the depth of discharge;
the proportion of each stress value to the total stress value is calculated by the following formula:
p=Snormalization/SGeneral assembly
Wherein p is the ratio of stress to total stress, SNormalizationIs a normalized stress value, SGeneral assemblyIs the total stress value.
9. The method of claim 8, wherein the upper bound of temperature is 70 ℃, the lower bound of temperature is-20 ℃, the upper bound of average state of charge is 100%, and the lower bound of average state of charge is 0%.
10. An in-vehicle management system for processing usage data of an in-vehicle battery and outputting an analysis result based on the battery attenuation attribution method of claim 1, the in-vehicle battery management system comprising:
a monitoring unit configured to acquire usage data of an on-vehicle battery, wherein the usage data includes temperature, depth of discharge, average state of charge, and usage duration data;
and the control unit is electrically connected with the monitoring unit, a data processing unit is arranged in the control unit, and the data processing unit is configured to process the acquired use data of the battery and output an analysis result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111311781.2A CN114114050A (en) | 2021-11-08 | 2021-11-08 | Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111311781.2A CN114114050A (en) | 2021-11-08 | 2021-11-08 | Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114114050A true CN114114050A (en) | 2022-03-01 |
Family
ID=80381477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111311781.2A Pending CN114114050A (en) | 2021-11-08 | 2021-11-08 | Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114114050A (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012166489A2 (en) * | 2011-05-31 | 2012-12-06 | Caterpillar Inc. | Methods and systems for estimating battery health |
CN112986828A (en) * | 2021-04-12 | 2021-06-18 | 上海电享信息科技有限公司 | Method for estimating health degree of battery |
-
2021
- 2021-11-08 CN CN202111311781.2A patent/CN114114050A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012166489A2 (en) * | 2011-05-31 | 2012-12-06 | Caterpillar Inc. | Methods and systems for estimating battery health |
CN112986828A (en) * | 2021-04-12 | 2021-06-18 | 上海电享信息科技有限公司 | Method for estimating health degree of battery |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108919129B (en) | Service life prediction method of power battery under time-varying working condition | |
CN103869256B (en) | Method for estimating SOH of power lithium ion battery based on alternating current impedance test | |
Tian et al. | Consistency evaluation and cluster analysis for lithium-ion battery pack in electric vehicles | |
CN108919137B (en) | A kind of battery aging status estimation method considering different battery status | |
CN104142477B (en) | The battery condition estimator of solid-state concentration models coupling experience equivalent-circuit model | |
Takyi‐Aninakwa et al. | A strong tracking adaptive fading‐extended Kalman filter for the state of charge estimation of lithium‐ion batteries | |
CN103424710B (en) | For monitoring the method and system that the performance of the aged monomer in set of cells changes | |
CN110795851B (en) | Lithium ion battery modeling method considering environmental temperature influence | |
US20110231122A1 (en) | Method and system for determining the kind of a battery | |
CN112098853A (en) | Capacity attenuation battery discharge power determination method and device | |
Soto et al. | Analysis of the main battery characterization techniques and experimental comparison of commercial 18650 Li-ion cells | |
CN110045291B (en) | Lithium battery capacity estimation method | |
CN108828448B (en) | Battery state of charge on-line estimation method based on charging voltage curve fusion Kalman filtering | |
CN104681851B (en) | A kind of vapour vehicle lithium-ion power battery method for group matching | |
CN112526352B (en) | SOH estimation method for retired lithium ion battery | |
CN115639481B (en) | Battery data preprocessing system and method based on big data prediction SOC | |
CN114264964B (en) | Method, device, equipment and medium for evaluating battery capacity | |
CN106180002B (en) | A kind of battery screening technique utilized for echelon | |
CN112485695A (en) | Detection method and device for power battery | |
CN115219918A (en) | Lithium ion battery life prediction method based on capacity decline combined model | |
Wang et al. | Parameter sensitivity analysis and parameter identifiability analysis of electrochemical model under wide discharge rate | |
CN114355200A (en) | Method and device for establishing battery allowable current model and battery management system | |
CN109716151B (en) | Method and apparatus for estimating battery voltage | |
CN112986828A (en) | Method for estimating health degree of battery | |
CN114114050A (en) | Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |