有色金属科学与工程  2021, Vol. 12 Issue (4): 82-87, 125
文章快速检索     高级检索
铝箔剪切机碎屑收集管优化设计[PDF全文]
刘跃1,2 , 常玲玲1 , 李会荣1 , 管小荣2     
1. 陕西国防工业职业技术学院智能制造学院,西安 710300;
2. 南京理工大学机械工程学院,南京 210094
摘要:传统铝箔剪切机碎屑收集装置存在2支管吸力不均匀问题,严重影响剪切机工作效率。文中使用计算流体力学方法(Computational Fluid Dynamics,CFD)对管路挡板结构进行了优化改进,在深入分析挡板布置方式对支管流态影响规律基础上,给出管路优化建议。数值结果表明:2个挡板高度相等时,第2个挡板会对第1个支管后部气流产生"阻碍"作用,形成低速回流区,无法形成负压及吸力;随着第2个挡板高度降低,第2个支管内负压及速度值降低,第1个支管内负压及速度值增大;当挡板比例为1/2时,2个支管物理量趋于相等,支管口速度相差约3%,可满足工业需求。此外,改变挡板高度及间距并保持挡板高度比例为1/2时,2个支管速度相差约5%,说明研究结果具有一定工程适用性。
关键词铝箔剪切机    铝箔碎屑    结构优化    负压    数值模拟    
Optimal design of scrap collection tube for aluminum foil shearing machine
LIU Yue1,2 , CHANG Lingling1 , LI Huirong1 , GUAN Xiaorong2     
1. School of Intelligent Manufacturing, Shaanxi Institute of Technology, Xi'an 710300, China;
2. School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
Abstract: The scrap collection device of the traditional aluminum foil shearing machine has the problem of uneven suction of the two pipes, which seriously affects the working efficiency.Applying Computational Fluid Dynamics (CFD), we optimized and improved the pipeline baffle structure upon an in-depth analysis on the effect of baffle layout on the branch pipe flow pattern. Suggestion for piping optimization are put forward. Numerical results show: when the heights of the two baffles are equal, the second baffle will "obstruct" the air flow behind the first branch pipe, forming a low-velocity recirculation zone. Thus, negative pressure and suction fail to be formed. The decreasing height of the second baffle reduces the negative pressure and velocity of the second branch pipe, which increases the negative pressure and velocity in the first branch pipe. When the baffle ratio is 1/2, the physical quantities of the two branch pipes tend to be equal, and the difference in the speed of the branch pipe mouth is about 3%. The industrial needs can be met. When changing the baffle height and spacing and keeping the baffle height ratio at 1/2, the speed difference between the two branch pipes is about 5%. It shows that the research results have certain engineering applicability.
Keywords: aluminum foil shearing machine    aluminum foil scraps    structure optimization    negative pressure    numerical simulation    

铝箔由于具有防潮、气密、耐磨蚀、无毒无味等优点,广泛应用于食品、医疗、电子等行业[1-4]中,我国铝箔产量更是在2008年开始便居全球首位,市场占比达56.7%[5]。经过长足发展,我国铝箔生产工艺、生产设备已逐步系统化、完善化,形成了专门的轧制、涂层、开卷、分切等模块[6]。其中,铝箔剪切机是在分切模块中将铝箔分剪成不同尺寸以适应不同行业需求的关键设备。查阅文献[7-9]可知,相关学者对剪切机关键零部件如碟形刀、套筒夹持装置等进行了持续优化,这些成果对改善铝箔生产质量起到了良好的促进作用。但在观察铝箔实际剪切流程时可以发现,剪切机工作时会不可避免地产生碎屑,其长时间运行时由于碎屑累积明显影响铝箔分切质量,目前大多数剪切机需要按时停机进行碎屑清扫,这对剪切机工作效率产生较大影响,现有针对铝箔碎屑收集的优化设计成果较少[10],且存在结构复杂、成本较高的缺陷。

