基于star-ccm+＆Ａmesim联合仿真的液冷电池包热管理

1年前浏览4055

今天看到一篇文章写的不错，来和大家分享一下。

以某车用锂离子电池组为研究对象，主要研究了爬坡工况、９０ｋｍ／ｈ匀速工况和ＮＥＤＣ三种工况下动力电池组的温升情况。利用ＳＴＡＲ－ＣＣＭ＋和Ａｍｅｓｉｍ软件联合对液冷电池包进行热管理仿真，分析流场和温度场的分布情况，预测综合工况下电池包模组的最高温度和模组间温差分布，并通过热管理试验验证三种工况下试验结果与仿真结果是否吻合，以提高仿真精度。

动力电池包内热量的累积不仅影响电池的使用效率及使用寿命，同时易造成动力电池系统故障并引发安全事故［，因此准确预测电池包内温度分布，并对温度场进行分析具有重要意义。动力电池包热管理系统设计中，通常结合仿真来预测电池包的温度分布、冷却系统的流量分配和压力分布等，从而预测热管理系统的性能。仿真一般分为３Ｄ仿真和１Ｄ仿真，３Ｄ仿真可用于电池包液冷板流场和压力场的仿真，以及模组温度场的仿真，以获得流场和温度场的细节，但３Ｄ仿真软件计算瞬态工况耗时较长，不便或无法用于系统级别仿真以及控制策略仿真；１Ｄ仿真从系统角度出发，模型从电池包扩展至包含整个冷却／加热系统外部环路等，由于建模中对各相应部件进行了简化，在对系统性能进行仿真的时候，能大大提高仿真速度，通常用于系统级别的瞬态循环工况仿真和制定电池包热管理控制策略等。

本文主要描述电池热管理设计中的３Ｄ＆１Ｄ联合仿真分析方法：利用３Ｄ稳态计算结果对１Ｄ模型进行流场和温度场标定，并制定冷却／加热控制策略，然后对整体电池热管理系统进行瞬态工况仿真分析，快速预测热管理系统性能，并通过试验检验仿真结果的准确性。

１模型简介

本文以某项目电池包为例，此电池包是一款新能源混合动力电池包。其平台电压为３５０Ｖ，共包含８个电池模组，每个电池模组包含１２片电芯，共计９６片电芯；单体电芯标称容量为５１Ａｈ，系统总容量４８９６Ａｈ，系统总能量１７ｋＷｈ。电池包采用箱体底部集成液冷系统的冷却方式，所有电池模组均放置于底部液冷板上，并通过导热硅胶垫进行热量传导。电池包结构布置如图１所示。

图１ＰＡＣＫ结构示意

基于star-ccm+的３Ｄ仿真

star-ccm+是ＣＤ－ａｄａｐｃｏ公司采用最先进的连续介质力学数值技术开发的新一代 CFD仿真软件，它将现代软件工程技术、最先进的连续介质力学数值技术和卓越的设计结合在一起，解决跨学科问题的综合工程解决方案。基于star-ccm+的３Ｄ仿真的主要目的是获得电池包的温度场以及液冷板的流场特性，具体实施方法如下。

２．１模型简化与网格划分

本文以某项目电池包为例，模型主要考虑模组和液冷板，并根据各自几何模型适当进行简化。如图２所示，其中，模组部分在简化了电芯的同时，保留极

耳、铜排、泡棉、绝缘膜和导热硅胶垫等细节，根据各个部分的几何特性和仿真需求差异采用了不同的网格参数设置，比如在铜排和导热硅胶垫等薄体几何区域采用拉伸薄体网格设置，在保证细节的同时也能降低网格数量；根据数模搭建液冷板仿真模型，抽取流体域并划分网格，根据实际需求在壁面生成边界层以更好的模拟真实流场。模型简化和网格划分都在Ｓｔａｒ－ＣＣＭ＋中进行，处理后的总体网格数量约1300W.

图２模组结构、水流道和网格划分

２２求解设置与边界条件

模型求解设置主要为：（模型采用稳态计算模式；

（流场和温度场耦合仿真；（湍流模型为ｋ－ε模型；（发热体形式为恒密度恒功率发热源；（流体为恒密度不可压缩流体。

边界条件设置主要为：①环境温度设置为２５℃；

②初始温度为２０℃；③冷却液流量分别为ｑｖ１、ｑｖ２和ｑｖ３，其中，ｑｖ１为５Ｌ／ｍｉｎ，ｑｖ２为１０Ｌ／ｍｉｎ，ｑｖ３为１５Ｌ／ｍｉｎ；冷却液温度设置为２０℃。

２．３仿真及结果分析

设置完成后分别求解三种流量下的流场和温度场，并根据不同流量下的流阻绘制阻力特性曲线图用于环路系统水泵性能匹配。

对于液冷电池包而言，并联支路的冷却液流量分配直接影响模组的温差，冷却液流动状态直接影响模组换热的强弱，进而影响电芯温度分布。因此，在进行热管理系统设计过程中需根据仿真结果，进行并联

支路流量分配计算，对于支路流量差别较大（＞５％）的情况进行优化设计以消除随之带来的模组温差；并对流道中出现回流区、分离区或流通截面急剧变化的区域进行优化设计，以降低局部流阻。

本例中，对流场和温度场仿真结果分析如下：

１）流场仿真结果分析

