首页 > 行业资讯 > 《储能科学与技术》推荐|彭可 等：​基于有限元的热力耦合场匣钵运动分析与优化

《储能科学与技术》推荐|彭可 等：​基于有限元的热力耦合场匣钵运动分析与优化

时间：2024-04-04 来源： 浏览：

《储能科学与技术》推荐|彭可 等：​基于有限元的热力耦合场匣钵运动分析与优化

原创 彭可 张志成 等 储能科学与技术
储能科学与技术

esst2012

中文核心、科技核心和cscd核心期刊，化学工业出版社和中国化工学会主办，主编黄学杰研究员。投稿及下载官网：http://esst.cip.com.cn/CN/2095-4239/home.shtml；欢迎给公众号投稿

作者： 彭可   张志成 胡有章 张旭辉 周稼辉 李彬

单位： 湖南师范大学工程与设计学院

引用： 彭可 , 张志成 , 胡有章 , 等 . 基于有限元的热力耦合场匣钵运动分析与优化 [J]. 储能科学与技术 , 2024, 13(2): 634-642. 

DOI ： 10.19799/j.cnki.2095-4239.2023.0496

本文亮点： 1.分析锂离子电池正极材料烧结过程中辊道窑内匣钵异常横向运动的影响因素并提供了分析数据。2.设计了一种改进机构，有助于推进锂离子电池正极材料生产设备的发展。

摘 要   为解决装有锂离子电池正极材料的匣钵在辊道窑烧结过程中易出现异常横向运动使匣钵破损的问题，本研究通过分析匣钵异常运动特性的影响因素对其进行优化。首先，利用Solidworks软件建立辊道-匣钵传动系统简化模型，对匣钵进行受力分析，推导出匣钵发生异常横向运动的原因是辊棒的弹性变形；其次，通过Abaqus软件进行有限元仿真，对比分析热力耦合场和重力场下匣钵的位移曲线以及辊棒的受力情况。结果表明，与重力场相比，热力耦合场的辊棒弹性变形量更大，且热力耦合场匣钵横向位移在0～1 s内相对重力场增大，而在1～5 s内相对重力场减小，说明在温度和重力共同作用下所引起的辊棒弹性变形是导致匣钵异常横向运动的主要原因。在此基础上，设计了一种匣钵传动机构，包括辊棒支撑机构、匣钵夹紧装置及匣钵挡齐装置，并通过对比生产验证了改进机构的可行性。本研究对在辊道窑烧结过程中匣钵异常横向运动的影响因素进行分析并提供了一种改进机构，有助于推进锂离子电池正极材料生产设备的发展。
关键词   锂离子电池； 正极材料； 辊道窑； 匣钵； 热力耦合； 有限元法
目前，锂电池产业已成为清洁能源产业的重要组成部分，广泛应用在储能、动力等领域。锂离子电池正极材料作为锂离子电池的关键材料，随着储能等行业的快速发展，市场对其质量要求越来越高，因此我国锂离子电池正极材料的生产工艺和生产设备的技术水平亟需得到进一步提升。在锂离子电池正极材料烧结自动化生产线上，烧结工艺是锂离子电池正极材料生产的关键一步。目前市场上主要的烧结反应设备是辊道窑，辊道窑是以旋转辊道为载体进行连续烧制的窑炉，辊道窑内的温度场是影响烧结产品质量的重要因素。锂离子电池正极材料在辊道窑烧结过程中，匣钵作为装载锂离子电池正极材料的载体，伴随着辊道窑中辊棒的转动向前移动。但匣钵在向前运动过程中容易出现横向运动现象，偏离其预期直线，从而使得匣钵在辊道窑烧结过程中受热不均，导致匣钵损坏和锂离子电池正极材料的质量不稳定。尤其当匣钵异常运动行为严重时，轻则窑内产生拥堵，重则导致匣钵破碎，损坏生产设备，增大生产成本，如图1所示。因此，研究辊道窑中匣钵产生异常横向运动的原因并对此加以改善，对锂电池产业的发展具有重要的意义。

