Background
The extremely-fast self-heating technology realized by the external short circuit of the power battery is an important solution for overcoming the poor low-temperature performance of the power battery and realizing the all-weather application at present. Because the lithium ion battery systems used in electric vehicles contain hundreds of lithium ion battery cells, a single battery failure can generate a large amount of heat, triggering thermal runaway of adjacent batteries, and causing the entire battery pack to fail. Particularly, when the battery is in an external short circuit, the temperature of the power battery rises rapidly, and the phenomena of air leakage, liquid leakage, even combustion and the like can occur in a short time, thus threatening the safety of the electric carrying tool. Therefore, the external short circuit behavior of the power battery needs to be accurately depicted, an accurate hot spot coupling model is established, the temperature of the battery with the external short circuit is accurately predicted, and the safety of the battery in the extremely-fast heating process is further improved.
At present, the research on the prediction of the battery temperature is mostly focused on the normal working state, and the research on the external short circuit of the battery is insufficient. The battery thermoelectric coupling model applied to the battery management system not only needs to have higher precision, but also needs to be convenient to operate and simple in parameter setting as much as possible. The existing electrochemical thermal coupling model based on partial differential equation has higher accuracy for predicting the battery temperature, but real-time calculation and parameterization are not easy to realize due to fussy parameter setting and larger calculation amount. Therefore, to balance complexity and temperature accuracy, lumped parameter thermal models based on equivalent circuit models are widely used in battery management systems. These models are applied to temperature prediction in a normal operating state, and have good accuracy in a certain temperature range, but are not applicable because the temperature span is large and the distribution is uneven when the external short circuit occurs in the battery. Therefore, how to establish a thermocouple die type specially aiming at the external short circuit characteristic of the power battery has important significance.
Disclosure of Invention
Aiming at the problems, the invention provides a modeling method for the thermoelectric coupling characteristics of a lithium ion power battery, which is suitable for the extremely-fast self-heating process of the power battery. The method specifically comprises the following steps:
establishing a power battery model considering the bidirectional coupling effect of an electric model and a thermal model of a lithium ion battery aiming at the power battery;
secondly, identifying the internal resistance, the polarization internal resistance and the polarization capacitance of the power battery based on the electric model, and calculating the heat generation of the power battery by combining the identification result with the thermal model;
respectively constructing a radial uniformly-distributed one-dimensional thermal model and an axial linearly-distributed one-dimensional thermal model of the cylindrical battery with the convection boundary condition by using the obtained heat production result of the power battery;
and step four, establishing a two-dimensional temperature distribution model of the power battery according to the radial uniformly distributed one-dimensional thermal model and the axial linearly distributed one-dimensional thermal model, thereby completing the modeling of the thermoelectric coupling characteristic.
Further, the electrical model employs a first order RC cell model. The battery is heated by a plurality of heating sources, activation, concentration and ohmic loss. Under the external short circuit condition of battery, the battery heat production condition is complicated, and except the joule heat that battery internal resistance produced under normal conditions, it produces heat to have a large amount of utmost point ears still. The heat production of the tabs is generated by the current passing through the battery tabs, and the tabs are made of metal materials and have certain resistance. When the ESC occurs in the battery, a large current passes through the tabs, which generates a large amount of heat. Therefore, in this first order RC cell model, tab heat generation is also considered. The temperature and SOC of the battery change greatly in a short time when an external short circuit occurs, so the battery internal resistance, polarization internal resistance, and polarization capacitance can be regarded as functions related to the temperature and SOC.
Further, the power battery heat generation is calculated by the following formula:
wherein Q isavgAverage heat production rate of battery, QtabFor the heat production rate of the battery tab, I is the current, ROIs the ohmic internal resistance, R, of the battery1For polarizing internal resistance, C1Is polarization capacitance, t is time, t0For the moment when the battery starts to discharge, RtabThe resistance of a tab, T the temperature of the battery and SoC the state of charge of the battery; f. ofRO、fR1、fC1Respectively representing functions corresponding to ohmic internal resistance, polarization internal resistance and polarization capacitance.
Of course, the electrical and thermal models described above are equally applicable to many different forms known in the art.
Constructing a one-dimensional thermal model of a radially uniformly distributed, axially linear distribution of cylindrical cells with convective boundary conditions details the thermal distribution at external short circuits of the cell. Through the disassembly of the battery structure and the test results, the model is properly simplified, and the battery reasonably and uniformly generates heat along the radial direction and linearly generates heat along the axial direction under the influence of the lugs. Due to the fact that the thermal conductivity is one or two orders of magnitude higher in the axial direction than in the radial direction. The heat conduction of the battery in the axial direction is neglected.
Therefore, further, the established one-dimensional thermal model with radially uniform distribution and the one-dimensional thermal model with axially linear distribution are based on the following assumptions on the physical properties of the battery:
(1) the internal material of the battery has uniform isotropic physical properties and consistent density;
(2) the internal heat transfer of the cell is mainly through heat conduction, neglecting thermal radiation and thermal convection effects;
(3) the specific heat capacity and the thermal conductivity of the battery are constants and do not change with temperature and SOC.
Further, establishing a two-dimensional temperature distribution model of the power battery specifically comprises:
control equation based on three-dimensional temperature distribution T (z, r, T) and boundary conditions:
wherein z and r are the position coordinates of a certain point on the battery under the cylindrical coordinate system, and ktIs the radial thermal conductivity of the cell, rho is the cell density, CpIs the average specific heat capacity of the battery, R is the radius of the battery, h is the surface convection heat transfer coefficient, T∞Is the temperature of the fluid surrounding the battery, VbIs the cell volume;
assuming that the battery generates heat in the axial direction and linearly changes from the negative electrode to the positive electrode, the heat generation rate Q (z, t) of the battery is expressed as:
wherein QavgFor average heat production inside the cell, QtabHeat is generated for the tab part in the battery;
due to uniform heat generation, neglecting axial heat transfer, the temperature distribution of the battery cell in any axial r direction is assumed to satisfy the following polynomial equation:
wherein a (t) b (t) d (t) is a time-varying constant;
the cell center temperature Tc, the surface temperature Ts, the volume average temperature
And temperature gradient
Can be expressed as:
Tc=a(t)
Ts=a(t)+b(t)+d(t)
through a series of calculations, a temperature distribution model of the form:
where τ is the time constant, and the matrices A, B, C, D are defined as follows:
the method provided by the invention establishes a relation between the electric model and the thermal model of the lithium ion power battery, considers the specific parameters influencing the temperature change under the condition of external short circuit of the battery, can be better suitable for the extremely-fast self-heating process of the power battery, and obviously improves the accuracy of battery temperature prediction during external short circuit.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. 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.
The method provided by the invention specifically comprises the following steps:
step one, establishing a power battery model considering the bidirectional coupling effect of an electric model and a thermal model of a lithium ion battery aiming at a power battery, as shown in figure 1;
secondly, identifying the internal resistance, the polarization internal resistance and the polarization capacitance of the power battery based on the electric model, and calculating the heat generation of the power battery by combining the identification result with the thermal model;
respectively constructing a radial uniformly-distributed one-dimensional thermal model and an axial linearly-distributed one-dimensional thermal model of the cylindrical battery with the convection boundary condition by using the obtained heat production result of the power battery;
and step four, establishing a two-dimensional temperature distribution model of the power battery according to the radial uniformly distributed one-dimensional thermal model and the axial linearly distributed one-dimensional thermal model, thereby completing the modeling of the thermoelectric coupling characteristic.
In one example, for the proposed method for constructing the external short-circuit thermocouple model of the battery, the model is applied at the ambient temperature of 20 ℃, and the surface and center temperature data of the external short-circuit test of the power battery with three different initial SOCs and the corresponding model surface and center temperature data are compared and analyzed.
Because the center of the power lithium battery is limited by space and only one thermocouple can be placed, the center temperature of the model cathode is compared with the test. The result of the comparison between the experiment and the model is shown in fig. 2, the surface temperature of the power battery obtained by the model is basically consistent with the actually measured temperature and is influenced by uneven heat transfer, the experimental temperature curve has certain postposition, and the positive electrode is influenced more obviously.
As can be seen from the temperature rise errors in the following table, the maximum error is 2.315%, and the accuracy is reliable. Therefore, the established thermoelectric coupling model applied to the external short circuit condition of the cylindrical lithium ion battery is closer to the actual temperature change rule of the power battery after external short circuit, the effectiveness of the model is verified, and the model can be used for predicting the external short circuit temperature of the power battery.
It should be understood that, the sequence numbers of the steps in the embodiments of the present invention do not mean the execution sequence, and the execution sequence of each process should be determined by the function and the inherent logic of the process, and should not constitute any limitation on the implementation process of the embodiments of the present invention.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.