CN113626946A - Method for determining boundary value of thickness and thermal conductivity of material for blocking battery thermal diffusion - Google Patents
Method for determining boundary value of thickness and thermal conductivity of material for blocking battery thermal diffusion Download PDFInfo
- Publication number
- CN113626946A CN113626946A CN202110778476.8A CN202110778476A CN113626946A CN 113626946 A CN113626946 A CN 113626946A CN 202110778476 A CN202110778476 A CN 202110778476A CN 113626946 A CN113626946 A CN 113626946A
- Authority
- CN
- China
- Prior art keywords
- thermal
- battery
- heat
- monomer
- thermal runaway
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000009792 diffusion process Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000000463 material Substances 0.000 title claims abstract description 31
- 230000000903 blocking effect Effects 0.000 title claims abstract description 23
- 239000000178 monomer Substances 0.000 claims abstract description 57
- 238000012360 testing method Methods 0.000 claims abstract description 35
- 238000009413 insulation Methods 0.000 claims abstract description 33
- 238000004088 simulation Methods 0.000 claims abstract description 25
- 239000012774 insulation material Substances 0.000 claims abstract description 21
- 230000008859 change Effects 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000013461 design Methods 0.000 claims description 22
- 238000010438 heat treatment Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 230000020169 heat generation Effects 0.000 claims description 8
- 230000007613 environmental effect Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000007707 calorimetry Methods 0.000 claims description 3
- 238000005485 electric heating Methods 0.000 claims description 3
- 238000004134 energy conservation Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 2
- 230000006872 improvement Effects 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 8
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 4
- 239000011810 insulating material Substances 0.000 description 4
- 229910001416 lithium ion Inorganic materials 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Secondary Cells (AREA)
Abstract
The embodiment of the invention provides a method for determining the thickness and the heat conductivity boundary value of a material for blocking thermal diffusion of a battery, which comprises the steps of obtaining material attribute parameters of battery monomers, thermal insulation thermal runaway temperature-time change test data of the battery monomers and temperature data of the central position of a contact surface between every two battery monomers and a thermal insulation plate in a module thermal diffusion process through tests, introducing a three-dimensional geometric model of a battery module into Ansys software, performing pretreatment and related setting, and obtaining the thickness and the heat conductivity boundary value of the thermal insulation plate material capable of blocking thermal diffusion through automatic batch simulation calculation. The embodiment of the invention can quickly determine the boundary value of the thickness and the thermal conductivity coefficient of the thermal insulation material capable of blocking thermal diffusion, determines the shape selection direction of the optimal thermal insulation material, saves a large amount of time and cost, is safer and more environment-friendly, can promote the improvement of the safe application technical level of the power battery, and has the condition of high engineering application.
Description
Technical Field
The invention relates to the technical field of batteries, in particular to a method for determining boundary values of thickness and thermal conductivity of a material for blocking battery thermal diffusion.
Background
Along with the rapid increase of the holding capacity of new energy automobiles and the popularization and application of new system high specific energy lithium ion batteries, the accidents of fire and even explosion of lithium ion power battery packs occur, and the attention and worry of consumers are more and more aroused. The safety problem of the lithium ion power battery can be finally attributed to thermal runaway of a battery cell and thermal diffusion of a battery system. The importance of effective protection against thermal runaway and diffusion of lithium ion batteries has been recognized by both government agencies, as well as by various large battery plants and whole vehicle plants.
At present, a heat insulation material with a certain thickness is arranged among different monomers in a battery pack, and the scheme is widely adopted in various large battery factories and whole automobile factories to delay or block battery heat diffusion. In order to determine a heat insulating material which can effectively block heat diffusion and has reasonable cost, various heat insulating materials are often selected by enterprises to perform corresponding heat diffusion tests on a battery system in the design and development process of the battery system. Since the solutions of any kind of material at different thicknesses need to be verified separately, the number of tests usually reaches dozens or even hundreds, which undoubtedly requires a lot of time and cost.
In addition, according to the above-described trial and error method, it is only known that one or some of the materials can block heat diffusion, and when the heat diffusion can not be blocked, the boundary value of a certain property parameter of the heat insulating material cannot be clarified. The thermal insulation material mainly relates to three material attribute parameters of thickness, density, specific heat and thermal conductivity, wherein the change of the thickness and the thermal conductivity plays a more critical role in the thermal insulation effect, so that the boundary value of the thickness and the thermal conductivity of the thermal insulation material capable of blocking thermal diffusion is quickly determined, and the subsequent further optimization of the thermal insulation material is facilitated.
Disclosure of Invention
The embodiment of the invention aims to overcome the defects that a large amount of time and money cost are consumed and the optimization and selection of a heat insulation material are not facilitated in the process of selecting the heat diffusion blocking material of the battery through a traditional repeated test mode, and the like, and provides an analysis method for rapidly determining the thickness and the heat conductivity boundary value of the heat insulation material for blocking the heat diffusion of the power battery.
The embodiment of the invention provides a method for rapidly determining the boundary value of the thickness and the thermal conductivity of a thermal insulation material for blocking the thermal diffusion of a battery, which comprises the following steps:
Step S100, obtaining material attribute parameters of single batteries, test data of thermal insulation thermal runaway temperature-time change of the single batteries and temperature data of the center position of a contact surface between every two single batteries and a thermal insulation plate through tests;
step S200, building a three-dimensional geometric model for the battery module, importing the three-dimensional geometric model into Ansys Fluent Meshing, setting grid parameters according to the three-dimensional geometric model, dividing a surface network, and then generating a body network;
step S300, switching to an Ansys Fluent solving mode, setting a physical model, fitting the thermal insulation thermal runaway temperature-time change test data of the battery monomer to obtain the relevant parameters of the thermal runaway reaction kinetic equation of the battery monomer, and endowing each battery monomer with a thermal runaway model; setting corresponding material attribute parameters for different areas of the volume grid; setting heat generation power and heat generation time for the heating plate; setting an environmental boundary condition for a surface in direct contact with an environment;
step S400, running simulation calculation, monitoring temperature values between every two battery monomers and at the center position of a contact surface of the heat insulation plate, comparing the temperature values with temperature data of different positions of the battery module obtained in the test, and verifying the effectiveness and accuracy of a simulation model;
Step S500, setting the thermal conductivity of a thermal insulation material as an Input Parameter in the Ansys Fluent based on the verified simulation model, and setting a Parameter alpha representing the thermal runaway reaction process of each monomer as an Output Parameter;
step S600, loading a Fluent version block in Ansys Workbench, importing the set simulation models of the Input and Output parameters, establishing a plurality of Design points according to needs, modifying the specific value of the Input Parameter in each Design Point, and selecting different heat conductivity values in a certain range corresponding to the heat insulation board material;
step S700, automatically solving all Design points in batch, judging whether each monomer is out of control thermally according to the Output Parameter value in each Design Point, and then determining the boundary value of the thermal conductivity coefficient capable of blocking thermal diffusion;
step S800, changing the thickness of the heat insulation plate, repeating the steps S200-S700, obtaining the heat conductivity coefficient boundary values corresponding to the heat insulation pads with different thicknesses, particularly, for the step S300, the heat insulation thermal runaway temperature-time change test data of the battery monomer does not need to be repeatedly fitted, and the obtained relevant parameters of the thermal runaway reaction kinetic equation of the battery monomer are directly given to the thermal runaway model of each battery monomer, and for the step S400, the validity and the accuracy of the simulation model do not need to be repeatedly verified; and sequentially obtaining the thickness and the thermal conductivity boundary value of the thermal insulation material required for blocking the thermal diffusion of the battery module.
Preferably, the material property parameters include density, specific heat capacity and thermal conductivity of the whole battery cell.
Preferably, the battery volume is obtained by adopting a drainage method, the weight is obtained by adopting a weighing method, and the material attribute parameter density is obtained according to the battery volume and the weight.
Preferably, the specific heat capacity of the whole battery cell is obtained by adiabatic calorimetry, that is, in an adiabatic environment created by an adiabatic calorimeter, an electric heating sheet with fixed heating power is used for heating the battery to raise the temperature of the battery, and the specific heat capacity change is obtained by the law of energy conservation.
Preferably, the overall thermal conductivity of the battery cell is obtained by a transient planar heat source method.
Preferably, the adiabatic thermal runaway temperature-time change test data of the whole battery monomer is obtained by performing an adiabatic thermal runaway test on the battery monomer in an acceleration calorimeter.
Preferably, the three-dimensional geometric model comprises a battery cell whole body of a battery monomer, a positive electrode lug, a negative electrode lug, a heat insulation plate, a heating plate and a clamping apparatus.
Compared with the prior art, the embodiment of the invention has the following specific beneficial effects:
in the product design or optimization improvement stage of the battery system, compared with the repeated thermal diffusion test method mainly adopted in the current stage, the method can quickly determine the boundary value of the thickness and the thermal conductivity coefficient of the thermal insulation material capable of blocking thermal diffusion, clearly determine the selection direction of the optimal thermal insulation material, save a large amount of time and cost, is safer and more environment-friendly compared with the traditional thermal diffusion test method, can promote the improvement of the safety application technical level of the power battery, and has the condition of high engineering application.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method according to an embodiment of the present invention;
FIG. 2 is a three-dimensional geometric model according to an embodiment of the present invention;
FIG. 3 illustrates a parameterized set of material thermal conductivity boundary values according to an embodiment of the present invention;
fig. 4 is a result of parametric analysis and calculation of the boundary value of the thermal conductivity of the material according to the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the 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.
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
The embodiment of the invention provides a method for determining a boundary value of the thickness and the thermal conductivity of a material for blocking thermal diffusion of a battery, and aims to quickly determine key attribute parameters and corresponding boundary thresholds of a thermal insulation material for blocking thermal diffusion of a power battery.
Referring to the attached drawing 1, taking a certain battery module composed of six square hard-shell battery monomers and a thermal baffle as an example, one battery monomer at the outermost side of the battery module is triggered by a heating plate to be thermally out of control, and the module and the heating plate are fixed by a clamp at the same time, the method of the embodiment of the invention is realized by the following steps:
and S100, obtaining material attribute parameters of the single batteries, test data of thermal insulation thermal runaway temperature-time change of the single batteries and temperature data of the center position of a contact surface between every two single batteries and the thermal insulation plate through tests.
Specifically, the battery monomer which is the same as the battery monomer in the module is selected, and material attribute parameters of the whole battery monomer, such as density, specific heat capacity and heat conductivity coefficient, are obtained through tests. Testing by adopting a drainage method to obtain the volume of the battery, obtaining the weight by adopting a weighing method, and obtaining the material attribute parameter density according to the volume and the weight of the battery, wherein the density is weight/volume; the specific heat capacity is obtained through adiabatic calorimetry, namely in an adiabatic environment created by an adiabatic calorimeter, an electric heating sheet with fixed heating power is adopted to heat the battery to enable the battery to rise in temperature, and the specific heat capacity change is obtained through the law of energy conservation; the thermal conductivity is obtained by a transient planar heat source method.
Specifically, an adiabatic thermal runaway test is carried out on the battery monomer in an acceleration calorimeter, and temperature-time change test data of the battery monomer in a complete thermal runaway process are obtained.
Specifically, a thermal diffusion test is carried out on the battery module, a thermocouple is arranged between every two battery monomers and at the center of the contact surface of the thermal insulation plate, and the heating plate is heated at a certain power until the adjacent battery monomers are triggered to be out of thermal control, and then the heating is stopped. And acquiring temperature data of different positions of the battery module captured by each thermocouple after the thermal diffusion process of the battery module is finished.
Step S200, building a three-dimensional geometric model for the battery module, importing the three-dimensional geometric model into Ansys Fluent Meshing, setting grid parameters according to the three-dimensional geometric model, dividing a surface network, and then generating a body network. Wherein, Ansys Fluent shifting is a Meshing tool for Ansys.
Specifically, a three-dimensional geometric model of a battery module in a thermal diffusion test is built through Ansys SCDM modeling software, the three-dimensional geometric model specifically comprises a battery monomer cell whole body 1, a positive electrode lug 2, a negative electrode lug 2, a heat insulation plate 3, a heating plate 4 and a fixture 5, only a fixing bolt with small influence on simulation is omitted, the details are shown in an attached drawing 2, the complete geometric model is subjected to shared topology, grid common nodes at subsequent interfaces are guaranteed, and the geometric model is stored. And (3) importing the three-dimensional geometric model into Ansys Fluent Meshing, setting reasonable grid parameters according to the size of the geometric model, reconstructing a surface grid, and generating a body grid.
Step S300, switching to an Ansys Fluent solving mode, setting a physical model, fitting the thermal insulation thermal runaway temperature-time change test data of the battery monomer to obtain the relevant parameters of the thermal runaway reaction kinetic equation of the battery monomer, and endowing each battery monomer with a thermal runaway model; setting corresponding material attribute parameters for different areas of the volume grid; setting heat generation power and heat generation time for the heating plate; environmental boundary conditions are set for surfaces in direct contact with the environment. Among them, Ansys Fluent is used for computational fluid dynamics simulation.
Specifically, a physical model setting is carried out by switching to an Ansys Fluent solving mode, a monomer thermal runaway model is applied to each battery monomer grid area, temperature-time test data obtained in a monomer thermal runaway test process are fitted through a built-in equation fitting tool of the Ansys Fluent, relevant parameters of a battery monomer thermal runaway reaction kinetic equation are obtained, and the thermal runaway model of each battery monomer is given. And setting corresponding material attribute parameters including density, specific heat capacity, heat conductivity coefficient and the like for different grid areas. And setting heat generation power for the heating plate, and setting heat generation time according to an actual heat diffusion test. A uniform natural convection environmental boundary condition is set for a surface in direct contact with the environment.
And S400, running simulation calculation, monitoring temperature values between every two battery monomers and at the center position of the contact surface of the heat insulation plate, comparing the temperature values with temperature data of different positions of the battery module obtained in the test, and verifying the effectiveness and accuracy of the simulation model.
Specifically, in the process of running simulation calculation, the temperature values of the central positions of the contact surfaces of every two battery monomers and the thermal insulation plate (the central positions are the same as the positions of all points measured by the thermocouples in the test process) are monitored, the temperature values are compared with the temperature data results of different positions of the battery module obtained in the test, and whether the simulation model meets the precision requirement is judged. Relevant parameters in the physical model setting can be properly adjusted to ensure the effectiveness and the accuracy of the simulation model.
Step S500, based on the verified simulation model, setting the thermal conductivity of the thermal insulation material as an Input Parameter in the Ansys Fluent, and setting the Parameter alpha of each monomer representing the thermal runaway reaction process as an Output Parameter.
Specifically, based on the simulation model which meets the precision requirement after verification, in the Ansys Fluent, the thermal conductivity of the thermal insulation board material is set as Input Parameter, a Report is defined to respectively monitor the Parameter α (the initial value of the Parameter is 1E-06, the Parameter α represents that no thermal runaway related reaction occurs in the monomer, in the thermal runaway process, the Parameter gradually changes from 1E-06 to 1, and 1 represents that the thermal reaction in the thermal runaway process of the monomer has completely occurred), and the Parameter α of six monomers is set as Output Parameter.
Step S600, loading a Fluent version block in an Ansys Workbench, importing the set simulation models of the Input Parameter and the Output Parameter, establishing a plurality of Design points (Design points) according to needs, modifying the specific value of the Input Parameter in each Design Point, and selecting different heat conductivity values in a certain range corresponding to the heat insulation board material. Wherein the Ansys Workbench is platform software of a simulation environment introduced by Ansys.
Specifically, in the Ansys Workbench, a Fluent version block is loaded, a simulation model with an Input Parameter and an Output Parameter set is imported, parallel solving is set, and an initialization mode is set. Entering a Parameter analysis setting interface, establishing a certain number of Design points according to the expected variation range and variation interval of the thermal conductivity of the thermal insulation material, modifying the specific value of the Input Parameter (corresponding to the thermal conductivity of the thermal insulation material) in each Design Point, and establishing 11 Design points in total as shown in figure 3, wherein the variation range of the thermal conductivity is from 0.135W/m K to 0.035W/m K, and the interval is 0.01W/m K.
Step S700, automatically solving all Design points in batch, judging whether each monomer is out of control thermally according to the Output Parameter value in each Design Point, and then determining the boundary value of the thermal conductivity coefficient capable of blocking thermal diffusion.
Specifically, the operation solution is run, that is, the operation of all Design points can be automatically and continuously completed, the simulation result is stored, and simultaneously the Output Parameter values (corresponding to the thermal runaway reaction process Parameter α) of the preset monomers in each Design Point are Output, and each Design Point can be determined according to the series of values, that is, whether thermal runaway occurs in each monomer under the thermal conductivity value of each thermal insulation material, and then the boundary value of the thermal conductivity coefficient capable of blocking thermal diffusion can be determined, as shown in fig. 4, when the thermal conductivity value is less than or equal to 0.065W/m K, except that thermal runaway occurs in the monomer directly contacting with the heating plate (the α value is 1), thermal runaway does not occur in the rest of monomers (the α value is infinitely close to the initial value 1E-06), that is, 0.065W/m K is the boundary value of the thermal conductivity coefficient capable of blocking thermal diffusion under the working condition, then the heat insulating material can be further optimized according to the value.
Step S800, changing the thickness of the heat insulation plate, repeating the steps S200-S700, obtaining the heat conductivity coefficient boundary values corresponding to the heat insulation pads with different thicknesses, particularly, for the step S300, the heat insulation thermal runaway temperature-time change test data of the battery monomer does not need to be repeatedly fitted, and the obtained relevant parameters of the thermal runaway reaction kinetic equation of the battery monomer are directly given to the thermal runaway model of each battery monomer, and for the step S400, the validity and the accuracy of the simulation model do not need to be repeatedly verified; and sequentially obtaining the thickness and the thermal conductivity boundary value of the thermal insulation material required for blocking the thermal diffusion of the battery module.
It should be noted that the english explanation mentioned in the present invention is related to ANSYS, and should not be understood as meaning other technical fields or industries.
It should be understood that various forms of the flows shown above, reordering, adding or deleting steps, may be used. For example, the steps described in the present application may be executed in parallel, sequentially, or in different orders, and are not limited herein as long as the desired results of the technical solutions disclosed in the present invention can be achieved.
The above-described embodiments should not be construed as limiting the scope of the invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions may be made in accordance with design requirements and other factors. Any modification, equivalent replacement, and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (7)
1. A method of determining a boundary value between a thickness of a material for blocking thermal diffusion of a battery and a thermal conductivity, the method comprising:
step S100, obtaining material attribute parameters of single batteries, test data of thermal insulation thermal runaway temperature-time change of the single batteries and temperature data of the central position of a contact surface between every two single batteries and a thermal insulation plate in a module thermal diffusion process through tests;
Step S200, building a three-dimensional geometric model for the battery module, importing the three-dimensional geometric model into Ansys Fluent Meshing, setting grid parameters according to the three-dimensional geometric model, dividing a surface network, and then generating a body network;
step S300, switching to an Ansys Fluent solving mode, setting a physical model, fitting the thermal insulation thermal runaway temperature-time change test data of the battery monomer to obtain the relevant parameters of the thermal runaway reaction kinetic equation of the battery monomer, and endowing each battery monomer with a thermal runaway model; setting corresponding material attribute parameters for different areas of the volume grid; setting heat generation power and heat generation time for the heating plate; setting an environmental boundary condition for a surface in direct contact with an environment;
step S400, running simulation calculation, monitoring temperature values between every two battery monomers and at the center position of a contact surface of the heat insulation plate, comparing the temperature values with temperature data of different positions of the battery module obtained in the test, and verifying the effectiveness and accuracy of a simulation model;
step S500, setting the thermal conductivity of a thermal insulation material as an Input Parameter in the Ansys Fluent based on the verified simulation model, and setting a Parameter alpha representing the thermal runaway reaction process of each monomer as an Output Parameter;
Step S600, loading a Fluent version block in Ansys Workbench, importing the set simulation models of the Input and Output parameters, establishing a plurality of Design points according to needs, modifying the specific value of the Input Parameter in each Design Point, and selecting different heat conductivity values in a certain range corresponding to the heat insulation board material;
step S700, automatically solving all Design points in batch, judging whether each monomer is out of control thermally according to the Output Parameter value in each Design Point, and then determining the boundary value of the thermal conductivity coefficient capable of blocking thermal diffusion;
step S800, changing the thickness of the heat insulation plate, repeating the steps S200-S700, obtaining the heat conductivity coefficient boundary value corresponding to the heat insulation plate capable of blocking heat diffusion under different thicknesses, particularly, for the step S300, no repeated fitting of the thermal runaway temperature-time change test data of the battery monomer is needed, directly endowing the obtained relevant parameters of the dynamic equation of the thermal runaway reaction of the battery monomer to each thermal runaway model of the battery monomer, and for the step S400, no repeated verification of the validity and the accuracy of the simulation model is needed; and sequentially obtaining the thickness and the thermal conductivity boundary value of the thermal insulation material required for blocking the thermal diffusion of the battery module.
2. The method of claim 1, wherein the material property parameters include density, specific heat capacity, and thermal conductivity of the battery cell as a whole.
3. The method of claim 2, wherein the cell volume is measured by a drainage method, the weight is measured by a weighing method, and the material property parameter density is obtained according to the cell volume and the weight.
4. The method as claimed in claim 2, wherein the specific heat capacity of the whole battery cell is obtained by adiabatic calorimetry, that is, in an adiabatic environment created by adiabatic calorimeter, an electric heating sheet with fixed heating power is used to heat the battery to raise the temperature of the battery, and the specific heat capacity change is obtained by the law of energy conservation.
5. The method of claim 2, wherein the thermal conductivity of the battery cell as a whole is obtained by a transient planar heat source method.
6. The method according to claim 1, wherein the adiabatic thermal runaway temperature-time variation test data of the battery cell as a whole is obtained by performing an adiabatic thermal runaway test on the battery cell in an acceleration calorimeter.
7. The method of claims 1-6, wherein the three-dimensional geometric model comprises a cell whole body of a battery cell, a positive electrode tab, a negative electrode tab, a heat insulation plate, a heating plate and a fixture.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110778476.8A CN113626946A (en) | 2021-07-09 | 2021-07-09 | Method for determining boundary value of thickness and thermal conductivity of material for blocking battery thermal diffusion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110778476.8A CN113626946A (en) | 2021-07-09 | 2021-07-09 | Method for determining boundary value of thickness and thermal conductivity of material for blocking battery thermal diffusion |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113626946A true CN113626946A (en) | 2021-11-09 |
Family
ID=78379430
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110778476.8A Pending CN113626946A (en) | 2021-07-09 | 2021-07-09 | Method for determining boundary value of thickness and thermal conductivity of material for blocking battery thermal diffusion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113626946A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115951242A (en) * | 2022-12-21 | 2023-04-11 | 国广顺能(上海)能源科技有限公司 | Battery testing method |
| CN116338473A (en) * | 2023-05-30 | 2023-06-27 | 苏州精控能源科技有限公司 | Cylindrical lithium battery testing method, device, medium and equipment for thermal runaway diffusion |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104008244A (en) * | 2014-05-29 | 2014-08-27 | 清华大学 | Design scheme for restraining thermal runaway extension of power battery module |
| JP2015215821A (en) * | 2014-05-13 | 2015-12-03 | 株式会社東芝 | Two-phase flow heat transfer analysis device and method |
| CN112231945A (en) * | 2020-09-15 | 2021-01-15 | 中国汽车技术研究中心有限公司 | Power battery system thermal diffusion joint simulation method based on star CCM + and Amesim |
| CN112446177A (en) * | 2020-11-16 | 2021-03-05 | 天华化工机械及自动化研究设计院有限公司 | Simulation method for heat insulation performance of external heat insulation material of high-temperature carbonization furnace |
| CN112711872A (en) * | 2021-03-29 | 2021-04-27 | 苏州清陶新能源科技有限公司 | Method for selecting type of battery module cell thermal insulation layer |
-
2021
- 2021-07-09 CN CN202110778476.8A patent/CN113626946A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015215821A (en) * | 2014-05-13 | 2015-12-03 | 株式会社東芝 | Two-phase flow heat transfer analysis device and method |
| CN104008244A (en) * | 2014-05-29 | 2014-08-27 | 清华大学 | Design scheme for restraining thermal runaway extension of power battery module |
| CN112231945A (en) * | 2020-09-15 | 2021-01-15 | 中国汽车技术研究中心有限公司 | Power battery system thermal diffusion joint simulation method based on star CCM + and Amesim |
| CN112446177A (en) * | 2020-11-16 | 2021-03-05 | 天华化工机械及自动化研究设计院有限公司 | Simulation method for heat insulation performance of external heat insulation material of high-temperature carbonization furnace |
| CN112711872A (en) * | 2021-03-29 | 2021-04-27 | 苏州清陶新能源科技有限公司 | Method for selecting type of battery module cell thermal insulation layer |
Non-Patent Citations (3)
| Title |
|---|
| 李煌: "三元锂离子电池热失控传播及阻隔机制研究", 《中国优秀博硕士学位论文全文数据库(电子期刊) 工程科技Ⅱ辑》, 15 January 2021 (2021-01-15), pages 1 - 140 * |
| 王晓娜;厉阳;侯德鑫;叶树亮;: "基于热成像的薄片材料热扩散率快速无损检测", 计量学报, no. 03 * |
| 虞跨海;李长浩;程永周;: "锂离子储能电池放电热行为仿真与实验研究", 电源技术, no. 01 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115951242A (en) * | 2022-12-21 | 2023-04-11 | 国广顺能(上海)能源科技有限公司 | Battery testing method |
| CN115951242B (en) * | 2022-12-21 | 2024-04-05 | 国广顺能(上海)能源科技有限公司 | Battery testing method |
| CN116338473A (en) * | 2023-05-30 | 2023-06-27 | 苏州精控能源科技有限公司 | Cylindrical lithium battery testing method, device, medium and equipment for thermal runaway diffusion |
| CN116338473B (en) * | 2023-05-30 | 2023-08-08 | 苏州精控能源科技有限公司 | Cylindrical lithium battery testing method, device, medium and equipment for thermal runaway diffusion |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | Extreme learning machine-based thermal model for lithium-ion batteries of electric vehicles under external short circuit | |
| Wang et al. | A unified modeling framework for lithium-ion batteries: An artificial neural network based thermal coupled equivalent circuit model approach | |
| Panchal et al. | Transient electrochemical heat transfer modeling and experimental validation of a large sized LiFePO4/graphite battery | |
| CN113626946A (en) | Method for determining boundary value of thickness and thermal conductivity of material for blocking battery thermal diffusion | |
| Jaguemont et al. | Characterization and modeling of a hybrid-electric-vehicle lithium-ion battery pack at low temperatures | |
| Hu et al. | A foster network thermal model for HEV/EV battery modeling | |
| Xu et al. | Modelling of thermal runaway propagation in lithium-ion battery pack using reduced-order model | |
| CN111475933A (en) | Lithium ion battery thermal runaway simulation method and system | |
| Xie et al. | An improved resistance-based thermal model for prismatic lithium-ion battery charging | |
| CN111597747B (en) | Multi-point triggering ternary lithium power battery module thermal runaway simulation and prediction method | |
| CN110633496B (en) | Method for determining thermal stress and temperature in lithium ion battery discharging process based on thermal-force coupling model | |
| CN111144029B (en) | Modeling method for thermoelectric coupling characteristics of lithium ion power battery | |
| Xie et al. | An improved electrothermal‐coupled model for the temperature estimation of an air‐cooled battery pack | |
| CN113794254A (en) | Thermal management strategy configuration method and device, computer equipment and storage medium | |
| Cao et al. | A full-scale electrical-thermal-fluidic coupling model for li-ion battery energy storage systems | |
| Samad et al. | Parameterization of battery electrothermal models coupled with finite element flow models for cooling | |
| Muratori et al. | A spatially-reduced dynamic model for the thermal characterisation of Li-ion battery cells | |
| CN111767625A (en) | Thermal simulation method for lithium ion battery pack | |
| CN117826615B (en) | Method for determining control parameters of cooling liquid of power battery of electric automobile | |
| Moralı | Computational modeling and statistical evaluation of thermal behavior of cylindrical lithium-ion battery | |
| CN114696005A (en) | Temperature difference testing method, system and device for liquid cooling system of power battery | |
| CN113835031B (en) | Information processing method, apparatus, electronic device and storage medium | |
| CN117272745A (en) | Thermal coupling simulation method combining thermal model and high-temperature creep analysis | |
| CN116679212A (en) | Simulation method, system, equipment and medium for battery quick charge | |
| Cicconi et al. | Modeling and thermal simulation of a PHEV battery module with cylindrical LFP cells |
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 |