图1     匣钵破碎及其碎片
近年来，为解决辊道窑烧结过程中的匣钵运动异常问题，有学者对其原因进行了一定的研究并提出了相关解决技术和方法。Darko Krmpotic通过在双通道陶瓷辊道窑温度不稳定区进行快速检测，发现在辊道窑温度高于1100 ℃的最热部分易出现生产故障，当温度下降超过5 ℃时，该位置会发生质量故障，因此温度的分布不均匀会导致辊道窑内出现异常情况。田力通过更改辊道窑的工艺配方参数来影响温度均匀性，并采用BP神经网络和遗传算法得到最优工艺配方参数，从而改善辊道窑温度场的均匀性。姜永正等基于柔性动力学分析柔性辊道对匣钵运动特性的影响并对辊棒进行刚柔性对比，得出了辊棒在重力场下产生的柔性变形是导致匣钵出现横向运动的主要原因，通过增大辊棒的刚度来减轻烧结过程中匣钵的异常横向运动现象，但未考虑辊道窑中的温度场对匣钵异常横向运动的影响。奚慧春等认为辊棒产生弯曲变形的原因是其受到了高温热载荷和压力载荷，通过对辊棒进行受力分析并提出了一种变截面辊棒，使得辊棒各截面的弯曲正应力一致，从而改善辊棒的抗弯曲能力。Michal等采用多体动力学方法对刚体在辊道上的运动进行数值模拟分析，结果表明造成刚体异常偏移运动的主要原因是承载刚体的辊棒的弯曲程度大于空载辊棒的弯曲程度，且辊棒中部的弹性变形量较外侧更大，使得不同位置刚体的接触力和摩擦力不同，导致辊棒中部刚体的轨迹差大于外侧刚体的轨迹差，通过在辊道中间增加支撑和提高辊棒刚性的方法，改善辊道输送过程中的异常偏移运动。
综上所述，以往对于匣钵异常运动行为的研究主要以单独温度场、单独重力场以及辊棒的受力分析为主，并通过提高辊道窑内温度场均匀性或改进辊棒的结构尺寸来改善匣钵异常运动现象。但仅单独考虑重力场或温度场存在较大的局限性，未能真实反映出辊棒的受力变形情况和匣钵的运动状态，使得辊棒挠度的调整不合理，解决匣钵异常运动问题的效果较差。因此，现需研究在温度场和重力场共同作用下对匣钵异常运动的影响，进一步改善匣钵异常运动现象。
本工作基于上述研究，利用有限元法对热力耦合场下的匣钵运动特性进行分析与优化。

1 辊道-匣钵传动模型

本工作以辊道-匣钵传动简化模型作为研究对象。辊道窑主体部分—烧结炉由于温度较高，为满足辊棒的耐热性要求，本工作辊棒的材料选取碳化硅，弹性模量和泊松比分别为100 GPa和0.2。采用的辊棒直径为50 mm，辊棒总长度为2545 mm，由于辊棒需安装在轴承上，辊棒实际有效支撑长度为2397 mm，相邻辊棒的轴心间距相同且均为110 mm。单个匣钵尺寸为330 mm×330 mm×330 mm，装满锂离子电池正极材料后单个匣钵总质量为40 kg，匣钵初始间隙为1 mm。为简化计算，对辊道-匣钵传动结构进行简化，选取4个匣钵以1 mm为间距作为一排并放置在辊道上方，再通过Solidworks进行建模，模型如图2、图3所示。

图2     辊道 - 匣钵传动简化模型

图3     传动模型放大图

2 匣钵受力分析

因匣钵质量分布均匀且视为刚体，可将匣钵视为一个质点。考虑辊棒的弹性变形，辊道两端与 a 点间的夹角分别为 和 ，匣钵的质量为 ，辊道与匣钵间的摩擦系数为 ，匣钵在辊道上所受的支持力为 ，匣钵在辊道上的受力正视图如图4所示。

