基于Fluent的单体锂离子电池热仿真分析

2021-10-08 来源：爱问旅游网

NEW ENERGY AUTOMOBILE | 新能源汽车

詹世东　李康伟　蔡友彬　石浩宇　吴彦鑫

广州城市理工学院汽车与交通工程学院广东省广州市 510000

摘要：本文采用Fluent软件中的MSMD电池模型进行仿真，分析了不同放电倍率下单体锂离子电池的电压和温度特性。首先，根据实际尺寸建立单体锂离子电池的三维模型，并进行网格划分及边界条件命名等前处理；其次，利用Fluent的MSMD电池模块对电池单体在300K室温条件下分别以1C、3C、5C、7C的放电倍率进行温度仿真，得到单体电池的温度特性曲线和电压特性曲线，分析结果表明，随着放电倍率的增大，电池的放电终止电压越低，耗电量越高，电池的温度和温差都随之增大。关键词：锂离子电池　MSMD　单体电芯Thermal simulation analysis of Lithium-ion cell based on FluentZhan Shidong　Li Kangwei　Cai Youbin　Shi Haoyu　Wu Yanxin

This paper uses the MSMD battery model simulation in the Fluent software to analyze the voltage and temperature characteristics of lithium-ion cell under Abstract： diﬀ erent discharge rates.First, a three-dimensional model of the lithium-ion battery cell according to the actual size is established, and pre-processing such as grid division and boundary condition naming are conducted. Secondly, the MSMD battery module in the ﬂ uent software is used to simulate the temperature of the battery cells at the discharge rates of 1C, 3C, 5C, and 7C at 300K room temperature. Finally, the temperature characteristic curve and voltage characteristic curve of lithium-ion cell are obtained. The analysis results show that as the discharge rate increases, the lower the discharge end voltage of the battery and the higher the power consumption, the temperature and temperature diﬀ erence of the battery increase accordingly.Key words：Lithium-ion battery; MSMD; Cell1　前言

目前，随着环境污染及能源短缺问题在全球范围内的出现，寻找更环保、低耗和安全的汽车能源显得尤为关键，因此推动电动汽车的发展成了人们普遍关注的重点。锂离子电池的高能量密度、循环寿命长、安全性高等优点被广泛应用于电动汽车行业。但是，锂离子电池的安全性能受充放电过程中的温度影响较大，锂离子电池最佳工作温度范围为20℃～40℃，最大温差不超过5℃。单体电池的温度和温差过大会导致锂离子电池的充放电效率降低，温差严重时还会产生电池的热失控和安全问题。

闫凯和楚金甫等人[1]利用CFD软件对18650单体电芯进行了1C放电工况下的热仿真分析，并使用热成像对比验证二者的测试结果，得出了CFD软件热仿真分析的数据相对准确的结论；北京建筑大学的周庆辉和陈展等人[2]利用Fluent软件进行了不同电流下锂离子电池的温度场仿真，得出电池电流越大对电池的温度和温差的影响越大；南昌大学的李甜甜等人[3]进行单体电池的温度场仿真分析，并对单体电池进行对应的放电实验得到了电池实际温度的相关数据，结合仿真

得到的数据，得出所建立的电池三维热模型具有真实性。

本文以研究单体锂离子电池热仿真为目的，采用基于Fluent的MSMD电池模块建立单体电池生热模型，对不同放电倍率下的单体电池的温度及电压情况进行研究。

图1　单体电池模型

2　单体电芯建模

2.1　单体锂离子电池模型建立

利用CATIA软件对格瑞普公司生产的软包钴酸锂16.8Ah单体电芯进行三维建模。因原电池模型较为复杂，本文将原模型的凸台、倒角、凹槽等与电池放电生热影响不大的细节删除，将原电池模型进行简化，得到简化后的电池模型尺寸，如图1所示。电芯的长度为154.5mm、宽度为85mm、厚度为10.5mm。极耳长度为11.5mm、宽度为25mm、两极耳间距为6.5mm。在SpaceClaim中对几何模型进一步进行简化，将电池的极耳与电芯部分设置为独立的几何结构，并对各部分进行命名。

2.2　电池模型的网格划分与前处理因电池模型简化后为规则的长方体几何结构，故在ANSYS Mesh中将网格类型设置

为六面体网格，网格尺寸控制为3mm，其余网格划分选项不作修改，对电池模型进行网格划分。最终网格平均质量为0.97831，网格单元数为4348，该网格质量可以保证在网格整体不失真的情况下进行温度仿真，且同时保证了数据精度与仿真速度。

