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 PDF

Info

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
Application number
CN202111311781.2A
Other languages
Chinese (zh)
Inventor
柯鹏
钱磊
朱卓敏
周玉
苏健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Powershare Information Technology Co ltd
Original Assignee
Shanghai Powershare Information Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Powershare Information Technology Co ltd filed Critical Shanghai Powershare Information Technology Co ltd
Priority to CN202111311781.2A priority Critical patent/CN114114050A/en
Publication of CN114114050A publication Critical patent/CN114114050A/en
Pending legal-status Critical Current

Links

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 system
Figure DDA0003342250380000011
PlossIs the percentage of the decay of the battery,
Figure DDA0003342250380000012
is a first parameter, N is the number of cycles,
Figure DDA0003342250380000013
as the second parameter, the parameter is,
Figure DDA0003342250380000014
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 values
Figure DDA0003342250380000015
Correcting 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

Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system
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
Figure BDA0003342250360000011
Figure BDA0003342250360000012
In the formula, PlossIs the percentage of the decay of the battery,
Figure BDA0003342250360000013
is a first parameter, N is the number of cycles,
Figure BDA0003342250360000014
as the second parameter, the parameter is,
Figure BDA0003342250360000015
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 parameter
Figure BDA0003342250360000021
Correcting to obtain a corrected attenuation reference curve, wherein fd,1Is a third parameter
Figure BDA0003342250360000022
Correction 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
Figure BDA0003342250360000023
Figure BDA0003342250360000024
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 obtain
Figure BDA0003342250360000025
And
Figure BDA0003342250360000026
further, the temperature correction value ST(T) is obtained by the following formula:
Figure BDA0003342250360000027
Figure BDA0003342250360000028
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,
Figure BDA0003342250360000029
to be at temperature
Figure BDA00033422503600000210
A third parameter value of; wherein f isd,t[T＝Tc]And
Figure BDA00033422503600000211
determined by the decay reference curve.
Further, the average state of charge correction value Sσ(σ) is obtained by the following formula:
Figure BDA00033422503600000212
Figure BDA00033422503600000213
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
Figure BDA0003342250360000031
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:
Figure BDA0003342250360000037
Figure BDA0003342250360000032
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,
Figure BDA0003342250360000033
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,
Figure BDA0003342250360000034
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,
Figure BDA0003342250360000035
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:
Figure BDA0003342250360000036
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:
Figure BDA0003342250360000041
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:
Figure BDA0003342250360000051
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 obtain
Figure BDA0003342250360000061
And
Figure BDA0003342250360000062
thus determining the attenuation reference curve for this model of battery:
Figure BDA0003342250360000063
in the formula, PlossIs the percentage of the decay of the battery,
Figure BDA0003342250360000064
is a first parameter, N is the number of cycles,
Figure BDA0003342250360000065
as the second parameter, the parameter is,
Figure BDA0003342250360000066
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 formula
Figure BDA0003342250360000067
Correction 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 parameter
Figure BDA0003342250360000068
Correction 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,1Or
Figure BDA0003342250360000069
Then, for the temperature correction value ST(T) has the following formula:
Figure BDA00033422503600000610
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,
Figure BDA00033422503600000611
to be at temperature
Figure BDA00033422503600000612
The third parameter value of. Wherein f isd,t[T＝Tc]And
Figure BDA00033422503600000613
f obtained from the attenuation reference curved,1Or
Figure BDA00033422503600000614
Is 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：
Figure BDA00033422503600000615
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
Figure BDA00033422503600000616
Similarly to the temperature correction value STTemperature parameter k in (T)TFor mean state of charge correction values Sσ(σ) has the followingThe following formula:
Figure BDA0003342250360000071
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,1Or
Figure BDA0003342250360000072
Determination 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σ：
Figure BDA0003342250360000073
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:
Figure BDA0003342250360000074
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,1Or
Figure BDA0003342250360000075
Determining; 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:
Figure BDA0003342250360000081
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;
Figure BDA0003342250360000082
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,1Or
Figure BDA0003342250360000083
Determined and then substituted into S as shown belowδSolving the equation of (delta) to obtain kδ,e1、
Figure BDA0003342250360000084
Figure BDA0003342250360000085
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,
Figure BDA0003342250360000086
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:
Figure BDA0003342250360000087
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:
Figure BDA0003342250360000091
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：
Figure BDA0003342250360000092
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
Figure FDA0003342250350000011
Figure FDA0003342250350000012
In the formula, PlossIs the percentage of the decay of the battery,
Figure FDA0003342250350000013
is a first parameter, N is the number of cycles,
Figure FDA0003342250350000014
as the second parameter, the parameter is,
Figure FDA0003342250350000015
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 parameter
Figure FDA0003342250350000016
Correcting to obtain a corrected attenuation reference curve, wherein fd，1Is a third parameter
Figure FDA0003342250350000017
Correction 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 cell
Figure FDA0003342250350000018
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 obtain
Figure FDA0003342250350000019
Figure FDA00033422503500000110
And
Figure FDA00033422503500000111
3. the battery degradation attributable method according to claim 1, wherein the temperature correction value ST(T) is obtained by the following formula:
Figure FDA00033422503500000112
Figure FDA0003342250350000021
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,
Figure FDA0003342250350000022
to be at temperature
Figure FDA0003342250350000023
A third parameter value of; wherein f isd，t[T＝Tc]And
Figure FDA0003342250350000024
determined 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:
Figure FDA0003342250350000025
Figure FDA0003342250350000026
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
Figure FDA0003342250350000027
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:
Figure FDA0003342250350000031
Figure FDA0003342250350000032
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,
Figure FDA0003342250350000033
as a second depth of discharge parameter, Sδ(δA) To a depth of discharge deltaAIs as followsA correction value for the depth of discharge,
Figure FDA0003342250350000034
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,
Figure FDA0003342250350000035
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:
Figure FDA0003342250350000036
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:
Figure FDA0003342250350000037
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.
CN202111311781.2A 2021-11-08 2021-11-08 Battery attenuation attribution method based on rain flow counting method and vehicle-mounted management system Pending CN114114050A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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