图4     匣钵受力正视图
若匣钵与辊道间没有发生相对滑动，匣钵所受的静摩擦力 的表达式如式(1)所示：
(1)
若匣钵与辊道间发生相对滑动，匣钵所受的滑动摩擦力 的表达式如式(2)所示， 的表达式如式(3)所示：
(2)
(3)
由于辊道弹性变形过大使得夹角 增大，导致 大于最大静摩擦力，匣钵会在辊道上朝着 a 点产生相对滑动。
匣钵在辊送时其受力分析如图5所示。由简化模型的辊棒轴心间距和匣钵尺寸可知，匣钵在辊送过程中绝大多数状态是由3根辊棒支撑，少数状态下由4根辊棒支撑，两种状态匣钵的受力分析相似。本工作为便于分析，示意图采用了3根辊棒支撑匣钵的情况。将匣钵视为刚体并考虑辊棒的弹性变形，因此匣钵与辊棒接触部分在左右侧边，左右两侧接触部分的支持力分别为 和 。辊棒在匣钵两侧的变形程度不同会导致支持力不同。由图5可知辊棒在匣钵右侧的弯曲程度更大，重力在垂直于辊棒方向上的分力会较小，因此 大于 ，数学表达式如式(4)所示：
(4)

图5     匣钵辊送受力分析
由于匣钵和辊道在两侧接触部分的摩擦系数相同，因此两侧动摩擦力关系式如式(5)所示：
(5)
此时，两侧摩擦力不相等，其绕匣钵质心 A 的力矩也不相等，所以会出现一个绕 Y 轴的力矩，该绕 Y 轴的力矩 表达式如式(6)所示：
(6)
由图5可以看出，在旋转力矩 的作用下，匣钵在辊道上有做顺时针旋转的趋势，从而导致匣钵出现向外的异常横向运动情况。

3 匣钵运动仿真分析

本工作在Abaqus仿真软件中分别模拟辊道-匣钵传动简化模型在单独重力场和热力耦合场环境下匣钵的运动状态并分析辊棒的受力情况和挠度。仿真平台使用惠普笔记本电脑，型号为HP ZBook Power 15.6 inch G8 Mobile Workstation PC，CPU型号为i9-11900H，主频为2.5 GHz，内存为32 GB。
因此，受仿真条件限制，为减少计算时间，匣钵在辊道上的运动仿真长度为1000 mm。在生产过程中，匣钵在辊道上运行速度较慢，且匣钵横向运动与辊棒转速无直接联系，因此本工作将辊棒的旋转角速度提高，使得匣钵通过1000 mm辊道所需的时间为5 s。对于其他的仿真参数，接触方式采用通用接触，摩擦系数设置为0.1，接触切向行为设为罚摩擦公式，匣钵间初始间隙设为1 mm，最后进行网格划分，节点总数为56742，单元总数为46944。

3.1　重力场匣钵运动仿真分析

在研究单一重力场对匣钵运动特性的影响时，为了避免重力加速度突变使得仿真结果异常，在仿真过程中，通过将重力加速度设为在1 s的时间内逐渐从0 m/s²增加到9.8 m/s²的方法来避免仿真结果出现异常情况。
设置好相关参数后开始仿真，在重力方向— Y 施加重力加速度后，辊棒开始出现弹性形变，匣钵跟随辊棒形变方向下沉，其匣钵中点 Y 向位移-时间曲线如图6所示。在0～1 s时间段内，由于重力加速度持续增大，匣钵中点下沉迅速，下沉至-1.27 mm左右。在1～5 s时间段内，匣钵中点 Y 向位移-时间曲线呈波浪状，此时匣钵中点在 Y 向位移变化较为稳定。经分析，波浪状曲线出现的主要原因是匣钵在运动过程中，每根辊棒承受的载荷不均匀，重力场下辊棒的应力分布情况如图7所示。

