CN113189508A - Reverse calculation method and device for internal thermal resistance parameters of lithium ion soft package battery - Google Patents
Reverse calculation method and device for internal thermal resistance parameters of lithium ion soft package battery Download PDFInfo
- Publication number
- CN113189508A CN113189508A CN202110486991.9A CN202110486991A CN113189508A CN 113189508 A CN113189508 A CN 113189508A CN 202110486991 A CN202110486991 A CN 202110486991A CN 113189508 A CN113189508 A CN 113189508A
- Authority
- CN
- China
- Prior art keywords
- cel
- gas
- soft package
- thermal resistance
- temperature
- 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.)
- Granted
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/389—Measuring internal impedance, internal conductance or related variables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
- G01K7/02—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using thermoelectric elements, e.g. thermocouples
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/378—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC] specially adapted for the type of battery or accumulator
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Battery Mounting, Suspending (AREA)
Abstract
The invention relates to the technical field of lithium ion battery thermal management, discloses a reverse method and a reverse device for solving internal thermal resistance parameters of a lithium ion soft package battery, and aims to provide an estimation tool for key parameters for accurately calculating the temperature of the lithium ion soft package battery. Firstly, defining internal thermal resistance parameters to be solved, and constructing a thermal resistance parameter reverse solving experimental device to implement a thermal induction experiment; secondly, physical equations such as heat generation, gas generation, heat transfer, expansion and the like are established and fused into a soft package battery temperature evaluation model; and thirdly, constructing a reverse model of the internal thermal resistance parameter of the soft package battery, and solving based on the test data and the physical model to obtain the internal thermal resistance parameter of the battery. Compared with the prior art, the thermal resistance parameter obtained by the method is suitable for the lithium ion soft package battery with the same electrolyte, and the application range is larger; based on experimental data and a physical model, the thermal resistance parameter is calculated by a reverse solving technology, so that the precision is high; clear steps, easy understanding, simple device, clear definition and good engineering practicability.
Description
Technical Field
The invention relates to the technical field of lithium ion battery thermal management, in particular to a method and a device for reversely solving internal thermal resistance parameters of a lithium ion soft package battery.
Background
The lithium ion soft package battery is widely applied to consumer electronic equipment due to the advantages of high energy density, continuous reduction of cost and the like, and is gradually expanded to the fields of energy storage systems and power automobiles. But safety problems due to thermal runaway of the battery are always present. Under the condition of thermal abuse, a chemical exothermic reaction occurs inside the soft package battery; when the heat release rate exceeds the heat release efficiency, the internal temperature rise accelerates the exothermic reaction, resulting in the continuous accumulation of heat and reaction gas inside the battery. Reaching a critical value of temperature or pressure will trigger thermal runaway of the cell, releasing a significant amount of heat and combustible gases.
In order to effectively predict the thermal behavior of the lithium ion soft package battery, a thermal model needs to be constructed to accurately estimate the internal temperature of the lithium ion soft package battery.
The accuracy of temperature estimation depends on modeling parameters to a great extent, some modeling parameters are easy to be measured by experiments or obtained by theoretical derivation, while some parameters are difficult to obtain accurate values directly, and internal thermal resistance parameters are one of the parameters. In the soft-package battery, the heat is realized from the battery cell to the soft package through the conduction, and the thermal resistance exists in the gas-liquid mixture of electrolyte and reaction gas in the clearance between the battery cell and the soft package. The most important thermal resistance parameter is the thermal conductivity of the gas-liquid mixture. The prior research shows that the heat conductivity of the gas-liquid mixture is related to both temperature and pressure, and the steady-state heat conductivity is controversial in the aspect of a theoretical calculation method. Under the condition of thermal abuse, a dynamic process of drastic change of temperature and pressure exists in the soft package battery, and the dynamic thermal conductivity of the gas-liquid mixture is more difficult to calculate by adopting a theory. On the other hand, the method is limited by the packaging form and the testing means of the soft package battery, and the dynamic curve of the thermal conductivity of the gas-liquid mixture obtained by the testing method is also challenging.
Therefore, aiming at the difficult problem of estimating the internal thermal resistance parameter of the lithium ion soft package battery, the internal thermal resistance parameter estimation method of the lithium ion soft package battery, which has the advantages of wide applicability, high precision and good usability, is developed, and the internal thermal resistance parameter testing device of the soft package battery, which has a simple structure and is easy to implement, is designed, so that the temperature prediction precision of the lithium ion soft package battery under the condition of thermal abuse can be effectively enhanced, and the method has important engineering significance for improving the reliability and the thermal management level of the lithium ion soft package battery.
Disclosure of Invention
The invention overcomes the defects of the prior art and provides a method and a device for reversely solving the internal thermal resistance parameter of a lithium ion soft package battery. Firstly, defining internal thermal resistance parameters of a lithium ion soft package battery to be solved, and constructing a thermal resistance parameter reverse solving experimental device to implement a thermal induction experiment; secondly, physical equations such as heat generation, gas generation, heat transfer, expansion and the like are established and fused into a soft package battery temperature evaluation model; and thirdly, constructing a reverse model of the internal thermal resistance parameter of the soft package battery, and solving based on the test data and the physical model to obtain the internal thermal resistance parameter of the battery. The invention provides a key parameter estimation tool for accurately calculating the temperature of the lithium ion soft package battery.
In order to achieve the purpose, the invention adopts the following technical scheme:
(1) defining internal thermal resistance parameters of the lithium ion soft package battery to be solved;
(2) constructing a reverse experimental device of thermal resistance parameters;
(3) implementing a soft package battery thermal induction experiment;
(4) establishing a soft package battery temperature evaluation model;
(5) and constructing a reverse model of the internal thermal resistance parameter of the soft package battery and solving to obtain the thermal resistance parameter.
Further, in the step (1), the internal thermal resistance parameter of the lithium-ion soft package battery to be obtained comprises initial thermal conductivityK 0Thermal conductivity temperature indexn(ii) a Initial thermal conductivityK 0The thermal conductivity of a gas-liquid mixture filled in a gap between a battery cell and a soft package in the soft package battery at the temperature of 273.15K is shown; temperature index of thermal conductivitynIs expressed asn=log(K w /K 0)/log(T cel /273.15), whereinK w Indicating the temperature of the gas-liquid mixtureT cel Thermal conductivity of the steel.
Further, in the step (2), the thermal resistance parameter reverse-solving experimental device comprises a soft package battery, a clamping module, a temperature change test box, a temperature measuring instrument and a coordinate system; the external dimension of the soft package batteryX cel ×Y cel ×Z cel Initial volume ofV 0The battery consists of a battery cell, a soft package and electrolyte; the soft package is thicknessd p The aluminum-plastic film of (1); the clamping module comprises a first clamping plate, a second clamping plate, a first bolt, a second bolt, a first nut and a second nut; the first clamping plate and the second clamping plate are completely consistent in structure, the epoxy glass cloth laminated plate is selected as a material, and the overall dimension isX b ×Y b ×Z b Is composed ofX b = 0.15·X cel ，Y b = Y cel +0.3·X cel ，Z b = 0.1·Y cel In aX b Directed away from the first endD b = 0.07·X cel At a distance from the second endD b Are respectively provided with diameters ofD b A first circular hole and a second circular hole; the first bolt and the second bolt are inner hexagonal flat round head screws with standard numbers of GB/T70.2-2015, and the pitch diameters of the threadsd b =int(D b ) Length of threadh b =int(Z cel +2·Z b +1.5·D b ) Where int (·) denotes an integer arithmetic; the first nut and the second nut are butterfly nuts with standard numbers of GB62-88, and the pitch diameter of the threads isd b . The temperature change test chamber is a general test device which can provide a rapid temperature change environment and a constant temperature environment; the thermometer is a common thermometer with at least two temperature measuring channelsThe thermocouple type thermodetector comprises a first thermocouple, a second thermocouple, a signal wire and a host.
Further, in the step (2), the assembling steps of the thermal resistance parameter reverse experiment device are as follows:
(2.1) aligning and sequentially stacking the centers of the first clamping plate, the thermocouple, the soft package battery and the second clamping plate;
(2.2) passing the first bolt through the first round holes of the first clamping plate and the second clamping plate, and fastening the first bolt with the first nut;
(2.3) passing the second bolt through the second round holes of the first clamping plate and the second clamping plate, and fastening the second bolt with the second nut;
(2.4) controlling the distance between the first and second clamping plates to be withinZ bb = 0.95·Z cel ；
(2.5) adhering the second thermocouple to the upper surface of the soft package battery;
(2.6) suspending the assembled soft package battery, the clamping module and the temperature measuring instrument at the central position inside the temperature change test box;
(2.7) placing the host of the temperature measuring instrument at a position which is convenient to observe outside the temperature change test box.
Further, in the step (3), the heat induction experiment of the soft package battery is implemented as follows:
(3.1) based on the thermal resistance parameter reverse experiment device, setting the temperature change process in the temperature change test box as follows: the initial temperature is 298K, the temperature is kept for 20min, the temperature is increased to 398K at the speed of 5K/min, and the temperature is kept for 20 min.
(3.2) reading the measured temperatures of the battery cell measured by the first thermocouple and the second thermocouple in the experimental processT’ cel And the measured temperature of the soft packT’ sur Reading the temperature of the air in the box output by the temperature change test boxT air Recording the current time every 1 minutetIs/are as followsT’ cel ，T’ sur ，T air (ii) a A series ofT’ cel ，T’ sur ，T air Can be written as vectors:
T’ cel = (T’ cel,1 , T’ cel,2, ... , T’ cel,t , ... , T’ cel,T )
T’ sur = (T’ sur,1 , T’ sur,2, ... , T’ sur,t , ... , T’ sur,T )
T air = (T air,1 , T air,2, ... , T air,t , ... , T air,T )
whereinTIndicating the time of the last recording.
Further, in the step (4), the step of establishing the soft package battery temperature evaluation model is as follows:
(4.1) establish the heat generation equation E1:q cel =m rea ·sum(Δh j ·R j )， j=1, 2; wherein sum (·) represents a summation operation;q cel which indicates the rate of heat generation,m gas which represents the mass of the released gas,R 1, R 2respectively represent the reaction rates of SEI decomposition and reconstruction,Δh 1, Δh 2respectively representing reaction enthalpy changes of SEI decomposition and reconstruction;
further, in step (4.1), the reaction rate of the SEI film was calculated using the lithium ion battery thermal model proposed by Hatchard and DahnR j Differential equation of (a):R j =-dt(x j )=f j ·x j ·exp(-E j /(R u ·T cel ) ); wherein dt (-) denotes timetThe derivative operation of (1), i.e., the rate of change; exp (·) represents an exponential operation,x 1, x 2respectively represent the amount of lithium ions contained and inserted in the SEI film;f 1, f 2respectively representing frequency factors of SEI decomposition and reconstruction reaction,E 1, E 2respectively representing activation energy of SEI decomposition and reconstruction reaction;R u is the universal gas constant;T cel representing a temperature of the cell; solving for the reaction rateR j The differential equation of (c) yields:
R j = f j ·x j 0·exp{ -E j /(R u ·T cel ) - f j ·t·exp[ -E j /(R u ·T cel )]}
wherein the content of the first and second substances,tis the current time;x 1 0, x 2 0respectively indicate the initial amount of lithium ions contained in and inserted into the SEI film.
(4.2) establishing the gas production equation E2: dt (m gas )=c gas ·q cel (ii) a Wherein the content of the first and second substances,m gas which represents the mass of the released gas,c gas the expression indicates a proportionality coefficient between a gas generation rate and a heat generation rate.
(4.3) based on the theory of heat transfer, the heat transfer equation E3 is established:
q cel = C cel ·dt(T cel )+(T cel - T sur )/R cel
wherein the content of the first and second substances,C cel the heat capacity of the cell is represented,R cel representing the equivalent thermal conduction resistance from the cell to the pouch,T sur indicating the temperature of the soft pack.
(4.4) based on the thermal convection theory, the thermal convection equation E4 is established:
(T cel - T sur )/R cel = (T sur - T air )/R air
wherein the content of the first and second substances,R air representing the convective thermal resistance of the air inside the temperature-dependent test chamber from the outer surface of the soft pack.
(4.5) establishing an equivalent thermal resistance equation E5 based on the parallel thermal circuit theory.
(4.6) establishing a soft package deformation equation E6 based on finite element analysis;
(4.7) establishing an equivalent clearance equation E7:w=0.5·V 0·C e (P gas )/(X cel ·Y cel )。
(4.8) equations E1-E7 are combined to obtain a soft package battery temperature evaluation model M1:
further, in step (4.5), the step of establishing the equivalent thermal resistance equation E5 is as follows:
(4.5.1) establishing a parallel thermal resistance equation:R cel =1/sum(K w ·A i /w i ), i1, 2., 6; wherein the content of the first and second substances,A i is shown asiA plurality of heat transfer surfaces, wherein the soft package battery has 6 heat transfer surfaces in 3 directions such as X, Y, Z,w i is to show toiA gap between the heat transfer surface and the soft bag,K w in relation to the cell temperature, write:K w =K 0·(T cel /273.15) n ；
(4.5.2) introducing an equivalent gapw=mean(w i ) Mean (-) denotes the averaging operation;
(4.5.3) based on the equivalent gapwThe parallel thermal resistance equation is rewritten into an equivalent thermal resistance equation E5:R cel = w·A sur -1·K 0 -1·(0.0127·T cel ) n-(ii) a Wherein the content of the first and second substances,A sur =sum(A i ) The surface area of the soft package battery is shown.
Further, in step (4.6), the step of establishing the soft package deformation equation E6 is as follows:
(4.6.1) establishing a reaction gas equation of state according to the ideal gas law:P cel ·V gas =m gas ·R gas ·T cel ；
wherein the content of the first and second substances,P cel represents the internal gas pressure of the pouch battery,V gas which represents the volume of the gas,R gas represents a gas constant of the reaction gas;
(4.6.2) based on the thermal resistance parameter reverse experiment device, establishing 1/4 type finite element structures including a second soft package through the existing commercial finite element analysis software;
(4.6.3) cell arrangement, the cells of the second soft pack being arranged to have a thickness ofd p The housing unit of (a);
(4.6.4) setting material properties including: elastic modulus 4106 MP, Poisson ratio 0.33, stress and plastic strain table;
(4.6.5) setting boundary conditions including: setting an X-direction symmetric boundary for the first area, a Y-direction symmetric boundary for the second area, a Z-direction symmetric boundary for the third area, and a clamped boundary for simulating the clamping module for the fourth area;
(4.6.6) providing a fifth area of the second soft packmA different pressure loadP cel,j (j=1,2,...,m) Establishing a corresponding solving condition;
(4.6.7) calling a general static solver of the commercial finite element analysis software, solving the solving conditions one by one to obtain gas volume responseV gas,j (j=1,2,...,m) Calculating the relative expansion coefficientC e, j =V gas,j /V 0；
(4.6.8) based on themGroup pressure loadP cel, j And the gas volume responseV gas, j Fitting to obtain the relationship between the pressure load and the gas volume as follows:V gas =V gas (P cel ) Substituting into the relation between pressure load and internal air pressureP cel =P gas -P 0Obtaining a pressure-volume equation:V gas =V gas (P gas ) WhereinP 0Represents standard atmospheric pressure;
(4.6.9) substituting the gas pressure-volume equation into the gas equation of state to obtain the soft pack deformation equation E6:P cel ·V gas (P gas ) = m gas ·R gas ·T cel 。
further, in the step (5), a reverse model of the internal thermal resistance parameter of the soft package battery is constructed, and the thermal resistance parameter is obtained by solving the model as follows:
(5.1) establishing a cell temperature function and a soft package temperature function, wherein the method comprises the following steps: at the initial thermal conductivityK 0And the thermal conductivity temperature indexnCalling a solver of the conventional commercial numerical analysis software to solve the temperature evaluation model for input variables to obtain the cell temperature and the soft package temperature at all moments as output variables, and writing the output variables into vectors
T cel = (T cel,1 , T cel,2, ... , T cel,t , ... , T cel,T )
T sur = (T sur,1 , T sur,2, ... , T sur,t , ... , T sur,T )
WhereinT cel,t =T cel (K 0,n,t) As a function of the cell temperature function,T sur,t =T sur (K 0,n,t) Is the soft package temperature function.
(5.2) establishing a reverse optimization objective function, and writing into:F obj (K 0,n)=||T cel -T’ cel ||+||T sur -T’ sur and | l, wherein | · | | represents vector modulo arithmetic.
(5.3) establishing a reverse model M2 of the internal thermal resistance parameter of the soft package battery, and writing the model into a standard optimized model form:
min F obj (K 0, n)
s.t. K 0>0, n>0
where min represents the minimization and s.t. represents the constraint.
(5.4) calling a solver of the commercial numerical analysis software to solve the thermal resistance parameter inverse model to obtain the optimal initial thermal conductivityK 0 * And an optimum thermal conductivity temperature indexn * Namely, the solution of the thermal resistance parameter to be solved.
Compared with the prior art, the invention has the advantages that:
firstly, in the aspect of applicability, the thermal resistance parameter obtained by the method is suitable for estimating the internal temperature of the lithium ion soft package battery with the same electrolyte, and has a larger application range. In the aspect of precision, the method is based on experimental data, combines physical models of heat generation, gas generation, heat transfer, expansion and the like, has theoretical guarantee and experimental basis on the precision of thermal resistance parameters obtained by a reverse calculation technology, and has more advantages compared with solutions of theoretical derivation methods or empirical estimation methods. Thirdly, the method does not relate to complicated and tedious coupled simulation of multiple physical fields and algorithm programming processes, has clear steps, is beneficial to understanding, and has excellent usability for general engineering technicians. Finally, the device for reversely solving the thermal resistance parameter in the lithium ion soft package battery has the advantages of simple structure, clear definition, no need of manufacturing complex tools and high-end instruments, and good engineering practicability.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention.
FIG. 2 is a schematic diagram of a reverse experimental apparatus for thermal resistance parameters in an embodiment of the present invention.
Fig. 3 is a schematic diagram of expansion finite element analysis of a pouch battery in an embodiment of the present invention.
Reference numerals: 20. a thermal resistance parameter reverse experiment device; 21. a pouch cell; 22. a clamping module; 23. a temperature change test chamber; 24. a temperature measuring instrument; 211. an electric core; 212. soft packing; 213. an electrolyte; 221. a first splint; 222. a second splint; 223. a first bolt; 224. a second bolt; 225. a first nut; 226. a second nut; 241. a first thermocouple; 242. a second thermocouple; 243. a signal line; 244. a host; 2211. a first end; 2212. a second end; 2213. a first circular hole; 2214. a second circular hole; 30. a finite element structure; 31. a second soft package; 32. a first region; 33. a second region; 34. a third region; 35. a fourth region; 36. and a fifth region.
Detailed Description
The invention will now be further described with reference to the following examples, which are not to be construed as limiting the invention in any way, and any limited number of modifications which can be made within the scope of the claims of the invention are still within the scope of the claims of the invention.
As shown in fig. 1-3, the invention provides a reverse method and device for obtaining internal thermal resistance parameters of a lithium ion soft package battery, wherein the method comprises the following processing steps:
step S1: and defining the internal thermal resistance parameter of the lithium ion soft package battery to be solved. The soft package battery consists of an electric core, a soft package and Electrolyte, wherein a Solid Electrolyte Interface (SEI) on the surface of an anode of the electric core generates decomposition-reconstruction reaction under a certain temperature condition, and heat and gas are released. The reaction gas and the electrolyte form a gas-liquid mixture which is filled in the gap between the battery cell and the soft package. The heat of reaction has to be conducted through the gap to the soft pack and finally released to the environment. Thermal conductivity of gas-liquid mixture and its initial thermal conductivity at 273.15KK 0Thermal conductivity temperature indexnAnd the like. Therefore, the internal thermal resistance parameter of the lithium ion soft package battery to be obtained is defined as: initial thermal conductivityK 0Thermal conductivity temperature indexn。
Step S2: and constructing a reverse experimental device of thermal resistance parameters. As shown in fig. 2, in the embodiment, the thermal resistance parameter inverse experiment apparatus 20 includes a pouch cell 21, a clamping module 22, a temperature change test box 23, a temperature measuring instrument 24, and a coordinate system 25. External dimension of soft package battery 21X cel ×Y cel ×Z cel 82mm by 63mm by 4.3mm, initial volumeV 0=22214mm3The weight is 20 g, the capacity is 3940 mAh, and the battery comprises a battery cell 211, a soft package 212 and electrolyte 213. The anode material of the battery core 211 is graphite, and the cathode material is lithium cobaltate; the soft bag 212 is an aluminum plastic film (thickness)d p =0.088 mm); the electrolyte 213 is prepared by dissolving lithium hexafluorophosphate in a solvent of ethylene carbonate and diethyl carbonate with a mass ratio of 1:1, wherein the mole fraction of lithium hexafluorophosphate is 1.0 mol/L. The clamping module 22 includes a first clamping plate 221, a second clamping plate 222, a first bolt 223, a second bolt 224, a first nut 225, and a second nut 226. The first clamping plate 221 and the second clamping plate 222 have the same structure, are made of epoxy glass cloth laminated plates and have the same overall dimensionX b ×Y b ×Z b Is composed ofX b = 0.15·X cel ，Y b = Y cel +0.3·X cel ，Z b = 0.1·Y cel In aX b Oriented away from the first end 2211D b = 0.07·X cel From the second end 2212D b Are respectively provided with diameters ofD b First circular aperture 2213, second circular aperture 2214. The first bolt 223 and the second bolt 224 are inner hexagonal flat round head screws with standard numbers GB/T70.2-2015 and the middle diameters of threadsd b =int(D b ) Length of threadh b =int(Z cel +2·Z b +1.5·D b ) Where int (·) denotes an integer arithmetic. The first nut 225 and the second nut 226 are wing nuts with standard number GB62-88, and the pitch diameter of the threads isd b . The temperature change test chamber 23 is a model XX temperature change test apparatus manufactured by the company XX. The thermometer 24 is a type XX thermocouple thermometer manufactured by XX corporation, and includes a first thermocouple 241, a second thermocouple 242, a signal line 243, and a host computer 244. The assembly process of the thermal resistance parameter reverse experiment device 20 is as follows: the first clamping plate 221, the thermocouple 241, the soft package battery 21 and the second clamping plate 222 are aligned in center and stacked in sequence; secondly, the first bolt 223 passes through the first circular hole 2213 of the first clamping plate 221 and the second clamping plate 222 to be coupled and fastened with the first nut 225, the second bolt 224 passes through the second circular hole 2214 of the first clamping plate 221 and the second clamping plate 222 to be coupled and fastened with the second nut 226, and the distance between the first clamping plate 221 and the second clamping plate 222 is controlled to be betweenZ bb = 0.95·Z cel (ii) a Sticking a second thermocouple on the upper surface of the soft package battery 21; thirdly, the assembled soft package battery 21, the clamping module 22 and the temperature measuring instrument 24 are suspended in the central position inside the temperature change test box 23, and the host 243 of the temperature measuring instrument 24 is placed outside the temperature change test box 23 to facilitate observation.
Step S3: a pouch cell heat induction experiment was performed. The experimental procedure was as follows:
(3.1) reversely solving the experimental device 20 based on the thermal resistance parameters, and setting the temperature change process in the temperature change test box as follows: the initial temperature is 298K, the temperature is kept for 20min, the temperature is increased to 398K at the speed of 5K/min, and the temperature is kept for 20 min.
(3.2) during the experiment, the measured temperature of the battery cell 211 measured by the first thermocouple 241 and the second thermocouple 242 is readT’ cel And the measured temperature of the soft pack 212T’ sur Reading the temperature of the air in the chamber output from the temperature change test chamber 23T air Recording the current time every 1 minutetIs/are as followsT’ cel ，T’ sur ，T air . A series ofT’ cel ，T’ sur ，T air Can be written as vectors:
T’ cel = (T’ cel,1 , T’ cel,2, ... , T’ cel,t , ... , T’ cel,T )
T’ sur = (T’ sur,1 , T’ sur,2, ... , T’ sur,t , ... , T’ sur,T )
T air = (T air,1 , T air,2, ... , T air,t , ... , T air,T )
whereinTIndicating the time of the last recording. Some of the test data are listed in table 1.
TABLE 1
Step S4: and establishing a soft package battery temperature evaluation model M1. The modeling steps are as follows:
(4.1) establish the heat generation equation E1:q cel = m rea ·sum(Δh j ·R j )， j = 1, 2
wherein sum (·) represents a summation operation;q cel which indicates the rate of heat generation,m gas which represents the mass of the released gas,R 1, R 2respectively represent the reaction rates of SEI decomposition and reconstruction,Δh 1, Δh 2respectively show the reaction enthalpy changes of SEI decomposition and reconstruction. Calculation of the reaction rate of SEI films Using the lithium ion Battery thermal model proposed by Hatchard and DahnR j ：
R j = - dt(x j ) = f j ·x j ·exp( -E j /(R u ·T cel ) )
Wherein dt (-) denotes timetExp (-) denotes an exponential operation,x 1, x 2respectively represent the amount of lithium ions contained and inserted in the SEI film;f 1, f 2respectively representing frequency factors of SEI decomposition and reconstruction reaction,E 1, E 2respectively representing activation energy of SEI decomposition and reconstruction reaction;R u is the universal gas constant;T cel indicating the temperature of the cell 211. The above equation is a differential equation, which is solved to obtain:
R j = f j ·x j 0·exp{ -E j /(R u ·T cel ) - f j ·t·exp[ -E j /(R u ·T cel )]}
wherein the content of the first and second substances,tis the current time;x 1 0, x 2 0respectively indicate the initial amount of lithium ions contained in and inserted into the SEI film.
(4.2) establishing the gas production equation E2: dt (m gas )= c gas ·q cel
Wherein the content of the first and second substances,m gas denotes the mass of released gas, dt: (m gas ) The gas production rate is shown,c gas the expression indicates a proportionality coefficient between a gas generation rate and a heat generation rate.
(4.3) based on the theory of heat transfer, the heat transfer equation E3 is established:
q cel = C cel ·dt(T cel )+(T cel - T sur )/R cel
wherein the content of the first and second substances,C cel the thermal capacity of the cell 211 is shown,R cel representing the equivalent thermal conduction resistance from the cell 211 to the soft pack 212,T sur indicating the temperature of the soft pack 212.
(4.4) based on the thermal convection theory, the thermal convection equation E4 is established:
(T cel - T sur )/R cel = (T sur - T air )/R air
wherein the content of the first and second substances,R air showing the convective resistance of the air inside the experimental box 23 as a function of temperature from the outer surface of the soft pack 212.
(4.5) establishing an equivalent thermal resistance equation E5 based on the parallel thermal circuit theory, and the steps are as follows.
(4.5.1) establishing a parallel thermal resistance equation:
R cel = 1 / sum( K w ·A i /w i ), i = 1, 2, ... , 6
wherein the content of the first and second substances,A i is shown asiThe soft package battery 21 has 6 heat transfer surfaces in 3 directions such as X, Y, Z,w i is to show toiThe gap between the heat transfer surface and the soft bag 212 from the battery cell 211,K w in relation to the cell 211 temperature, write:K w =K 0·(T cel /273.15) n ；
(4.5.2) introducing an equivalent gapw=mean(w i ) Mean (-) indicates the averaging operation.
(4.5.3) based on equivalent gapwThe parallel thermal resistance equation is rewritten into an equivalent thermal resistance equation E5:
R cel =w·A sur -1·K 0 -1·(0.0127·T cel ) n-
wherein the content of the first and second substances,A sur =sum(A i ) Representing the surface area of pouch cell 21.
(4.6) based on finite element analysis, the soft pack deformation equation E6 is established, as follows.
(4.6.1) establishing a reaction gas equation of state according to the ideal gas law:
P cel ·V gas = m gas ·R gas ·T cel
wherein the content of the first and second substances,P cel indicating the internal gas pressure of the pouch cell 21,V gas which represents the volume of the gas,R gas represents a gas constant of the reaction gas;
(4.6.2) based on the thermal resistance parameter back-solving experiment device 20, establishing 1/4 type finite element structure 30 in the finite element analysis software ABAQUS, wherein the finite element structure comprises a second soft package 31;
(4.6.3) setting the thickness of the second soft pack 31 to be a thicknessw p The shell element of (a), corresponding to the soft pack 212 in fig. 2;
(4.6.4) setting material properties including: modulus of elasticity 4106 MP, poisson's ratio 0.33, stress and plastic strain table (as in table 2);
TABLE 2
(4.6.5) setting boundary conditions including: setting an X-direction symmetric boundary for the first region 32, a Y-direction symmetric boundary for the second region 33, a Z-direction symmetric boundary for the third region 34, and a clamped boundary for simulating the clamping module 22 shown in FIG. 2 for the fourth region 35;
(4.6.6) providing a fifth area 36 of the second soft pack 31m=6 different pressure loadsP cel,j (j=1,2,...,m) Establishing a corresponding solving condition;
(4.6.7) calling a general static solver of ABAQUS to respectively solve each solution condition to obtain gas volume responseV gas,j (j=1,2,...,m) The calculation results are shown in table 3;
TABLE 3
(4.6.8) based onmGroup pressure loadP cel, j And gas volume responseV gas,j Fitting to obtain the relationship between the pressure load and the gas volume as follows:V gas =0.658·V 0·P cel substituting into the relation between pressure load and internal air pressureP cel =P gas -P 0Obtaining a pressure-volume equation:V gas =0.658·V 0·(P gas -P 0)，P 0indicating standard atmospheric pressure.
(4.6.9) substituting the pressure-volume equation into the gas equation of state to obtain the deformation equation for soft pack E6:
0.658·P cel ·V 0·(P gas -P 0) = m gas ·R gas ·T cel
(4.7) establishing an equivalent clearance equation E7:w=0.5·V gas /(X cel ·Y cel )。
(4.8) equations E1-E7 are combined to obtain a soft package battery temperature evaluation model M1:
in this example, the parameters of the temperature evaluation model M1 are shown in table 4.
TABLE 4
Step S5: constructing a reverse model M2 of the internal thermal resistance parameter of the soft package battery and solving to obtain the thermal resistance parameter, wherein the modeling steps are as follows:
and (5.1) establishing a cell temperature function and a soft package temperature function. At an initial thermal conductivityK 0And temperature index of thermal conductivitynCalling an ODE solver in MATLAB (matrix laboratory) of commercial numerical analysis software to solve the temperature evaluation model M1 to obtain the cell temperature and the soft package temperature at all moments as output variables, and writing the output variables into vectorsT cel =(T cel,1,T cel,2,...,T cel,t ,...,T cel,T ) AndT sur =(T sur,1,T sur,2,...,T sur,t ,...,T sur,T ) (ii) a WhereinT cel,t = T cel (K 0, n, t) For the temperature of the cellThe function of the degree function is that,T sur,t = T sur (K 0, n, t) As a function of the temperature of the soft pack.
(5.2) establishing a reverse optimization objective function, and writing into:F obj (K 0,n)=||T cel -T’ cel ||+||T sur -T’ sur and | l, wherein | · | | represents vector modulo arithmetic.
(5.3) establishing a reverse model M2 of the internal thermal resistance parameter of the soft package battery, and writing the model into a standard optimized model form:
min F obj (K 0, n)
s.t. K 0>0, n>0
where min represents the minimization and s.t. represents the constraint.
(5.4) calling an optimization problem solver 'fmincon' of MATLAB to solve the model M2 to obtain the optimal initial thermal conductivityK 0 * =0.055W/(m.K) and optimum thermal conductivity temperature indexn * =1.84, the solution of the thermal resistance parameter is required.
The embodiment shows that the method and the device for reversely solving the internal thermal resistance parameter of the lithium ion soft package battery have the following advantages. Firstly, in the aspect of applicability, the thermal resistance parameter obtained by the method is suitable for estimating the internal temperature of the lithium ion soft package battery with the same electrolyte, and has a larger application range. For example, the sample electrolyte adopted in the embodiment is widely applied to the lithium ion soft package battery, and the obtained thermal resistance parameter is suitable for a series of similar battery products, so that the development cost of a battery manufacturer can be greatly saved, and the product development iteration cycle can be shortened. In the aspect of precision, the method is based on experimental data, combines physical models such as heat generation, gas generation, heat transfer, expansion and the like (such as a soft-package battery temperature evaluation model M1 in step 4.8), has theoretical guarantee and experimental basis on the precision of thermal resistance parameters obtained by a reverse solving technology, and has more advantages compared with solutions of theoretical derivation methods or empirical estimation methods. Thirdly, the method does not relate to complicated and tedious coupled multi-physical field simulation and algorithm programming processes (for example, models M1 and M2 can be solved by MATLAB), and the steps are clear, easy to understand and easy to use for general engineering technicians. Finally, the reverse thermal resistance parameter solving device in the lithium ion soft package battery provided by the invention is simple in structure and clear in definition, does not need to manufacture complex tools, does not need high-end instruments (such as adiabatic acceleration calorimeter/ARC), and has good engineering practicability.
Claims (5)
1. The method for reversely calculating the internal thermal resistance parameter of the lithium ion soft package battery is characterized by comprising the following processing steps of:
(1) defining internal thermal resistance parameters of the lithium ion soft package battery to be solved;
(2) constructing a reverse experimental device of thermal resistance parameters;
(3) implementing a soft package battery thermal induction experiment;
(4) establishing a soft package battery temperature evaluation model;
(5) constructing a reverse model of the internal thermal resistance parameters of the soft package battery and solving to obtain thermal resistance parameters;
in the step (1), the internal thermal resistance parameter of the lithium ion soft package battery to be solved comprises the following steps: initial thermal conductivityK 0Thermal conductivity temperature indexnThe initial thermal conductivityK 0The thermal conductivity of a gas-liquid mixture filled in a soft package battery between a battery core and a soft package gap is 273.15K, and the temperature index of the thermal conductivity isnIs expressed asn=log(K w /K 0)/log(T cel /273.15), whereinK w Indicating the temperature of the gas-liquid mixtureT cel Thermal conductivity of the steel.
2. The reverse method for the internal thermal resistance parameter of the lithium ion soft package battery according to claim 1, in the step (2), the reverse experimental device for the thermal resistance parameter comprises: the device comprises a soft package battery, a clamping module, a temperature change test box, a temperature measuring instrument and a coordinate system; what is needed isThe external dimension of the soft package batteryX cel ×Y cel ×Z cel Initial volume ofV 0The battery consists of a battery cell, a soft package and electrolyte; the soft package is thicknessd p The aluminum-plastic film of (1); the clamping module comprises a first clamping plate, a second clamping plate, a first bolt, a second bolt, a first nut and a second nut; the first clamping plate and the second clamping plate are completely consistent in structure, the epoxy glass cloth laminated plate is selected as a material, and the overall dimension isX b ×Y b ×Z b Is composed ofX b =0.15·X cel ，Y b =Y cel +0.3·X cel ，Z b =0.1·Y cel (ii) a In the above-mentionedX b Directed away from the first endD b =0.07·X cel At a distance from the second endD b Are respectively provided with diameters ofD b A first circular hole and a second circular hole; the first bolt and the second bolt are inner hexagonal flat round head screws with standard numbers of GB/T70.2-2015, and the pitch diameters of the threadsd b =int(D b ) Length of threadh b =int(Z cel +2·Z b +1.5·D b ) Where int (·) denotes an integer arithmetic; the first nut and the second nut are butterfly nuts with standard numbers of GB62-88, and the pitch diameter of the threads isd b (ii) a The temperature change test chamber is a general test device which can provide a rapid temperature change environment and a constant temperature environment; the temperature measuring instrument is a general thermocouple type temperature measuring instrument with at least two temperature measuring channels and comprises a first thermocouple, a second thermocouple, a signal wire and a host;
in the step (2), the assembling steps of the thermal resistance parameter reverse experiment device are as follows:
(2.1) aligning and sequentially stacking the centers of the first clamping plate, the thermocouple, the soft package battery and the second clamping plate;
(2.2) passing the first bolt through the first round holes of the first clamping plate and the second clamping plate, and fastening the first bolt with the first nut;
(2.3) passing the second bolt through the second round holes of the first clamping plate and the second clamping plate, and fastening the second bolt with the second nut;
(2.4) controlling the distance between the first and second clamping plates to be withinZ bb =0.95·Z cel ；
(2.5) adhering the second thermocouple to the upper surface of the soft package battery;
(2.6) suspending the assembled soft package battery, the clamping module and the temperature measuring instrument at the central position inside the temperature change test box;
(2.7) placing the host of the temperature measuring instrument at a position which is convenient to observe outside the temperature change test box.
3. The method for reversely calculating the internal thermal resistance parameter of the lithium-ion soft package battery according to claim 2, in the step (3), the step of implementing the soft package battery thermal induction experiment is as follows:
(3.1) based on the thermal resistance parameter reverse experiment device, setting the temperature change process of the temperature change test box as follows: keeping the initial temperature at 298K for 20min, heating to 398K at the rate of 5K/min, and keeping the temperature for 20 min;
(3.2) reading the measured temperatures of the battery cell measured by the first thermocouple and the second thermocouple in the experimental processT’ cel And the measured temperature of the soft packT’ sur Reading the temperature of the air in the box output by the temperature change test boxT air Recording the current time every 1 minutetIs/are as followsT’ cel ，T’ sur ，T air (ii) a A series ofT’ cel ，T’ sur ，T air Can be written as vectors:
T’ cel = (T’ cel,1, T’ cel,2, ... , T’ cel,t , ... , T’ cel,T )
T’ sur = (T’ sur,1, T’ sur,2, ... , T’ sur,t , ... , T’ sur,T )
T air = (T air,1 , T air,2, ... , T air,t , ... , T air,T )
whereinTIndicating the time of the last recording.
4. The method for reversely solving the internal thermal resistance parameter of the lithium-ion soft package battery according to claim 2 is characterized in that in the step (4), the step of establishing the soft package battery temperature evaluation model is as follows:
(4.1) establish the heat generation equation E1:q cel = m rea ·sum(Δh j ·R j )， j=1, 2; wherein the content of the first and second substances,j=1 represents an SEI decomposition reaction,j=2 represents SEI reconstruction reaction; sum (·) represents a summation operation;q cel which indicates the rate of heat generation,m gas which represents the mass of the released gas,R j respectively represent the reaction rates of SEI decomposition and reconstruction,Δh j representing the reaction enthalpy change of SEI decomposition and reconstruction;
(4.2) establishing the gas production equation E2: dt (m gas )=c gas ·q cel (ii) a Wherein the content of the first and second substances,m gas which represents the mass of the released gas,c gas representation ofA proportionality coefficient between a gas production rate and a heat generation rate;
(4.3) based on the theory of heat transfer, the heat transfer equation E3 is established:
q cel = C cel ·dt(T cel )+(T cel - T sur )/R cel
wherein the content of the first and second substances,C cel the heat capacity of the cell is represented,R cel representing the equivalent thermal conduction resistance from the cell to the pouch,T sur representing the temperature of the soft pack;
(4.4) based on the thermal convection theory, the thermal convection equation E4 is established:
(T cel - T sur )/R cel = (T sur - T air )/R air
wherein the content of the first and second substances,R air representing the convective thermal resistance of the air inside the temperature-dependent test box from the outer surface of the soft pack;
(4.5) establishing an equivalent thermal resistance equation E5 based on a parallel thermal circuit theory;
(4.6) establishing a soft package deformation equation E6 based on finite element analysis;
(4.7) establishing an equivalent clearance equation E7:w=0.5·C e (P gas )/(X cel ·Y cel )；
(4.8) equations E1-E7 are combined to obtain a soft package battery temperature evaluation model M1:
in the step (4.1), the thermal model of the lithium ion battery proposed by Hatchard and Dahn is adopted to calculate the SEI decomposition and reconstruction reaction rate：R j =-dt(x j )=f j ·x j ·exp(-E j /(R u ·T cel ) ); wherein dt (-) denotes timetThe derivative operation of (1), i.e., the rate of change; exp (·) represents an exponential operation,x 1, x 2respectively represent the amount of lithium ions contained and inserted in the SEI film;f 1, f 2respectively representing frequency factors of SEI decomposition and reconstruction reaction,E 1, E 2respectively representing activation energy of SEI decomposition and reconstruction reaction;R u is the universal gas constant;T cel representing a temperature of the cell; solving the reaction rate of SEI decomposition and reconstruction to obtain:
R j = f j ·x j 0·exp{ -E j /(R u ·T cel ) - f j ·t·exp[ -E j /(R u ·T cel )]}
wherein the content of the first and second substances,tis the current time; in thatx j 0In (1),j=1 denotes the initial amount of lithium ions in the SEI film,j=2 represents an initial amount of lithium ion intercalation into the SEI film;
in step (4.5), the step of establishing the equivalent thermal resistance equation E5 is as follows:
(4.5.1) establishing a parallel thermal resistance equation:R cel =1/sum(K w ·A i /w i ), i1, 2., 6; wherein the content of the first and second substances,A i is shown asiA plurality of heat transfer surfaces, wherein the soft package battery has 6 heat transfer surfaces in 3 directions such as X, Y, Z,w i is to show toiA gap between the heat transfer surface and the soft bag,K w in relation to the cell temperature, write:K w =K 0·(T cel /273.15) n ；
(4.5.2) introducing an equivalent gapw= mean(w i ) Mean (-) denotes the averaging operation;
(4.5.3) based on the equivalent gapwThe parallel thermal resistance equation is rewritten into an equivalent thermal resistance equation E5:R cel =w·A sur -1·K 0 -1·(0.0127·T cel ) n-(ii) a Wherein the content of the first and second substances,A sur =sum(A i ) Representing the surface area of the soft package battery;
in step (4.6), the step of establishing the soft package deformation equation E6 is as follows:
(4.6.1) establishing a reaction gas equation of state according to the ideal gas law:P cel ·V gas =m gas ·R gas ·T cel ；
wherein the content of the first and second substances,P cel represents the internal gas pressure of the pouch battery,V gas which represents the volume of the gas,R gas represents a gas constant of the reaction gas;
(4.6.2) based on the thermal resistance parameter reverse experiment device, establishing 1/4 type finite element structures including a second soft package through the existing commercial finite element analysis software;
(4.6.3) cell arrangement, the cells of the second soft pack being arranged to have a thickness ofd p The housing unit of (a);
(4.6.4) setting material properties including: elastic modulus 4106 MP, Poisson ratio 0.33, stress and plastic strain table;
(4.6.5) setting boundary conditions including: setting an X-direction symmetric boundary for the first area, a Y-direction symmetric boundary for the second area, a Z-direction symmetric boundary for the third area, and a clamped boundary for simulating the clamping module for the fourth area;
(4.6.6) providing a fifth area of the second soft packmA different pressure loadP cel,j (j=1,2,...,m) Establishing a corresponding solving condition;
(4.6.7) calling a general static solver of the commercial finite element analysis software, solving the solving conditions one by one to obtain gas volume responseV gas,j (j=1,2,...,m) Calculating the relative expansion coefficientC e, j =V gas,j /V 0；
(4.6.8) based on themA different pressure loadP cel, j And the gas volume responseV gas, j Fitting to obtain the relationship between the pressure load and the gas volume as follows:V gas =V gas (P cel ) Substituting into the relation between pressure load and internal air pressureP cel =P gas -P 0Obtaining a pressure-volume equation:V gas =V gas (P gas ) WhereinP 0Represents standard atmospheric pressure;
(4.6.9) substituting the gas pressure-volume equation into the gas equation of state to obtain the soft pack deformation equation E6:P cel ·V gas (P gas )=m gas ·R gas ·T cel 。
5. the reverse calculation method for the internal thermal resistance parameter of the lithium ion soft package battery according to claim 1 is characterized in that in step (5), the steps of constructing a reverse calculation model for the internal thermal resistance parameter of the soft package battery and solving to obtain the thermal resistance parameter are as follows:
(5.1) establishing a cell temperature function and a soft package temperature function, wherein the method comprises the following steps: at the initial thermal conductivityK 0And the thermal conductivity temperature indexnCalling a solver of the conventional commercial numerical analysis software to solve the temperature evaluation model for input variables to obtain the cell temperature and the soft package temperature at all moments as output variables, and writing the output variables into vectors
T cel = (T cel,1 , T cel,2, ... , T cel,t , ... , T cel,T )
T sur = (T sur,1 , T sur,2, ... , T sur,t , ... , T sur,T )
WhereinT cel,t = T cel (K 0, n, t) As a function of the cell temperature function,T sur,t =T sur (K 0,n,t) Is the soft package temperature function;
(5.2) establishing a reverse optimization objective function, and writing into:F obj (K 0,n)=||T cel -T’ cel ||+||T sur -T’ sur | |, where | · | | | represents vector modulo arithmetic;
(5.3) establishing a reverse model M2 of the internal thermal resistance parameter of the soft package battery, and writing the model into a standard optimized model form:
min F obj (K 0, n)
s.t. K 0>0, n>0
where min represents the minimization and s.t. represents the constraint;
(5.4) calling a solver of the commercial numerical analysis software to solve the thermal resistance parameter inverse model to obtain the optimal initial thermal conductivityK 0 * And an optimum thermal conductivity temperature indexn * Namely, the solution of the thermal resistance parameter to be solved.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110486991.9A CN113189508B (en) | 2021-05-03 | 2021-05-03 | Reverse calculation method and device for internal thermal resistance parameters of lithium ion soft package battery |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110486991.9A CN113189508B (en) | 2021-05-03 | 2021-05-03 | Reverse calculation method and device for internal thermal resistance parameters of lithium ion soft package battery |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113189508A true CN113189508A (en) | 2021-07-30 |
| CN113189508B CN113189508B (en) | 2022-06-21 |
Family
ID=76983624
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110486991.9A Active CN113189508B (en) | 2021-05-03 | 2021-05-03 | Reverse calculation method and device for internal thermal resistance parameters of lithium ion soft package battery |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113189508B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106872898A (en) * | 2017-02-06 | 2017-06-20 | 中国第汽车股份有限公司 | Electrokinetic cell monomer interface thermal resistance method for rapidly testing |
| JP2018147680A (en) * | 2017-03-03 | 2018-09-20 | 住友電気工業株式会社 | Temperature abnormality determination device, temperature abnormality determination method, and computer program |
| CN109188285A (en) * | 2018-07-24 | 2019-01-11 | 合肥工业大学 | External thermal resistance estimation method in a kind of lithium battery |
| CN109613055A (en) * | 2018-12-27 | 2019-04-12 | 上海工程技术大学 | A kind of the stable state measuring method and measurement device of cylindrical battery radial direction thermal coefficient |
| CN109884527A (en) * | 2019-02-13 | 2019-06-14 | 深圳市比克动力电池有限公司 | A kind of lithium ion battery thermal runaway quantity of heat production calculation method |
| CN111044930A (en) * | 2019-12-06 | 2020-04-21 | 清华大学 | Method and system for testing heat transfer resistance of lithium ion battery module |
| CN111597719A (en) * | 2020-05-19 | 2020-08-28 | 华南理工大学 | Power battery simulation method based on electric heating and thermal runaway coupling model |
-
2021
- 2021-05-03 CN CN202110486991.9A patent/CN113189508B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106872898A (en) * | 2017-02-06 | 2017-06-20 | 中国第汽车股份有限公司 | Electrokinetic cell monomer interface thermal resistance method for rapidly testing |
| JP2018147680A (en) * | 2017-03-03 | 2018-09-20 | 住友電気工業株式会社 | Temperature abnormality determination device, temperature abnormality determination method, and computer program |
| CN109188285A (en) * | 2018-07-24 | 2019-01-11 | 合肥工业大学 | External thermal resistance estimation method in a kind of lithium battery |
| CN109613055A (en) * | 2018-12-27 | 2019-04-12 | 上海工程技术大学 | A kind of the stable state measuring method and measurement device of cylindrical battery radial direction thermal coefficient |
| CN109884527A (en) * | 2019-02-13 | 2019-06-14 | 深圳市比克动力电池有限公司 | A kind of lithium ion battery thermal runaway quantity of heat production calculation method |
| CN111044930A (en) * | 2019-12-06 | 2020-04-21 | 清华大学 | Method and system for testing heat transfer resistance of lithium ion battery module |
| CN111597719A (en) * | 2020-05-19 | 2020-08-28 | 华南理工大学 | Power battery simulation method based on electric heating and thermal runaway coupling model |
Non-Patent Citations (1)
| Title |
|---|
| 曹涛锋: "接触热阻对燃料电池温度分布影响的数值模拟", 《工程热物理学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113189508B (en) | 2022-06-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Integrated computation model of lithium-ion battery subject to nail penetration | |
| CN109900734B (en) | Method for identifying thermophysical parameters in square lithium ion battery | |
| Oh et al. | Characterization and modeling of the thermal mechanics of lithium-ion battery cells | |
| Xia et al. | Multiphysical modeling for life analysis of lithium-ion battery pack in electric vehicles | |
| Xie et al. | An improved resistance-based thermal model for prismatic lithium-ion battery charging | |
| Sheng et al. | Quasi steady state method to measure thermophysical parameters of cylindrical lithium ion batteries | |
| Xie et al. | A resistance-based electro-thermal coupled model for an air-cooled battery pack that considers branch current variation | |
| Xie et al. | Improving battery thermal behavior and consistency by optimizing structure and working parameter | |
| Cao et al. | A full-scale electrical-thermal-fluidic coupling model for li-ion battery energy storage systems | |
| Sheng et al. | Experimental and numerical approach for analyzing thermal behaviors of a prismatic hard-cased lithium-ion battery | |
| CN114692244A (en) | Lithium battery pack heat abuse safety risk assessment method based on multi-physical-field simulation | |
| CN114186437A (en) | Multi-physical-field coupling degradation model order reduction method for power system reliability simulation analysis | |
| CN113189508B (en) | Reverse calculation method and device for internal thermal resistance parameters of lithium ion soft package battery | |
| CN113835031B (en) | Information processing method, apparatus, electronic device and storage medium | |
| Zhai et al. | Thermal model and experimental verification of lithium-ion battery based on heat generation rate | |
| Yu et al. | A novel heat generation acquisition method of cylindrical battery based on core and surface temperature measurements | |
| Same et al. | Effect of thermal parameters on behaviour of A lithium-ion battery: simulation study | |
| Nguyen et al. | Deformation behavior of single prismatic battery cell cases and cell assemblies loaded by internal pressure | |
| Huang et al. | Simulation of thermal behavior evolution in lithium-ion pouch cells using a coupled reaction-thermal-expansion model | |
| CN113076518B (en) | Method for estimating internal temperature and pressure of single lithium ion soft package battery | |
| Yu et al. | Thermal parameters of cylindrical power batteries: Quasi-steady state heat guarding measurement and thermal management strategies | |
| Lyu et al. | Thermal Simulation for a 48V/50Ah Lithium-ion Battery Pack | |
| CN116451412B (en) | Soft package lithium battery multi-field coupling analysis model and solving method for edge calculation | |
| Shi et al. | Experimental determinations of thermophysical parameters for lithium-ion batteries: A systematic review | |
| Wu et al. | Influence of capacity and energy density of lithium-ion battery on thermal reaction during high rate charging |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210730 Assignee: Hunan Hongyuyun Technology Co.,Ltd. Assignor: Hunan City University Contract record no.: X2022980020750 Denomination of invention: A method and device for reverse calculation of internal thermal resistance parameters of lithium ion flexible battery Granted publication date: 20220621 License type: Exclusive License Record date: 20221109 |