对单体锂离子电池模型的边界进行命名，在ANSYS Mesh中将电芯外表面定义为wall-cell，将极耳顶面定义为tap-n与tap-p，将极耳的其余面定义为wall-n与wall-p。

128

AUTO TIME

时代汽车 www.cnautotime.com3　锂离子电池热仿真基础

3.1　锂离子电池热特性分析

锂离子电池生热主要是指充放电过程中因电池内部发生化学反应使得锂离子在电池正负极材料之间来回运动所产生的反应热Qa，这部分热量占电池生热量的绝大部分，除外，锂离子电池的生热量还包括欧姆内阻的焦耳热Qb、极化反应的极化热Qc。所以锂离子电池在充放电过程中所产生的总热量应为[5]：

(1)式中，n为电池的数量；m为电池极耳的质量；Q为电池的电化学反应热；I为充放电电流；M为摩尔质量；F为法拉第常数；Re为欧姆内阻；Rp为极化内阻。

3.2　MSMD电池模型的微分方程在Fluent中的MSMD模块，采用微分方程求解电池热场和电场：

[6]

[4]

表1　钴酸锂电芯材料参数

参数平均密度/（kg·m）-3导热系数不同，所以电池最高温度位于电芯中间，最低温度位于极耳附近。由表2可知，

数值2040随着放电倍率的增大，单体电芯最高温度、最低温度、温差都随之相应增大。在7C放电倍率下，单体电芯温差最大，为1.81K。

平均比热容/（J·（kg·K））-11200-1X方向平均导热系数/（W·（m·K））Y方向平均导热系数/（W·（m·K）-1）Z方向平均导热系数/（W·（m·K）-1）扩散系数uds0扩散系数uds114.20.7514.21e72e8图3　单体电芯电压特性曲线

3.83.6电池电压(V)3.43.23.02.80

200

4006008001000120014001600

放电时间(S)

1C3C5C7C

时间步数设为50，单位步长为30s，分别以1C、3C、5C、7C的放电倍率进行单体电芯温度场仿真，经计算得出各放电倍率下单体电芯的壁面温度云图、电压特性曲线和温度特性曲线，如图2、3、4所示。

4.2　不同放电倍率下的仿真结果

(2)

　　　 (3)

(4)

上式中，σ+和σ−为正负极的有效电导率，

图2　不同放电倍率下单体电芯壁面温度云图

由图3可知，随着放电倍率的逐渐增大，电池放电终止电压逐渐降低，放电终止时间逐渐减小，放电电压曲线的斜率不断增大，电池电量耗尽时间逐渐减小。由图4可知，随着放电倍率的逐渐增大，电池放电终止温度逐渐升高，斜率增大，电池升温速度加快。

ECh为电ϕ+和ϕ−为正负极的相电势，jECh和q化学反应产生的体积电流传递速率和电化学

short分别为电池内部短路引反应热，jshort和qabuse为热滥用条起的电流传递率和产热率，q件下热失控反应产生的热量。

(a)1C放电倍率放电(b)3C放电倍率放电

图4　单体电芯温度特性曲线

325320壁面温度(K)3153103053000

200

4006008001000120014001600

放电时间(S)

1C3C5C7C

4　单体电芯放电生热仿真

4.1　单体电芯温度仿真

本文采用ANSYS Fluent的MSMD 电池模块建立单体电芯放电生热温度场仿真模型，在Fluent软件中打开能量方程，设置为瞬态求解。在Fluent求解器中分别输入以下UDF指令打开MSMD电池模型[7]：

define/models/addon-module/8，在MSMD电池模型中选择 NTGK 模型，设置电芯容量为16.8 A.h，在Conductive Zones中分别定义电芯及正负极，在Electric Contacts中定义正负极极耳连接面。

设置正负极极耳材料为铝，负极极耳设置扩散系数uds0与uds1，电芯材料属性设置如表1所示，设置初始温度为300K，电芯与极耳的壁面在自然条件下对流换热系数设置为5W/(m·K)，设置电池电压和温度监测，

从图2可以得到在1C、3C、5C和7C放电倍率下的电池最低温度和最高温度以及温差的对比数据，如表2所示。

表2　不同放电倍率下的温度数据

放电倍率1C3C5C7C最低温度/K302.23307.49311.57314.71最高温度/K302.46308.31312.91316.52温差/K0.230.821.341.814.3　仿真结果分析