图6     重力场匣钵中点 Y 向位移-时间曲线

图7     重力场下辊棒的应力分布
在仿真过程中，选取四个匣钵的正面上边缘中点BOX-1、BOX-2、BOX-3、BOX-4四点作为位移监测点，如图8所示。仿真结束后得到4个监测点的 Z 向位移-时间曲线，如图9所示。图中纵坐标表示匣钵在 Z 轴方向的位移量，四条曲线总体呈外扩样式，并且随着仿真的进行，外扩行为的变化较为稳定。经分析，其原因是辊棒达到一定的弹性变形量时，匣钵一方面受式(1)中 的影响往辊道中间靠拢，另一方面受式(6)中 的作用往外位移，两者相对平衡使得匣钵的横向位移的变化相对稳定。仿真结束后，由图9计算可得5 s时匣钵最大间隙由初始间隙1 mm增大至1.71 mm。

图8     匣钵监测点分布

图9     重力场 4 个匣钵监测点的 Z 向位移-时间曲线

3.2　热力耦合场匣钵运动仿真分析

为分析热力耦合场对匣钵横向异常位移的影响，本工作在上述匣钵运动仿真模型中耦合重力场和温度场。为保证仿真结果的可靠性，本工作选取匣钵在烧结炉中烧结时的最大温度1150 ℃作为仿真温度，同样将重力加速度设为在1 s的时间内逐渐从0 m/s²增加到9.8 m/s²。仿真结束后，匣钵中点 Y 向位移-时间曲线如图10所示。在0～1 s时间段内，由于重力加速度持续增大，匣钵中点下沉迅速，下沉至-1.37 mm左右。在1～5 s时间段内，匣钵中点 Y 向位移-时间曲线为波浪状，与单独重力场下的曲线变化类似，热力耦合场下的辊棒应力分布情况如图11所示。

图10     热力耦合场匣钵中点 Y 向位移-时间曲线

图11     热力耦合场下辊棒的应力分布
在仿真过程中，监测点的选取与单重力场情况下一致。仿真结束后，得到4个匣钵监测点的 Z 向位移-时间曲线，如图12所示。曲线仍呈外扩形式，随着时间的推移，外扩行为的变化同样较为稳定。热力耦合场下5 s时匣钵最大间隙由初始间隙1 mm增大至1.55 mm，与单独重力场相比5 s时匣钵最大间隙减小了0.16 mm。

图12     热力耦合场 4 个匣钵监测点 Z 向位移 - 时间曲线

3.3　单重力场与热力耦合场对比分析

本工作通过对比单独重力场与热力耦合场下的匣钵中点 Y 向位移-时间曲线以及4个匣钵监测点Z向位移-时间曲线，分析热力耦合引起的辊道的弹性变形对匣钵异常横向运动的影响程度。匣钵中点 Y 向位移-时间曲线的对比如图13所示，BOX- Y -Middle代表单独重力场下的匣钵中点 Y 向位移-时间曲线，BOX- Y ’-Middle代表热力耦合场下的匣钵中点 Y 向位移-时间曲线。从图13中可看出，由于热力耦合场比单独重力场多施加了温度场，导致匣钵中点 Y 向位移程度变大，匣钵下沉的速度更快，从而说明热力耦合场使得辊棒的变形加剧。

图13     匣钵中点 Y 向位移 - 时间曲线对比
4个匣钵监测点 Z 向位移-时间曲线的对比如图14所示，BOX-1- Z 代表单独重力场下第一个匣钵监测点 Z 向位移-时间曲线，BOX-2- Z 、BOX-3- Z 、BOX-4- Z 同理。BOX-1’- Z 代表热力耦合场下第一个匣钵监测点的 Z 向位移-时间曲线，BOX-2’- Z 、BOX-3’- Z 、BOX-4’- Z 同理。