设定环境温度为２０℃，冷却液为５０％乙二醇水溶液，进口冷却液温度为２０℃，以冷却液流量ｑｖ１、ｑｖ２和ｑｖ３对冷却系统内流场进行仿真分析。当系统内冷却液流量为ｑｖ３时，模拟水冷板内各支路流量分配情况如图３所示。可以看出，各流道内冷却液流量分配相对均匀，两支路流量占比分别为５２％和４８％，其各支路流道内流量差为±２％（＜５％），满足系统支路流量分配设计要求；图４为模拟冷却液流量为ｑｖ３时水冷板底部截面冷却液流速分布示意图，由于液冷板采用铝挤成型加堵头搅拌摩擦焊的工艺，在流道转弯处形成直角转弯，导致此处冷却液流速变化较大，并且在转弯位置容易形成回流区，造成局部流阻较高，因此在直角转弯处流道存在优化空间，以降低因冷却液流动急剧变化对电芯温度分布的影响。

图３流量分配示意

图４水冷板底部截面流速示意

图５为液冷系统最终绘制的流阻特性曲线图，随着系统内冷却液流量的增加，系统内压降呈快速递增趋势，当系统内冷却液流量为ｑｖ１时，系统压降为２３４４Ｐａ，当系统内冷却液流量增加至ｑｖ２和ｑｖ３时，系统压降分别增加了２２１％和５７２％，为７．５ｋＰａ和１５．７ｋＰａ。

２）温度场结果分析

通过仿真不同流量下电池系统的温度场，得到电池温度与流量的关系。图６显示为假定环境温度为０℃，单体电芯发热量为恒定值５．２Ｗ，冷却液入口温度为２０℃时，在ｑｖ１、ｑｖ２和ｑｖ３三种流量下电池包内模组温度分布云图。可以看出，系统内冷却液流量不同时电池模组温度也不相同。冷却液流量越小，则在相同时间内冷却液可带走的模组热量越少，电池模组的温度相对越高。当冷却液流量为ｑｖ１时，电池最高温度为３６℃；当冷却液流量为ｑｖ３时，电池最高温度为３２℃，可见在一定程度上增大冷却液流量可以有效降低电池温度。为进一步了解单体电芯温度分布情况，对ｑｖ１、ｑｖ２和ｑｖ３三种流量下单电芯温度分布情况进行模拟，得到温度分布云图如图７所示。分析可知，单体电池温度分布呈从上至下递减趋势，靠近水冷板侧温度较低，远离水冷板一侧电池温度较高，不同流量下电池上下温差差别不大，一般为４℃左右。不同流量下单体电芯温度也不相同，冷却液流量越大，电芯温度越低，不同流量下电芯之间的最大温差为２．１℃。

图５液冷板阻力特性曲线

图６不同流量下电池温度云图

图７不同流量下单电芯温度分布

图８所示为冷却液流量为ｑｖ３情况下，８个模组中电芯和铜排的平均温度。从温度分布看，各电芯间和铜排间温度分布相对均匀，不同模组之间电芯和铜排的温差均小于０．５℃，可见液冷板均温性较好。其中，电芯的平均温度约为２９．８℃，铜排的平均温度高于电芯温度约２．４℃，为３２．２℃，一方面是因为电芯整体热容较大，在热量传递过程中温升较小，而铜排由于其自身内阻产热且热容相对较小，因此温升较高；另一方面是因为电芯通过水冷板传导部分热量，导致其温度下降较快。当冷却液流量为ｑｖ３时冷却液与冷板交界面换热系数分布如图９所示，交界面上平均换热系数为５８４．９Ｗ／ｍ２（Ｋ，进口处流道内冷却液流速相对较快，温度相对较低，其换热系数高于平均数，随着冷却液的持续流动和换热，冷却液温度逐渐升高，换热系数也随之降低。其中，在冷板底部流道转弯处由于其方形结构造成冷却液局部流速较快，其换热系数也高于平均数。

图８电芯和铜排温度分布

图９冷却液与冷板交界面换热系数分布

３基于ＡＭＥｓｉｍ的１Ｄ仿真

ＡＭＥＳｉｍ是一款多学科领域的复杂系统建模与仿真平台，其智能求解器能够根据用户所建模型的数学特性自动选择最佳积分算法，并根据不同仿真时刻

的系统特点动态地切换积分算法和调整积分补偿以缩短仿真时间和提高仿真精度。基于ＡＭＥｓｉｍ的１Ｄ仿真的主要目的是获得整个电池包瞬态循环工况的温度分布以及用于控制策略仿真，具体实施方法如下：

３．１１Ｄ模型搭建

本例中由于模组温度分布并不是主要关心问题，故将模组离散成电芯本体和极耳及汇流排两个热容体，热容体之间用热传导元件连接，根据电池包的几何尺寸和材料属性设置各元件参数，并且施加３Ｄ仿真相同的边界条件进行流场和温度场标定，然后分别计算相应的平均温度。底部采用换热管道与模组进行换热，且只计算模组法向方向上的传热量，通过调整换热元件壁面的相对粗糙度（ｒｅｌａｔｉｖｅｒｏｕｇｈｎｅｓｓ）和管件的努赛尔数表达式（Ｎｕｓｓｅｌｔｎｕｍｂｅｒｅｘｐｒｅｓ－ｓｉｏｎ），使１Ｄ模型的阻力特性和电芯温度与３Ｄ结果保持一致。具体如图１０所示。