负压管内流动是常见的介质输送通道[11-13]。本文拟根据铝箔剪切机工作特点设计一种由双支管组成的铝箔碎屑收集管路系统,由鼓风机向管内输送气流,通过管内局部尺寸变化产生负压[14-15]。采用目前工程中常用的计算流体方法(CFD)[16-18]对不同挡板设计时的压力-速度耦合流场进行三维数值计算,在分析挡板结构改变管内压力-速度分布规律及流动现象产生原因基础上,重点探讨了获得支路均匀负压及吸力的挡板布置方法,并对设计方法的可靠性进行分析,为铝箔剪切机设备改进提供参考。

1 几何模型

图 1所示为常用的LT1350剪切机分切铝箔流程[5],主要由2个碟形刀完成铝箔切割工作。为了完成铝箔碎屑收集工作,初步设计的圆形管路样品如图 2所示,管路外观由主管及2个支管组成,支管口置于2个碟形刀附近,预想通过管内气流流动产生负压并在支管口形成吸力使得碎屑通过支管进入主管,并通过主管出口统一收集。管路气流输送动力来源于进口处鼓风机,铝箔碎屑厚度为0.1~0.2 mm,比较轻薄,为了获得支管口理想吸力,气流输送功率参考常用工业吸尘器功率进行设置[19],取2.2 kW,对应气流输送量为320 m3/h。

图 1 铝箔剪切机碟形刀 Fig. 1 Foil slitter

图 2 铝箔碎屑收集圆管 Fig. 2 Foil scrap collection pipe with round model

2 计算模型 2.1 计算方程

输送介质为气体,整体流动马赫数Ma < 0.2,主要求解黏性不可压流体方程组,其矢量形式如下[20]

(1)

式(1)中:ρ为密度,U为速度矢量,f为体积力,p为压力,μ为气体动力黏度。该方程可较好地反映管流黏性切应力分布及流动压力损失。

此外,管流雷诺数Re=1.0×105(由入口尺寸及速度确定),属于湍流运动,经比较,选择在近壁及远场均有较好表现的SST湍流模型(Shear-stress transport)。其具体方程构造如下[21]

(2)
(3)

式(2)、式(3)中:模型常数β*=0.09,卡门系数k=0.41。此外,混合函数F1、湍动能生成项Pk、动力黏度μt运动黏度νt及混合模型参数αβσkσω定义及取值详见文献[22]。

2.2 数值模型

考虑到方管对比于圆管具有加工工艺简单、数值模型构造精度高的优点,本文中管路采用方管设计。以设置2个等高度挡板且考虑挡板厚度为例,三维方形管路外形如图 3上部所示,主管边长为90 mm,支管截面尺寸为90 mm×40 mm,支管标准间距为1 000 mm。三维坐标系如图 3所示,管路延轴向往后为x轴正方向,纵向为y方向,横向为z方向。管路进口处为速度进口,出口为压力出口,2个支管P1、P2为压力进口。网格划分如图 3下部所示,整体为结构网格,为了充分考虑壁面作用对气流流态的影响,单独绘制壁面层网格,并对较关心的支管处网格进行加密处理。此外,文中不同结构管路计算网格均相同,以保证计算结果的可比性。使用Fluent软件求解不可压流体N-S方程组,方程离散采用二阶迎风格式[23]

图 3 数值模型及局部网格 Fig. 3 Foil scrap collection pipe and numerical grid

2.3 模型验证

为了验证文中所选数值边界及湍流模型的适用性,课题组设计加工出等高度挡板方形管路样品(具体模型及尺寸见3.1节),并使用转子流量计测量了2个支管的气流流量,经过换算可得P1、P2支管纵向速度v的测量值Exp.(v=Q/AQ为气流流量,A为截面面积),约定v的方向为流入管内为正,流出管外为负。2个支管口平均纵向速度的模拟值与实验值对比如图 4所示。模拟结果具体为(3.29 m/s,-8.17 m/s),实验结果具体为(3.13 m/s,-7.78 m/s),通过比较可以发现数值结果与测量值误差相差在5%左右,这说明了本文中使用数值方法是可行的。为了便于比较,将等高度挡板计算条件下支管口平均纵向速度值及测量值列于表 1中。