图14     4 个匣钵监测点 Z 向位移 - 时间曲线对比
在0～1 s时间段内，热力耦合场下匣钵的 Z 向位移大于重力场下匣钵的 Z 向位移。经分析，主要原因是匣钵在热力耦合场下比在单独重力场下的 Y 向位移和下沉速度更大，导致旋转力矩 更大，从而使得热力耦合场下的匣钵 Z 向位移增大。
在1～5 s时间段内，热力耦合场下匣钵的 Z 向位移小于重力场下匣钵的 Z 向位移。为了更好地分析其原因，本工作结合了1.5 s时的辊棒挠度对比曲线，如图15所示。图中Roller-1表示重力场下第1根辊棒，Roller-2、Roller-3、Roller-4以此类推。Roller-1’表示热力耦合场下第1根辊棒，Roller-2’、Roller-3’、Roller-4’以此类推。图中Roller-4和Roller-4’辊棒在1.5s时由于未承载匣钵，其曲线弧度较小。从图15中可以看出辊棒在热力耦合场下的弹性变形量比在单独重力场下的弹性变形量更大。由于此时匣钵在 Y 向的位置相对稳定，不再沿 Y 轴快速下沉，相对0～1 s时间段而言，热力耦合场和重力场下匣钵所受旋转力矩的差值减小。而通过第2节匣钵受力分析可知，当辊棒的弹性变形量达到一定的程度，匣钵反而会克服最大静摩擦力向内滑动，阻止匣钵因受旋转力矩向外运动。此阶段热力耦合场下辊棒的弹性变形量更大，匣钵克服最大静摩擦力向内滑动的位移量增加，导致在图14中1～5 s时间段内，热力耦合场与重力场相比，匣钵向外的 Z 向位移量减小。综上所述，在辊道窑烧结过程中，由于温度和重力耦合所引起的辊棒弹性变形导致匣钵在辊道上产生了异常横向位移。

图15     1.5s 时辊棒挠度曲线对比

4 匣钵传动机构设计及改进验证

4.1　匣钵传动机构设计

由上述研究可知，匣钵在辊道上运行时，由于热力耦合引起了辊棒的弹性变形，从而导致匣钵出现异常横向位移。为了降低此现象的影响程度，辊棒的弹性变形量可简化为挠度概念，考虑到辊棒工作环境为高温状态，辊棒挠度虽与其自身截面尺寸等有关，但为了避免辊棒过大导致辊道窑舱体气密性差及加工困难等情况的出现，本工作通过合理设计辊棒支撑机构以减小辊棒挠度，从而减小辊棒的弹性变形量，并且通过在辊道窑中设置匣钵夹紧装置和匣钵挡齐装置进一步减小匣钵的横向位移量。
本工作辊棒支撑机构如图16所示，采取在辊棒下方安装4个支撑机构来减小辊棒的弹性变形。支撑机构通过支撑机构垫板与辊道窑舱体连接，为便于支撑机构的安装，采取L形支撑板作为主体结构。支撑板通过6个螺钉起到固定作用，同时在支撑板上安装支撑机构轴承，以便后续调整辊道的平面度和减小辊棒在支撑机构上滚动时的摩擦力。

图16     辊棒支撑机构图
匣钵夹紧装置如图17所示，舱体外侧通过安装板左右对称安装夹紧装置。匣钵夹紧装置的主要驱动元件为夹紧气缸，夹紧气缸活塞杆通过螺栓与对夹机构连接板连接，对夹机构连接板与挡板导柱固定连接，挡板导柱的另一端连接夹紧板垫板。匣钵夹紧装置工作时，夹紧气缸动作，带动挡板导柱和夹紧板垫板朝舱体内部运动，完成夹紧动作。