从锂离子电池在不同的放电倍率下进行的放电生热分析仿真结果看出，当放电倍率增大时，电池的电压初始电压和终止电压减低，电池电压下降速度加快，电池内阻增大，电池内部生热量随之增大，导致电池各部分产生的温度和温差增大。若以更大放电倍率对电池进行放电，或将电池串联组成电池包放电使用时，将可能导致电池温度超过电池的合理工作温度范围，如果不及时采取适当

因为电池放电时大部分热量由电池内部电化学反应产生，小部分热量由内阻热产生，热量以热传导的形式向外扩散，电芯各方向

AUTO TIME

129

的散热措施，将影响电池的正常放电工作，并且会危害到电池的使用寿命，甚至出现热失控等安全隐患。

基金项目：广东省攀登计划科技创新战略专项（项目编号：pdjh2020a0820）。

[5]方雄灿. 电动汽车锂离子电池组热特性分析及电池箱体散热结构优化[D].合肥工业大学， 2017.

[6]薛宁.动力电池热行为及其过热性分析[D].吉林大学， 2017.

5　结束语

1）随着锂离子电池放电倍率的增大，电池的终止电压降低，电量下降速度加快，电池的最低与最高温度增大，温差也增大。

2）对锂离子电池进行不同放电倍率下的温度仿真得出，电池中心部分温度最高，极耳部分温度最低，由此可根据不同放电倍率下的温度变化进行对电池包的结构散热进行设计及优化。

3）在进行电池1C、3C、5C的放电倍率放电时，电池温度逐渐升高，但都低于40℃，以7C放电倍率进行放电时电池最低温度和最高温度都超过电池放电的合理工作范围，此时需及时采取散热措施，保证电池使用寿命。

参考文献：

[1]闫凯，楚金甫，陈西山，古伟鹏，常乐.纯电动车用动力电池系统温度场特性的数值研究[J].时代汽车，2018（02）：35-37.[2]周庆辉，陈展，刘少楠.基于Fluent的锂离子动力电池的热分析[J].北京建筑大学学报， 2017，33（03）：33-37+65.[3]李甜甜. 电动汽车锂离子电池热特性分析及散热优化[D]. 南昌大学， 2016.

[4]田华，王伟光，舒歌群，严南华.基于多尺度、电化学-热耦合模型的锂离子电池生热特性分析[J].天津大学学报（自然科学与工程技术版）， 2016， 49（07）：734-741.

[7]郭磊.对FLUENT辐射模型的数值计算与分析[J].制冷与空调（四川）， 2014，28（03）：358-360.

作者简介

詹世东：（1999.03—），男，汉族，本科。研究

方向：锂离子电池热管理。

（上接第127页）

根据此次实验结果，在现实生活中选择合适的空气间隙，能够使变压器的传输效率得到改善。

参考文献：

[1]翟鹏伟. 电动汽车无线充电系统设计[D]. 北京, 北方工业大学， 2015.

[2]李松林. 基于电磁感应耦合的无线电能传输研究[D] 成都电子科技大学， 2012. [3]陈珂睿，王泽忠，刘胜南. 非接触式电能传输系统功率及效率影响因素[J]. 电网技术， 2014， 38（3）：807-811. [4]杨民生王耀南新型无接触感应耦合电能传输技术研究综述。湖南文理学院学报.[5]Buja G，BertoluzzoM，Dashora HK.Lumped track layout design for dynamic wireless charging of electric vehicles[J]. IEEE

Transactions on Industrial Electronics， 2016， 63（10）： 6631-6640.

[6]朱华松耦合无线能量传输系统相关技术研究[D] 南京南京邮电大学， 2013.

5　总结

本文针对电动汽车松耦合变压器存在的原副边耦合系数低的问题，建立变压器磁路模型，通过补偿方案提高变压器传输效率。通过仿真与实验结果表明：所设计的串串补偿的松耦合 LLC-DCT 和全桥 DCT 相比，可以提高能量传输效率，优化后输出效率体提高了1.86%。

作者简介

吴金华：（1983.10—），男，汉族、江西鹰潭人，

汽车学院，硕士学位、副教授。研究方向：新能源汽车技术、控制理论及自动化技术。

130

AUTO TIME

因篇幅问题不能全部显示，请点此查看更多更全内容

查看全文

全部栏目

基于Fluent的单体锂离子电池热仿真分析