图 4 数值边界及湍流模型验证 Fig. 4 Validation of numerical boundary and turbulence model

表 1 不同挡板高度比例(h2/h1)时支管纵向速度比较/(m·s-1) Table 1 Longitudinal velocity comparison of different pipe structures/(m·s-1)
点击放大

3 结果分析

文中为准定常计算,为了获得收敛的分析数据,方程残差标准设置为10-5,并且监测点数据(P1P2出口中心点压力值)不再变化时认为计算结束。

3.1 等高度挡板计算

决定支管吸力大小的负压分布取决于管内局部压力-速度耦合变化结果,为了探讨双支管气流流态特征及获得管内合理负压分布的挡板设置方法,文中首先计算设置2个等高度挡板(第2个挡板高度与第1个挡板高度之比h2/h1=1)时的管路流场,并将该算例命名为"等高度挡板"计算。等高度挡板管路设计如图 5所示,2个挡板高度为45 mm。计算开始前需要设置初值条件,根据鼓风机流量及进口边界确定初始风速为U(10.97,0,0),单位为m/s,P1P2支管口及出口压力初值为0。

图 5 挡板等高管路模型 Fig. 5 Two baffles with same height

图 6所示为等高度挡板计算中心截面的压力分布,可以看到第2个支管附近出现了预想的负压区域,气流流向管内,并可在支管口产生吸力;但在第1个支管附近却没有出现理想的负压,压力值大于0,且2个挡板之间也出现了较大的正压区域,这便导致第1个支管不能产生吸力,进而无法满足铝箔碎屑收集要求。图 7给出了等高度挡板计算时2个支管附近的流线及纵向速度分布。可以看到与压力分布相对应的结果:P2支管纵向速度为负,平均速度为-8.17 m/s,气流运动较为顺畅,仅在支管与主管连接处由于方向改变产生低速涡旋区;P1支管口平均速度为3.29 m/s,气流流出管外,可以看到,2个挡板之间形成了较大的反向涡流区域,导致P1支管附近流态较复杂,这说明第2个挡板会对第1个挡板后部区域气流产生"阻碍"作用,这是导致P1支管无法产生吸力的主要原因。

图 6 压力云图 Fig. 6 Pressure contours

图 7 支管附近流线及纵向速度 Fig. 7 Longitudinal velocity and local streamline near branch pipes

3.2 挡板优化设计

由以上分析可知简单设置2个等高度挡板无法满足铝箔碎屑收集要求,考虑到第2个挡板对第1个挡板后部气流的"阻碍"影响,一种简单的办法是尝试改变第2个挡板高度来增大气流过流面积,进而平衡2个支管负压分布及吸力水平。为了验证该思路的可行性,计算了不同挡板高度比例(h2/h1)时的管内流场。图 8所示为挡板高度比例分别为2/3和1/2时的方管模型,第1个挡板保持45 mm高度不变,第2个挡板高度分别设置为30 mm和22.5 mm,管路中心截面压力分布云图如图 9所示,相比于图 5,降低第2个挡板高度后,第1个支管附近压力值呈下降趋势并出现了需要的负压,进一步观察可发现挡板高度比例为2/3时,第1个支管附近负压值仍小于第2个支管附近负压值,而1/2比例时2个支管负压观察值较相近。图 10给出了1/2比例时支管附近流线及纵向速度分布,可以看到第1个支管附近流态得到明显改善,不利回流区范围减小,气流流向管内且2个支管纵向速度分别为-6.18,-6.01 m/s,相差不到3%。

图 8 不同挡板高度比例时管路模型 Fig. 8 Pipe models under different baffle proportions

图 9 压力云图 Fig. 9 Pressure contours

图 10 支管附近流线及纵向速度 Fig. 10 Longitudinal velocity and local streamline near branch pipes