图17     匣钵夹紧装置图
为避免多个匣钵夹紧时匣钵之间相互错开，通过增加匣钵挡齐装置使匣钵在同一直线上，匣钵挡齐装置如图18所示。采用气缸作为其驱动装置，气缸固定支座安装在舱体外侧机架上，同时挡齐气缸通过螺栓与气缸固定支座连接，挡齐气缸伸缩端与旋转块连接，旋转块与旋转轴固定连接。当匣钵运行到靠近旋转轴位置时，挡齐气缸开始动作，带动旋转块旋转，从而带动旋转轴和挡板固定条旋转，挡板固定条挡住多个匣钵使其在同一直线上，完成匣钵挡齐动作。

图18     匣钵挡齐装置图

4.2　改进验证

由于在辊道窑完成烧结后，匣钵间隙需要停产才可测量，同时使用直接测量匣钵间隙方法的测量精度较差。因此，考虑到生产经济效益、验证成本以及验证结果的可信度，本工作以某工厂一个月的匣钵破损率来验证安装了辊棒支撑机构、匣钵夹紧装置及匣钵挡齐装置后匣钵在辊道上的横向位移是否满足生产要求。
经生产验证，该工厂在安装改进机构前的一个月内，锂离子电池正极材料烧结工位中的匣钵破损率为1.53%，而在安装改进机构后的一个月内，锂离子电池正极材料烧结工位中的匣钵破损率为0.69%。结果表明通过合理设计辊棒支撑机构、匣钵夹紧装置及匣钵挡齐装置可以大幅度改善匣钵在辊道窑烧结过程中的异常横向运动，减小匣钵破损率。因此，该匣钵传动机构设计合理有效，能较好地满足生产要求。

5 结 论

本工作分析了辊道窑烧结过程中匣钵异常横向运动的原因，并通过建立辊道-匣钵传动简化模型，对在热力耦合场下的辊棒受力和匣钵位移进行有限元分析。最后，本工作设计了一种匣钵传动机构，进行了生产对比验证，得到如下结论：
（1）在辊道窑烧结过程中热力耦合引起的辊棒的弹性变形是导致匣钵异常横向运动的主要原因。
（2）热力耦合场相比于重力场引起的辊棒弹性变形更大，导致热力耦合场下匣钵的异常横向位移会先因匣钵下沉速度更快而相对增大，接着又因辊棒挠度增加使得匣钵向中部滑动的趋势更大而相对减小。
（3）通过设计匣钵传动机构，包括辊棒支撑机构、匣钵夹紧装置和匣钵挡齐装置，能大幅度减轻匣钵的异常横向运动现象。生产对比验证的结果表明设计的匣钵传动机构可有效降低匣钵破损率，能较好地满足生产要求。
匣钵在辊道窑内的运动优化虽能减少匣钵碰撞现象，提高匣钵使用次数，但由于辊道窑内环境较为复杂，仍可能因意外发生碰撞事故。因此可从匣钵材料和窑内环境控制等方面进一步减少匣钵的破损。目前，对匣钵材料的煅烧性能和腐蚀机理的研究较多，但同时也应考虑匣钵材料的强度、轻量化和抗疲劳性，选取合适的匣钵材料，提高匣钵的使用寿命。烧结温度对锂离子电池正极材料的结构和电化学性能有重要的影响，但也应考虑其对烧结设备的影响，保证产品质量的同时，尽量降低辊道窑温度并控制辊道窑内温度均匀，减少烧结设备的变形断裂现象以及匣钵运动时产生的意外碰撞事故。

通讯作者： 彭可，教授，主要研究方向为机械制造自动化与机电一体化控制，E-mail： 。

第一作者： 彭可（1973—），男，博士，教授，主要研究方向为机械制造自动化与机电一体化控制，E-mail：；

本刊推荐

（点击图片可跳转到相应文章合集）

邮发代号：80-732
联系热线：010-64519601/9602/9643
投稿网址：http://esst.cip.com.cn/CN/2095-4239/home.shtml

下一条：返回列表
版权：如无特殊注明，文章转载自网络，侵权请联系cnmhg168#163.com删除！文件均为网友上传，仅供研究和学习使用，务必24小时内删除。
相关推荐
热门信息