2. 青岛软控机电工程有限公司, 青岛 266300
2. Qingdao Soft Control Mechanical and Electrical Engineering Co., Ltd., Qingdao 266300, China
冷却鼓是对压延机压延成型胶片进行降温、冷却定形的组件。混炼胶片从冷喂料挤出机挤出,经过胶片输送机被传送至压延机挤压成型,通过压延成型的胶片脱离压延机的鼓筒后,胶片的温度还很高,所以需经过一组冷却鼓为胶片充分降温冷却,使其定形,消除成型过程中胶片的内应力,然后才能把胶片卷曲成捆。使用冷却鼓对胶片进行冷却的工作原理是向冷却鼓内部通入冷却水,冷却水通过冷却鼓与胶片传热,对胶片进行冷却。
在连续大批量生产中,胶片温度的控制是影响生产效率及产品质量的重要因素。影响胶片温度的参数有很多,包括冷却水流量、胶片厚度、宽度、移动速度、冷却鼓结构、数量、材料种类和表面性能等。近年来很多学者对冷却鼓或类似结构的设计制造及冷却效果进行了研究。杨云等[1]介绍了冷却辊加工的技术关键和工艺措施。陈洪冰等[2]根据冷却辊的使用工况,对自行研制的材料喷焊层的化学成分、硬度分布、相结构与组织结构进行了分析与测试。杨宏伟等[3]将冷却牵引装置冷却辊由6个增加到10个,并对压延机后端张力控制装置进行改进,改进后不但解决了产能瓶颈,且张力控制稳定,提升了产品质量。韩亚娜等[4]从冷却辊筒的结构、材料性能等方面出发,分析了冷却辊筒的机械加工工艺及加工过程中易出现的问题,并提出解决措施。郭茜等[5]采用有限元分析方法对冷却鼓进行了稳态温度场数值模拟,获得了冷却鼓温度场及其热流的分布特点。陈北荣[6]通过分析冷却辊的故障,对冷却辊进行了局部改进设计,大大提高了某凹印机冷却辊的工作可靠性。张羽玲等[7]利用计算流体动力学(computational fluid dynamics,CFD)方法对冷却辊表面、轴向和周向的温度场分布进行仿真分析,研究了喷淋式冷却辊和传统式冷却辊稳定工作状态下的温度分布散点图及云图。
尽管冷却鼓已经广泛应用于胶片冷却,但目前其设计主要还是依据经验,缺乏完善的传热分析和计算方法。由于胶片和冷却鼓接触面不断变化,属于固体移动边界传热问题,目前难以进行有效的数值模拟,为此本文将固体胶片假设为流体胶片,以流体胶片的流动来代替固体胶片递进移动,通过热-流-固耦合换热计算,模拟胶片与冷却鼓之间的换热,并与实测数据进行比较,以验证本文传热模拟方法的合理性。
1 胶片-冷却鼓装置及冷却过程分析 1.1 胶片-冷却鼓装置及几何模型本文研究的胶片冷却装置包括冷却架、驱动电机和多个冷却鼓。冷却鼓以上下双层的结构安装在冷却架上,驱动电机通过带轮和皮带带动多个冷却鼓共同旋转,并与胶片移动速度保持一致。胶片与冷却鼓的三维模型如图 1所示,为保证胶片向前移动,上下层的冷却鼓旋转方向相反,若鼓1为左旋,则鼓2为右旋,其他鼓的旋向依次类推。
冷却鼓结构包括一个双层钢制筒体及转轴、进水/出水管道及进水、出水口,转动机构及冷却鼓端板,其中通道的进水与出水管道在其转轴中,冷却鼓内层上有带有转轴通孔的纵向圆形隔板以支撑其鼓筒体不变形,冷却鼓内部为夹套结构,即由内外鼓构成,内鼓筒是空心的,其上焊接螺旋式夹板,夹板的法向截面形状为矩形,外鼓、内鼓、夹板之间形成冷却水通道。冷却鼓内部结构如图 2所示,结构参数见表 1。冷却鼓中螺旋流道增加了冷却水在鼓中的流通路径和停留时间,目的是提高冷却效率,节约冷却水。
在实际运行中,胶片通过冷却鼓与冷却水进行换热,开始时,胶片、冷却鼓和冷却水处于非稳态阶段,温度一直发生改变。当运行一段时间后,胶片、冷却鼓和冷却水三者之间换热达到动态热平衡,即冷却鼓的温度基本不再发生变化,胶片放出的热量等于冷却水吸收的热量,达到稳定状态。本文采用稳态计算, 即
1) 连续性方程
$ {{\rm div}(\mathit{\boldsymbol{U}}) = 0} $ | (1) |
2) 动量方程
$ \left\{ \begin{array}{l} {\rm div}(u\mathit{\boldsymbol{U}}) = \upsilon {\rm div}({\rm grad}\ u) - \frac{1}{\rho } \cdot \frac{{\partial p}}{{\partial x}} + {S_u}\\ {\rm div}(v\mathit{\boldsymbol{U}}) = \upsilon {\rm div}({\rm grad}\ v) - \frac{1}{\rho } \cdot \frac{{\partial p}}{{\partial y}} + {S_v}\\ {\rm div}(w\mathit{\boldsymbol{U}}) = \upsilon {\rm div}({\rm grad}\ w) - \frac{1}{\rho } \cdot \frac{{\partial p}}{{\partial z}} + {S_w} \end{array} \right. $ | (2) |
3) 能量方程
$ {\rm div}(T\mathit{\boldsymbol{U}}) = {\rm div}\left( {\frac{\lambda }{{\rho {c_p}}}{\rm grad}\ T} \right) + \frac{{{S_h} + \varphi }}{{\rho {c_p}}} $ | (3) |
式中,
4) 在冷却鼓上胶片与空气接触表面进行对流换热,采用第三类边界条件,其换热方程[9]为
$ {\left. { - \lambda \frac{{\partial T}}{{\partial n}}} \right|_{n = b}} = {h_1}({T_{n = b}} - {T_{\rm air,f}}) $ | (4) |
式中,T为胶片的温度, K;n为胶片沿厚度方向上的坐标, mm;b为胶片的厚度, mm;h1为胶片与空气的对流换热系数, W/(m2·K);Tn=b为与空气接触的胶片表面温度, K;Tair, f为靠近胶片的空气温度, K。
5) 冷却鼓与空气接触表面进行对流换热,其换热方程采用牛顿冷却公式[9]
$ {q_2} = {h_2}({T_{\rm out}} - {T_{\rm air,c}}) $ | (5) |
式中,q2为冷却鼓壁面与空气进行对流的换热量, W/m2;h2为冷却鼓壁面与空气的对流换热系数, W/(m2·K);Tout为与空气接触的冷却鼓壁面温度, K;Tair,c为靠近冷却鼓壁面的空气温度, K。
1.3 胶片冷却过程的技术条件本文中有两种胶片规格参数,宽度分别为680 mm和760 mm,厚度分别为5.15 mm和6.4 mm,分别对这两种规格的胶片冷却过程进行数值模拟,共有6种工况,具体参数见表 2,本文以工况1为例进行分析。
在整个冷却过程中, 胶片的几何形状变化很小, 可忽略因牵引力而引起的胶片在拉伸方向上的形变; 由于摩擦引起的胶片内能的变化相对比较小, 因此忽略冷却鼓表面与胶片之间的摩擦热;胶片在宽度方向上的温度变化很小, 故取胶片出鼓截面处的面平均温度为胶片出鼓温度。
在冷却过程中,胶片移动,冷却鼓转动,胶片与冷却鼓接触面不断变化,属于固体递进移动边界传热问题,目前难以进行有效的数值模拟。为此本文将固体胶片假设为流体胶片,以流体胶片的流动来代替固体胶片的递进移动,通过热-流-固耦合换热求解胶片与冷却鼓之间的传热。通过模拟值与实测值对比确定胶片的导热系数。
由于胶片厚度很小,移动速度较慢,因此认为流体胶片的流动为层流,其传热形式只有传导。在工程实际中,胶片导热系数随胶片材料组分不同而不同。因此,在模拟中先在胶片导热系数范围内选取流体胶片的导热系数,根据某一个工况的模拟结果与实测结果对比是否满足误差要求来确定导热系数值,并假设该值不受胶片厚度和温度的影响,也即可用于其他5种工况。事实上由于胶片很薄,导热路径很短,因此导热系数影响较小。
此外计算中还作了以下假设:1)假定流体的传热过程为稳定状态;2)流动状态为定常且不可压缩状态;3)由于密度差异引起的浮升力和重力忽略不计;4)不考虑流体流动时的黏性耗散所产生的热效应。
2.1 网格独立性验证由于整个胶片-冷却鼓模型大、网格数量多,为减少计算量,提高计算效率,采用1个冷却鼓及其鼓上胶片来进行网格独立性验证,结果如图 3所示,可以看出网格数量超过313万后,计算结果变化较小。因此,本文选择单鼓网格数量为313万的网格进行计算, 其有限元网格模型如图 4所示。
本文采用多参考系模型(multiple reference frame model, MRF)进行稳态计算。当胶片、冷却鼓和冷却水三者之间换热达到稳定状态时,在MRF模型中对三者进行稳态温度场分析。
2.3 湍流模型选择管道入口处流速为
$ {v_{\rm m}} = \frac{Q}{A} = \frac{{137.29 \times {{10}^{ - 3}}/60}}{{10 \times \frac{{\rm \pi} }{4} \times {{0.03}^2}}} = 0.325\ 356\ {\rm m/s} $ | (6) |
管道入口处雷诺数为
$ Re = \frac{{\rho vd}}{\mu } = \frac{{998.2 \times 0.325\ 356 \times 0.03}}{{0.001\ 003}} = 9\ 713.98 $ | (7) |
螺旋流道内雷诺数为
$ Re = \frac{{\rho \omega {H^2}}}{\mu } = \frac{{998.2 \times 0.416\ 67 \times {{0.0864}^2}}}{{0.001\ 003}} = 12\ 960 $ | (8) |
式中,ρ为冷却水的密度,kg/m3;vm为入口边界的平均速度,m/s;d为入口管道的特征长度,m;H为螺旋流道的特征长度,m;μ为冷却水的动力黏度,kg/(m·s);ω为冷却鼓转动角速度,rad/s。临界雷诺数随流体区域形状的不同而有所差别,一般情况下,根据雷诺数大小,可以将流动形态作如下划分:
1) Re<2 000为层流区,流动总是表现为层流;
2) 2 000≤Re≤4 000为过渡区,流动既可能是层流,也可能是湍流,与外界扰动情况有关;
3) Re>4 000为湍流区,流动一般呈现为湍流。
综上分析,冷却水流动模型选择湍流模型。通过查阅相关资料[10],选用Realizable k-ε模型,该模型适用于涉及旋转、漩涡及局部过渡流的复杂剪切流动。壁面函数选择Scalable Wall Function。
2.4 材料属性设置在Fluent材料数据库中定义相关材料,材料属性如表 3所示。如果使用k-ε模型或Spalart-Allmaras模型来模拟湍流,则可以在指定的流体区域关闭湍流模拟,“Laminar Zone”选项会将湍流黏度设置为零,并禁用流体区域中的湍流产生[11]。为了进一步验证并提高计算效率,本文采用单个冷却鼓计算胶片在不同动力黏度下冷却水和胶片的出鼓温度,结果如表 4所示,从表 4可以看出动力黏度的改变对此计算结果没有影响。
1) 冷却水入口采用速度入口边界条件,其入口速度由冷却水流量与入口管径计算得出。入口温度15.6 ℃。湍流特性使用湍流强度I与水力直径dhydro定义,其计算公式如下。
$ \begin{array}{*{20}{l}} {I = 0.16\cdot R{e^{ - \frac{1}{8}}}} \end{array} $ | (9) |
$ {d_{\rm hydro}} = \frac{{4{A_{润湿截面}}}}{{{L_{润湿截面}}}} = {D_{入}} $ | (10) |
2) 出口为压力出口边界。
3) 胶片与冷却鼓外表面、内表面和冷却水的交界面为interface交界面,通过耦合求解胶片、冷却鼓和冷却水之间的换热。
4) 胶片与空气接触的外表面及胶片侧面设置为对流换热边界条件。生产过程中, 胶片表面温度为54 ℃, 空气温度为14 ℃, 则取平均温度为29 ℃, 取此温度下的一系列参数来计算胶片与空气的对流换热系数。查阅资料得到此平均温度下空气的基本物性参数[12]:空气的密度ρair=1.128 kg/m3, 空气的动力黏度μair=1.91×10-5 kg/(m·s), 空气的导热系数λair=2.76×10-2 W/(m·K)。普朗特数Pr=0.699, 胶片速度v=0.15 m/s,胶片宽度L=0.76 m,故雷诺数
$ Re = \frac{{{\rho _{\rm air}}vL}}{{{\mu _{\rm air}}}} = \frac{{1.128 \times 0.15 \times 0.76}}{{1.91 \times {{10}^{ - 5}}}} = 6\ 732.56 $ | (11) |
胶片与空气的对流换热系数h1公式为[13]
$ \begin{array}{l} \ \ \ \ \ \ \ {h_1} = 0.402\frac{{{\lambda _{air}}}}{L}R{e^{1/2}}\frac{{P{r^{1/3}}}}{{{{[1 + {{(0.033\ 6/Pr)}^{2/3}}]}^{1/4}}}} + \\ 0.45 = 0.402 \times \frac{{2.76 \times {{10}^{ - 2}}}}{{0.76}} \times {6\ 732.56^{1/2}} \times \\ \frac{{{{0.699}^{1/3}}}}{{{{[1 + {{(0.033\ 6/0.699)}^{2/3}}]}^{1/4}}}} + 0.45 = \\ 1.48\ {\rm W/({m^2} \cdot K)} \end{array} $ | (12) |
5) 冷却鼓上下端面与空气接触的表面及内筒体内侧与空气接触的表面均设置为对流换热边界条件,由于胶片在冷却过程中主要热量被冷却水带走的同时,一部分热量被冷却鼓吸收,使鼓体温度升高。鼓筒表面温度为35 ℃, 空气温度为14 ℃, 则取平均温度25 ℃,查阅资料得到此平均温度下空气基本物性参数[12]:空气的运动黏度νair=1.65×10-5 kg/(m·s);空气的导热系数λair=2.71×10-2 W/(m·K);α为体积变化系数,对于理想气体即等于绝对温度的倒数;D为冷却鼓外径;Nu为努塞尔数,Gr为格拉晓夫数,普朗特数Pr=0.701, C取0.53,n取0.25,冷却鼓与空气的对流换热系数h2由以下公式确定[9]
$ \left\{ \begin{array}{l} Gr = \frac{{\alpha \cdot g \cdot \Delta T \cdot {D^3}}}{{{\upsilon ^2}}}\\ Nu = C \cdot {(Gr \cdot Pr)^n}\\ h = \frac{{{\lambda _{\rm air}}}}{D}Nu \end{array} \right. $ | (13) |
$ \begin{array}{l} \ \ \ \ \ \ \ {h_2} = \frac{{{\lambda _{\rm air}}}}{D} \cdot C \cdot {\left( {\frac{{\alpha \cdot g \cdot \Delta T \cdot {D^3}}}{{{\upsilon ^2}}} \cdot Pr} \right)^n} = \\ \frac{{2.71 \times {{10}^{ - 2}}}}{{0.72}} \times 0.53 \times \\ {\left( {\frac{{\frac{1}{{298}} \times 9.81 \times (308 - 287) \times {{0.72}^3}}}{{{{(1.65 \times {{10}^{ - 5}})}^2}}} \times 0.701} \right)^{1/4}} = \\ 3.2\ {\rm W/({m^2} \cdot K)} \end{array} $ | (14) |
选择SIMPLEC算法,此算法适用于大多数常规的不可压缩流动问题,并且在四边形或者六面体网格上更具有优势。求解过程采用稳态求解,其空间离散格式均采用二阶迎风格式进行计算。松弛因子及残差控制均保持默认设置。
3 结果与讨论 3.1 温度场分布胶片表面温度分布如图 5所示,从胶片入鼓至出鼓,温度分布逐渐降低;胶片通过鼓1的温度分布如图 6所示,胶片两端受对流换热的冷却作用,两端的温度略低于中间温度。另外,由于冷却水沿冷却鼓轴向流动,因此胶片靠近冷却水出口侧的温度略高于入口侧的温度。
冷却鼓外表面温度分布如图 7所示,鼓1至鼓10的温度分布梯度逐渐降低。鼓1外表面温度分布如图 8所示,沿冷却鼓轴向,在和胶片接触区域温度高,对于冷却鼓端部温度来说入口端略低于另外出口端。
取胶片出鼓截面处的面平均温度记为胶片出鼓温度,冷却水出口截面处的面平均温度记为冷却水的出鼓温度。
表 5~12为工况1~工况6的模拟结果以及模拟结果与实测数据的对比,可以得到如下分析。
1) 从鼓1至鼓10的冷却过程中,胶片出鼓温度、冷却水出鼓温度和冷却鼓平均温度都在逐渐下降,但下降幅度不断减小。这是由于胶片经过每一个鼓后,温度降低,胶片与冷却水之间温差逐渐减小,导致换热速率降低,进而使温度下降幅度减小。
2) 冷却水出鼓温度的模拟结果与实测数据对比,相对误差均低于6%,如图 9所示。胶片离开第十个冷却鼓的温度的模拟结果与实测数据对比,相对误差均低于8%,说明模拟方法是可行的。
3) 如图 10所示,从鼓1至鼓10,冷却水进出鼓温差变小,说明前鼓比后鼓的冷却效率高。
4) 胶片与冷却鼓在换热过程中通过各换热面的能量见表 6,胶片放出的能量主要由冷却水吸收,向空气中散发的能量可以忽略不计。由此可见冷却鼓在胶片冷却过程中起着至关重要的作用。
4 结论(1) 本文将固体胶片假设为流体胶片,以流体胶片的流动来代替固体胶片的递进移动,应用Fluent软件建立了胶片-冷却鼓三维耦合换热模型,对胶片在冷却鼓上的冷却过程进行热-流-固耦合仿真分析, 得到胶片及各冷却鼓的温度分布规律,结果表明胶片区域整体换热均匀,不会出现局部与整体的换热速率差距较大的情况。
(2) 对6组不同的实际胶片冷却过程进行了仿真模拟,并将各冷却鼓出水温度以及胶片的出鼓温度分别与实测温度进行对比,得到冷却水出鼓温度相对误差均低于6%,胶片出鼓温度的相对误差均低于8%,说明针对与胶片冷却相类似的耦合强度较弱且相对移动稳定的工程问题,本文提出的流动传热模拟方法可以用于工程计算,并有望解决类似的固体递进移动边界传热模拟难题。
[1] |
杨云, 董春光. 水冷却辊的加工[J]. 江苏冶金, 2004, 32(3): 40-41. YANG Y, DONG C G. Water cooling roller processing[J]. Jiangsu Metallurgy, 2004, 32(3): 40-41. (in Chinese) DOI:10.3969/j.issn.1005-6068.2004.03.014 |
[2] |
陈洪冰, 眭向荣, 孟晓霞, 等. 层流冷却辊喷焊制造技术[J]. 江苏冶金, 2008, 36(3): 50-53. CHEN H B, SUI X R, MENG X X, et al. Spraywelding manufacturing technology of laminar cooling roller[J]. Jiangsu Metallurgy, 2008, 36(3): 50-53. (in Chinese) DOI:10.3969/j.issn.1005-6068.2008.03.017 |
[3] |
杨宏伟, 孙凯, 田永林. 全钢压延生产线冷却辊及张力控制装置的改进[J]. 橡塑技术与装备, 2009, 35(3): 35-38. YANG H W, SUN K, TIAN Y L. Improvement of cooling roller and tension control device in whole steel rolling production line[J]. China Rubber/Plastics Technology and Equipment, 2009, 35(3): 35-38. (in Chinese) DOI:10.3969/j.issn.1009-797X.2009.03.007 |
[4] |
韩亚娜, 范勋铭, 曲荣旭, 等. 大型绿色低温一次法混炼生产线冷却辊筒机械加工工艺设计[J]. 橡塑技术与装备(橡胶版), 2015, 41(1): 54-60. HAN Y N, FAN X M, QU R X, et al. Process design of cooling roller machining for large green single stage mixed production line at low temperature[J]. China Rubber/Plastics Technology and Equipment(Rubber edition), 2015, 41(1): 54-60. (in Chinese) |
[5] |
郭茜, 严密. 平板流铸工艺中冷却辊的温度场数值模拟[J]. 稀有金属材料与工程, 2015, 44(8): 2048-2052. GUO X, YAN M. Thermal analysis for cooling rolls in planar-flow melt spinning[J]. Rare Metal Materials and Engineering, 2015, 44(8): 2048-2052. (in Chinese) |
[6] |
陈北荣. 凹印机冷却辊的故障分析与局部改造设计[J]. 常州工学院学报, 2016, 29(1): 31-33. CHEN B R. Fault analysis and partial innovation design of the cooling roll of gravure press[J]. Journal of Changzhou Institute of Technology, 2016, 29(1): 31-33. (in Chinese) DOI:10.3969/j.issn.1671-0436.2016.01.007 |
[7] |
张羽玲, 陈允春, 邵明月, 等. 基于MRF模型的冷却辊冷却性能的CFD分析[J]. 包装工程, 2017, 38(7): 11-14. ZHANG Y L, CHEN Y C, SHAO M Y, et al. CFD analysis of cooling performance of cooling roller based on MRF model[J]. Packaging Engineering, 2017, 38(7): 11-14. (in Chinese) |
[8] |
唐家鹏. ANSYS FLUENT 16.0超级学习手册[M]. 北京: 人民邮电出版社, 2016. TANG J P. ANSYS FLUENT 16.0 super learning manual[M]. Beijing: Posts & Telecom Press, 2016. (in Chinese) |
[9] |
傅秦生. 热工基础与应用[M]. 3版. 西安: 西安交通大学出版社, 2016. FU Q S. Fundamentals and applications of thermal engineering[M]. 3rd ed. Xi'an: Xi'an Jiaotong University Press, 2016. (in Chinese) |
[10] |
胡坤, 顾中浩, 马海峰. ANSYS CFD疑难问题实例详解[M]. 北京: 人民邮电出版社, 2017. HU K, GU Z H, MA H F. Detailed examples of ANSYS CFD problems[M]. Beijing: Posts & Telecom Press, 2017. (in Chinese) |
[11] |
段中喆. ANSYS FLUENT流体分析与工程实例[M]. 北京: 电子工业出版社, 2015. DUAN Z Z. ANSYS FLUENT fluid analysis and engineering examples[M]. Beijing: Publishing House of Electronics Industry, 2015. (in Chinese) |
[12] |
任泽霈, 蔡睿贤. 热工手册[M]. 北京: 机械工业出版社, 2002. REN Z P, CAI R X. Thermal engineering manual[M]. Beijing: China Machine Industry Press, 2002. (in Chinese) |
[13] |
CHURCHILL S W, OZOE H. Correlation for laminar forced convection with uniform heating in flow over a plate and in developing and fully developed flow in a tube[J]. Journal of Heat Transfer, 1973, 95(1): 78-84. DOI:10.1115/1.3450009 |