图 11给出了不同挡板高度比例时特征线上(坐标位置(0,40,0)~(1 800,40,0),单位为mm)压力和速度量化值,可以看到挡板高度比例为1时,2个挡板之间均为正压分布,对应第1个支管纵向速度为正。随着挡板比例降低(从2/3至1/3),第1个支管附近负压区域增加,纵向速度负值逐渐增加,产生吸力;与之对应的是第2个支管附近负压区域减小,速度负值减小。可见随着挡板比例下降,2个支管附近负压区域及纵向速度差值逐渐减小,并趋于均匀。同时,由图 11(b)速度量化值也可以看到,当第2个挡板高度与第1个挡板高度为1/2时,对应支管附近特征线上纵向速度值较为接近,这与图 10结果也是一致的。

图 11 特征线上的压力及速度沿程分布 Fig. 11 Pressure and velocity comparison along typical line

由以上分析可知,通过改变挡板高度比例来实现2个支管产生均匀吸力是可行的,为了便于设计及生产,根据计算结果,推荐h2/h1=1/2比例设计。同样为了便于比较,不同挡板高度比例时支管纵向速度值列于表 1中。

3.3 设计可靠性分析

为了探讨随设计参数变化时,挡板h2/h1=1/2条件下支管口速度均匀的一致性,比较了保持挡板高度比例为1/2,改变挡板高度及减小挡板距离时的管路计算结果。图 12给出了第1个挡板高度为40 mm,第2个挡板高度为20 mm时的中心截面压力及速度分布云图,可以看到2个支管负压及纵向速度均较相近,速度值分别为-5.54,-5.24 m/s;图 13给出了2个挡板距离缩小为500 mm时中心截面的压力及速度分布云图,2个支管负压及纵向速度同样较相近,速度值分别为-7.78,-7.39 m/s。

图 12 变挡板高度管路压力及纵向速度分布 Fig. 12 Pressure and velocity comparison between different baffle heights

图 13 变挡板距离管路压力及纵向速度分布 Fig. 13 Pressure and velocity comparison between different baffle distances

图 14给出了变挡板高度及变挡板距离时特征线上压力和速度量化值,可以看到计算规律与前文一致,2个支管附近均出现负压区域且纵向速度值非常接近。

图 14 特征线上的压力及速度沿程分布 Fig. 14 Pressure and velocity comparison along typical line

由以上分析可知,保持挡板高度1/2比例不变,改变管路结构尺寸设计时2个支管速度差值在5%左右,可满足双支管吸力的均匀性,进而满足铝箔剪切机碎屑收集要求。

4 结论

使用CFD方法探讨了利用负压管完成铝箔剪切机碎屑收集任务的可行性,在分析双支管流态特征的基础上,重点探讨了如何通过简易改变挡板设置获得2个支管均匀吸力。通过本文数值结果主要得出以下结论:

1) 2个挡板高度相等时,第2个支管附近可产生负压及吸力,但第1个支管附近气流由于第2个挡板的"阻碍"作用形成较大的低速涡旋回流区,气流流出管外,无法形成吸力。

2) 第1个挡板高度不变,随着第2个挡板高度降低,第1个支管附近负压值增加,气流流向管内速度及吸力增大;第2个支管附近负压值减小,气流流向管内速度及吸力减小。当比例为1/2时,2个支管物理量趋于平衡分布,速度值相差3%左右,可满足工业要求。

3) 固定第2个挡板与第1个挡板高度比例为1/2,改变挡板高度值及间距时,2个支管口向内速度差值在5%左右,说明挡板1/2比例随管路结构参数影响较小,可推广使用。

参考文献
[1]
王伟东. "禁塑令"为包装用铝带来发展新机遇[J]. 中国金属通报, 2020(6): 20–21. DOI: 10.3969/j.issn.1672-1667.2020.06.013.
[2]
王盈来, 李艳红, 黄燕山. 铝箔涂碳层厚度对锂离子电池性能的影响[J]. 电池, 2021, 51(1): 50–53. DOI: 10.3969/j.issn.1008-7923.2021.01.011.
[3]
PANIN S, BURKOV M, LYUBUTIN P, et al. Application of aluminum foil for "strain sensing" at fatigue damage evaluation of carbon fiber composite[J]. Science China: Physics, Mechanics and Astronomy, 2014, 57(1): 59–64. DOI: 10.1007/s11433-013-5368-y.
[4]
PANGESTI G G, PANDIANGAN K D, SIMANJUNTAK W, et al. Synthesis of zeolite-Y from rice husk silica and food grade aluminum foil using modified hydrothermal method[J]. Journal of Physics: Conference Series, 2021, 1751(1): 012089. DOI: 10.1088/1742-6596/1751/1/012089.
[5]
张宇树, 王祝堂. 中国已成为世界铝箔初级强国[J]. 轻合金加工技术, 2018, 46(12): 6–9.
[6]
辛达夫. 当代铝箔生产工艺及装备[M]. 北京: 冶金工业出版社, 2013: 206-211.
[7]
李会荣, 张永军. 铝箔分切机碟形刀结构优化[J]. 机械设计与制造工程, 2017, 46(2): 74–76. DOI: 10.3969/j.issn.2095-509X.2017.02.015.
[8]
江开伦, 商雨. 1700mm铝箔分切机剪切装置的改进[J]. 有色金属加工, 2020, 49(3): 47–49. DOI: 10.3969/j.issn.1671-6795.2020.03.012.
[9]
田磊. 高强度套筒自动送料及成型执行系统关键技术研究[D]. 石家庄: 河北科技大学, 2018.
[10]
郑晓杰, 胡海庆, 时超凡. 铝箔废边剪切机吸收打碎装置[P]. 中国: CN209682296U, 2019-11-26.
[11]
胡盛翔. 基于Fluent的气力输送弯管流场仿真模拟[J]. 水泥工程, 2020(增刊1): 34–36.
[12]
严天曈, 谢果, 孙立成, 等. 低温热源驱动两级负压管式淡化装置性能研究[J]. 太阳能学报, 2020, 41(6): 357–362.
[13]
张博文, 王海宁, 张迎宾. 公路隧道空气幕与射流风机通风方法对比及优化[J]. 有色金属科学与工程, 2020, 11(3): 80–89.
[14]
STUDENT O Z, SVIRS'KA L M, DZIOBA I R. Influence of the long-term operation of 12Kh1M1F steel from different zones of a bend of steam pipeline of a thermal power plant on its mechanical characteristics[J]. Materials Science, 2012, 48(2): 239–246. DOI: 10.1007/s11003-012-9498-6.
[15]
祁金胜, 曹洪振, 石岩, 等. 虾米腰弯管内置导流板优化[J]. 山东大学学报(工学版), 2020, 50(5): 64–69.
[16]
邓玉伟, 潘作峰, 侯杭生. 简易风道气动噪声仿真分析[J]. 机械设计, 2020, 37(增刊2): 176–180.
[17]
刘志林, 何凯, 文继有, 等. 大口径薄壁白铜管无外模扩径拉拔数值模拟[J]. 有色金属科学与工程, 2021, 12(1): 90–98.
[18]
郭学益, 闫书阳, 王双, 等. 数值模拟氧气底吹熔炼工艺参数优化[J]. 有色金属科学与工程, 2017, 8(5): 21–25.
[19]
明乐乐, 胡小文, 张辉, 等. 吸尘器用高速离心风机的优化设计[J]. 风机技术, 2018, 60(增刊1): 7–13.
[20]
林建忠, 阮晓东, 陈邦国, 等. 流体力学[M]. 2版. 北京: 清华大学出版社, 2019: 72-108.
[21]
MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8): 1598–1605. DOI: 10.2514/3.12149.
[22]
LIU Y, GUAN X R, XU C. A production limiter study of SST-SAS turbulence model for bluff body flows[J]. Journal of Wind Engineering & Industrial Aerodynamics, 2017, 170: 162–178.
[23]
吴子牛. 计算流体力学基本原理[M]. 北京: 科学出版社, 2003.