基于滑移边界的饼干面团双辊挤压成型的数值模拟.pdf

己引田格式甲文马献昌龙继兰竿其干语移边界的饼而团压成型的数估桔圳O0O2A5/Q1.5A1.550力第45 卷第3 期2023年6 月实践学基于滑移边界的饼干面团双辊挤压成型的数值模拟1)田文马戴晨龙杜继芸袁方洋2)（江苏省食品先进制造装备技术重点实验室，江苏无锡2 1412 2)（江南大学机械工程学院，江苏无锡2 1412 2)摘要饼干是广受欢迎的零食之一，饼干面团的量化生产要求各个面团的挤出重量保持稳定性和一致性。饼干面团属于非牛顿流体，采用Polyflow软件对饼干面团在双辊挤压成型过程进行数值模拟，分析滑移系数和辊子转速之间的关系以及饼干面团挤出流量的均匀性。结果表明：面团出口流量具有不均匀性，五个孔道出口流量呈现中间大两边小的分布。滑移系数和辊子转速呈线性正相关，数值模拟可预测出口平均质量流量，且预测精度不会随转速的增加而降低。关键词非牛顿流体，挤出成型，滑移系数，流量不均匀性中图分类号：O373，T S2 0 3文献标识码：Adoi:10.6052/1000-0879-22-421NUMERICAL SIMULATION OF DOUBLE ROLL EXTRUSION OF BISCUITDOUGH BASED ON SLIP BOUNDARY CONDITIONS!)TIANWenmaDAI Chenlong DU JiyunYUANFangyang2)(Jiangsu Key Laboratory of Advanced Food Manufacturing Equipment&Technology,Wuxi 214122,Jiangsu,China)(School of Mechanical Engineering,Jiangnan University,Wuxi 214122,Jiangsu,China)Abstractt Biscuits are one of the most popular snacks,and the mass production of biscuit dough requires theextrusion weight of each dough to remain stable and consistent.The biscuit dough is a non-Newtonian fluid.The Polyflow software was used to simulate the biscuit dough extrusion process in a double-roller extruder.Therelationship between the slip coefficient and the roller speed,and the uniformity of the outlet flow rate of thebiscuit dough were analyzed.The results showed that the dough outlet flow rate was ununiform.For the fiveoutlet orifices configuration,the extrusion mass flow rate on both sides was lighter than that inside.The slipcoefficient increased linearly with the roller rotation speed.Numerical simulation can predict the averageextrusion mass flow rate,and the prediction accuracy does not decrease with the increase of the rotationalspeed.Keywordsnon-Newtonian fluid,extrusion molding,slip coefficient,flow nonuniformity饼干是常见的烘焙食品之一，因其货架期长、口味多样、质地多样、消费范围广等吸引人的特点而占据着市场显著地位。如今，饼干生产已经从劳动密集型的手工业发展到高度机械化的程度。饼干面团的成型是饼干生产的重要环节，其方法包括辊压、冲印、挤条成型等。软面团或含有巧克力片或坚果的面团越来越受到消费者的欢迎，它们必须通过挤压（结合线切割）而不是模压成型，常用于制作高品质的饼干产品。饼干面团的挤压成型在机器展向存在重量分布不均的问题，2022-07-21收到第1稿，2 0 2 2-0 9-0 6 收到修改稿。1）国家自然科学基金（12 17 2 152）和江苏省食品先进制造装备技术重点实验室（FMZ202206）资助项目。2）袁方洋，副教授，研究方向为多相流体力学与食品机械。E-mail：f y y u a n j i a n g n a n.e d u.c n1面团双批折压成至的效值模拟：万字一买域，2 0 2 5，45（5）：544-50Tian Wenma,Dai Chenlong,Du Jiyun,et al.Numerical simulation of double roll extrusion of biscuit dough based on slipboundary conditions.Mechanics in Engineering,2023,45(3):544-550545田文马等：基于滑移边界的饼干面团双辊挤压成型的数值模拟第3 期虽然通过引入隔板来限制过流面积可对挤出流量进行调节，但隔板的进给量、辊子转速的设置依然需要经验或是大量调试确定。面团是具备非牛顿特性的复杂流体，对面包面团挤出所开展的研究较多。Sofou等 2 对室温下的面包面团进行了流变性表征，发现当时间效应不显著时，可以使用Herschel-Bulkley黏塑性模型描述该面团。面包面团水平辊压过程的实验与有限元模拟结果表明，面团挤出过程表现出明显的壁面滑移，需确定适当的滑移规律 3-4。不同于面包面团，饼干的制备原料采用短筋面团，具有较强黏性但缺乏延展性和弹性。在屈服范围内，面团的应变呈现出黏弹性；在大变形和超过屈服值时，由于缺乏连续的蛋白质相网络，面团表现为黏塑性 5。Peck等 开发了饼干面团水平辊压的模型并与实验对照，发现基于幂律模型及Levine润滑假设的面团流变学建模给出了偏低的力学预测。如果将数据拟合到Levine模型，表观幂律流变学建模可预测饼干面团水平压片过程。Chak-rabarti-Bell 等 7 使用Abaqus模拟了面团压片过程，预测的面团厚度与实验数据非常吻合；对无酵面团的挤压建模表明，非线性黏弹性模型与Mullins模型结合可以很好地捕捉连续体的流变特性8。Padilla等 9 采用ANSYS数值模拟了所设计的四孔道活塞挤出机，以避免较大的流动死区。周克 10 基于Polyflow仿真平台对单螺杆挤压工艺进行数值研究，发现转速对各参数正相关。饼干面团双辊挤压成型过程中，壁面的滑移效应不可忽视。文献 2-411研究发现，滑移的存在会使挤出压力、辊分离力和输入功率等操作参数减小。张迅等 13 在研究同向啮合三螺杆挤出机混合性能时考虑滑移效应的影响，发现滑移边界条件对数值模拟的结果影响很大。Malik等 14研究了压力对双螺杆挤出机壁面滑移的影响，发现机筒表面的滑移相较螺杆表面具有更显著的影响，滑移系数是压力相关的函数。Lewandowski 等 15基于Polyflow对单螺杆挤出过程中具有滑移效应的聚合物熔体流动进行了有限元模拟，发现在正压力梯度范围内，滑移会降低螺杆的泵送能力。本文对饼干面团的挤压成型过程进行数值研究，研究展向存在的重量分布不均规律。基于Polyflow平台，构建非牛顿流动模型，着重考虑双辊挤压时的壁面滑移特性。通过对预测的挤出流量与实验数据相比较，验证数值模型的有效性。1饼干面团双辊挤压成型的问题描述图1为双辊面团挤出机的结构简图，从上至下依次为进料斗、一对旋转方向相反的柱状圆辊、压力室、模具。挤出机通常位于烤箱带上方，依靠双辊的持续转动，面团被不断地挤压出模具，被张紧的钢丝或刀片切断，形成的面团薄片落在皮带或传送带上，送至烤箱烤制成型。feedhopperrollersGfillerblockmold图1双辊面团挤出机Fig.1Dough extruder with double rollers双辊面团挤出机的优点是在流程上能连续生产，即面团通过辊子的转速来实现面团重量调节，且能在展向增加挤出压力室出口的数量来提高生产效率。小型挤出机的压力室有数个出口，大型产线则可达数十个甚至更多。压力室的设置有助于面团在腔室流道的均匀流动，但双辊挤压成型力546实2023年第45卷践学方法存在着机器展向面团重量分布不均的问题，即机器宽度范围内挤压速率的差异一两侧的重量往往比中间轻。面团块的重量差异会显著影响面团后续的烘焙制作。为解决该问题，通常采用在压力室内引入隔板，通过限制面团在模具中的过流截面面积来调整各孔道的面团挤出流量。隔板的进给位置通常依靠经验在生产时实时调节，既浪费时间又浪费物料。然而，不同批次面团的性质、面团的堆积高度等都会影响面团的挤压速率；并且，在短时间内获得各个孔道挤出的面团重量也相当困难。因此，获得精确的面团成型孔道分布模型对饼干的标准化生产有重要意义。2数值模型2.1面团流体的控制方程将面团视为连续介质，考虑面团为非牛顿流体，则面团经过辊子挤压流动的质量和动量守恒方程为p+V.(pu)=0(1)tdu=f-Vp+v2u(2)dtPP式中p为密度，t是时间，u是速度矢量，fb是单位质量的体积力，p是压力，是动力黏度。2.2黏度模型饼干面团是固液糊状物，具有复杂的流变性。低筋饼干面团具有假塑性流体特征即剪切变稀的流体特性。由于挤压过程中线速度不大，剪切速率导致的黏度改变并不剧烈，选择应用广泛的Bird-Carreau黏度模型ln-1=o1+()22(3)式中为黏度，o为零剪切黏度，入为松弛时间，为剪切率，n为非牛顿指数。其中o，入和n分别取值10 0 0 0 Pa.s16,1.5 s17,0.718。2.3滑移边界条件面团挤出作业一段时间后，辊子沟槽填充了大量面团物料，如图2 所示。此时，双辊挤压成型的表面是由辊子齿间与沟槽面团构成的复合接触面（图2(b)），导致两辊间面团产生明显的滑移。滑移对辊压挤出有重要影响，单纯设置无滑移条件无法描述实际的流动现象。此处采用广义纳维线性定律描述滑移效应 13 1,19 Ts=(uw-ut)(4)式中T。为壁面切向滑移力，为滑移系数，uw为壁面的切向速度，t为流体的切向速度。当=0时，完全滑移；当=时，无滑移。滑移系数越大，面团与壁面间的摩擦力越大。滑移长度与物料性质、辊子结构（特别是沟槽的尺寸）、壁面粗糙度相关，目前尚无法直接确定饼干面团挤出过程的辊子表面滑移系数值，因此对滑移系数的表征是获得该过程滑移特性的核心。（a）工作过程中的辊子表面(b）填充有面团物料的辊子截面(a)Rollersurfaceduringwork(b)Roller section filled with dough图2 挤压成型过程的辊子表面特征Fig.2Roller surface characteristics in extrusion process2.4计算设置计算模型为小规模生产所用的有5个压力腔出口的小型挤出机。采用四面体网格将模型离散化，网格最大尺寸为5mm，如图3 所示。面团547田文马等：基于滑移边界的饼干面团双辊挤压成型的数值模拟第3 期密度为1150 kg/m，认为初始时刻进料斗充满面团原料。双辊以相反方向向下挤压面团。0250500 mm图3 五孔面团挤出模型网格划分Fig.3Mesh generation of the dough extrusion modelwithfiveoutlets本文基于Polyflow仿真平台进行数值模拟。Polyflow具有网格重置技术，适合转动流场的模拟。它专门用于黏弹性材料模拟，适用于塑料、树脂等高分子材料的挤出成型、吹塑成型、拉丝、层流混合、涂层过程中的流动及传热和化学反应问题。3结果与讨论3.1面团挤出过程分析选取辊子转速n=0.625r/min的工况，对挤出通道的剪切力、压力和速度云图进行分析。如图4(a)所示，沿辊子周围的面团受到剪切力较为明显；如图4(b)所示，压力腔的边缘位置对面团的剪切力较大。图5所示为面团挤出过程中所受压力分布，根据云图显示，面团在压力腔的入口位置及通道的中间位置所受压力较大，到通道后半段，由于出口处面团的不断挤出，压力变小。图6 所示为面团挤出过程中的速度分布，在辊子周边，由于辊子的转动，面团的速度较大，同时由于辊子带动，在上方料斗中的面团也以较小的速度运动。面团的最大速度出现在两个辊子中间位置，即压力腔的入口位置，此处由于辊子挤压效果最明显，面团速度较大。local shearratelocal shearrate4.540e+004.540e+004.301e+004.301e+004.062e+004.062e+003.823e+003.823e+003.584e+003.584e+003.345e+003.345e+003.106e+003.106e+002.868e+002.868e+002.629e+002.629e+002.390e+092.390e+002.151e+002.151e+001.912e+001.912e+001.673e+001.673e+001.434e+081.434e+001.195e+001.195e+009.558e-019.558e-017.169e-017.169e-014.779e-014.779e-012.390e-012.390e-0100S-1(a)纵向剖面(b)横向剖面(a)Longitudinal section(b)Cross section图4面团挤出通道的剪切率分布Fig.4Shear rate distribution of the dough extrusion channelpressurepressure9.121e+049.121e+048.313e+048.313e+047.505e+047.505e+046.697e+046.697e+045.888e+045.888e+045.080e+045.080e+044.272e+044.272e+043.646e+043.646e+042.656e+042.656e+041.848e+041.848e+041.040e+041.040e+042.320e+042.320e+04-5.760e+03-5.760e+031.384e+04-1.384e+04-2.192e+04-2.192e+04-3.000e+04-3.000e+04-3.808e+04-3.808e+04-4.616e+04-4.616e+04-5.424e+04-5.424e+04-6.232e+04-6.232e+04PaPa(a)纵向剖面(b)横向剖面(a)Longitudinal section(b)Cross section图5面面团挤出通道的压力分布Fig.5Pressure distribution of the dough extrusion channel力548实2023年第45卷践学velocitiesvelocities1.092e-021.092e-021.035e-021.035e-029.773e-039.773e-039.198e-039.198e-038.623e-038.623e-038.049e-038.049e-037.474e-037.474e-036.899e-036.899e-036.324e-036.324e-035.749e-035.749e-035.174e-035.174e-034.599e-034.599e-034.024e-034.024e-033.449e-033.449e-032.874e-032.874e-032.300e-032.300e-031.725e-031.725e-031.150e-041.150e-045.749e+005.749e+0000ms-1 ms-(a)纵向剖面(b)横向剖面(a)Longitudinal section(b)Cross section图6面团挤出通道的速度分布Fig.6Velocity distribution of the dough extrusion channel3.2滑移系数的确定滑移系数的取值目前没有具体的参考和经验值。选取辊子转速分别为0.6 2 5r/min，1.5r/m i n，2.4r/min,3.64r/min,4.55r/min,6.135 r/min,并且Polyflow软件模拟的滑移系数值相应地分别设置为50 0，12 0 0，16 50，2 3 0 0，3 2 0 0，43 0 0。图7 给出了辊子转速与出口平均质量流量的关系。出口平均质量流量为挤出机稳定工作时五孔道的出口质量流量平均值。辊子转速相同情况下，实验和模拟的出口平均速度基本吻合（平均偏差为0.61%），且随着辊子转速的增加，出口平均速度也相应地线性增大。图8 给出了滑移系数随辊子转速n的变化关系。辊子转速越大，滑移系数也越大。采用线性方程拟合两者之间的函数关系得=701.5n(5)3.5exp.simu.3.02.52.01.51.00.5001234567roller speed/(r-min-1)图7辊子转速和出口平均质量流量的关系Fig.7 Relationship between roller speed and average massflowoftheoutlets可以看到，滑移系数与辊子转速成正比，即与辊子表面线速度成正比。又有研究表明滑移系数与表面压力成正比14，因此壁面表面压力与速度成正相关，与图5和图6 的结果一致。3.3展向流量不均匀性分析采用式（5）给出的模型，模拟获得了不同辊子转速时的五孔道挤出质量流量分布，并与实验数据对比。如图9 所示，为n=1.5r/min和2.4r/min时的五孔道出口质量流量分布曲线，以及模拟结果基于实验数据的相对误差。总体而言，实验与模拟结果相差不大，并且两者都显示中间孔道的挤出流量更大，而靠近侧边孔道的流量较小，这与Manleyl研究结果一致。模拟所得的流量不均匀性更明显，即中间孔道的流量较实验更大，而两边的孔道流量较实验更小，这可能是由于压力室壁面的滑移特性与辊子表面存在差异 14。X1035simu.-=701.5n4321001234567roller speed/(r:min-1)图：车辊子转速和滑移系数的关系Fig.8Relationship between roller speed and slip coefficient549田文马等：基于滑移边界的饼干面团双辊挤压成型的数值模拟第3 期(r-ur.)/eimo51.4n=1.5 r/min4n=2.4 r/min1.231.020.8sseu0.60.4simulation n=2.4 r/min0experiment n=2.4 r/min0.2o-simulation n=1.5 r/min1-experiment n=1.5 r/min021234512345serial number of orificesserial numberof orifices(a）平均出口质量流量（b）预测流量的相对误差(a)Average outlet mass flowrate(b)Relative error of predicted flowrate图9 五孔道出口质量流量分布Fig.9Mass flow distribution of five orifice outlets基于实验结果比较不同转速下的数值预测结果，发现相对误差均可以控制在土5%范围内，满足工程应用要求。并且，五孔道的相对误差呈现相同的分布，相对误差大小似乎不受转速改变的影响，即模型的预测精度不会随转速的增加而降低。但是，实验与模拟均未能给出较好的展向分布对称性，可能来自油脂的影响 2 0 及面团密度的非均一、老化等因素 6。式（5）给出的模型较为简单，精确的预估模型需进一步考虑这些因素的影响。4结论本文主要探究了在双辊挤压条件下，饼干面团通过五孔道时出口流量的分析，采用Polyflow仿真软件，基于滑移边界条件和选取的黏度模型进行了数值模拟计算，结论如下。（1）挤压过程中，辊子周边的面团受到剪切力明显，面团在压力室的附近受剪切力较大。面团在压力室的入口位置及通道的中间位置受压力较大，到通道后半段，由于出口处面团的不断挤出，压力变小。面团的最大速度出现在两个辊子中间位置，即压力室的入口位置。（2）在相同的辊子转速下，通过对比实验和模拟的出口流量，确定了辊子转速和滑移系数的线性关系。该模型对五孔道面团挤出流量的预测误差在土5%范围内，满足工程应用要求，且模型的预测精度不会随转速的增加而降低。（3）在两个不同辊子转速下，孔道面团展向流量表现出不均匀性，可能是压力室表面滑移、油脂流失或者面团密度非均一、老化的影响。本研究不仅对饼干的标准化生产提供技术指导，还可对3 D食品打印等技术提供理论依据。参考文献1 Manley D.Manleys Technology of Biscuits,Crackers andCookies,Fourth edition.UK:Woodhead Publishing Lim-ited,20112 Sofou S,Muliawan EB,Hatzikiriakos SG,et al.Rheologic-al characterization and constitutive modeling of breaddough.Rheologica Acta,2008,47(4):369-3813 Mitsoulis E.Numerical simulation of calendering visco-plastic fluids.Journal of Non-newtonian Fluid Mechanics,2008,154(2-3):77-884 Mitsoulis E,Hatzikiriakos SG.Rolling of bread dough:ex-periments and simulations.Food and Bioproducts Pro-cessing,2009,87(2):124-1385 Arepally D,Reddy RS,Goswami TK,et al.Biscuit baking:a review.LWT-Food Science and Technology,2020:1097266 Peck MC,Rough SL,Barnes J,et al.Roller extrusion ofbiscuit doughs.Journal of Food Engineering,2006,74(4):431-4507 Chakrabarti-Bell S,Bergstrom JS,Lindskog E,et al.Com-putational modeling of dough sheeting and physical inter-pretation of the non-linear rheological behavior of wheatflour dough.Journal of Food Engineering,2010,100(2):278-2888 Mohammed MAP,Wanigasooriya L,Chakrabarti-Bell S,etal.Extrusion of unleavened bread dough:experiments andsimulations.Journal of Rheology,2017,61(1):49-659 Padilla C,Vivanco A,Vinces L,et al.Design of a multi-hole cylindrical extruder,driven by a linear actuator and（责任编辑：王永会）力实5502023年第45卷践学used for the formation of bakery dough.2020 IEEE XXVIIInternational Conference on Electronics,Lima,Peru,202010周克.单螺杆螺压过程推进剂流变参数及物料混合特性的数值模拟研究.博士论文.南京：南京理工大学，2 0 15Zhou Ke.Numerical simulation on the rheological paramet-er variation and mixing ability of single-screw extrusionprocess of propellant.PhD Thesis.Nanjing:Nanjing Uni-versity of Science&Technology,2015(in Chinese)11 Mitsoulis E,Sofou S.Calendering pseudoplastic and visco-plastic fluids with slip at the roll surface.Journal of AppliedMechanics,2006,73(2):291-29912 Karapetsas G,Mitsoulis E.Some experiences with the slipboundary condition in viscous and viscoelastic flows.Journ-al of Non-Newtonian Fluid Mechanics,2013,198:96-10813张迅，朱向哲，王阁.滑移边界条件对同向啮合三螺杆挤出机混合性能的影响.石油化工高等学校学报，2 0 12,2 5(4)：6 8-7 2Zhang Xun,Zhu Xiangzhe,Wang Ge.Numeral simulationof the effects of wall slip on an intermeshing co-rotating tri-screw extruder in the cross-section.Journal of Petrochemic-al Universities,2012,25(4):68-72(in Chinese)14 Malik M,Kalyon DM,Golba JC.Simulation of co-rotatingtwin screw extrusion process subject to pressure-dependentwall slip at barrel and screw surfaces:3D FEM analysis forcombinations of forward-and reverse-conveying screw ele-ments.International Polymer Processing,2014,29(1):51-6215 Lewandowski A,Wilczynski K.Global modeling of singlescrew extrusion with slip effects.International Polymer Pro-cessing,2019,34(1):81-9016 Thacker D.The characteristics and processing require-ments of wheat for biscuit making.HGCA Project Report(No.83),UK,199317 Maache-Rezzoug Z,Bouvier JM,Allaf K,et al.Effect ofprincipal ingredients on rheological behaviour of biscuitdough and on quality of biscuits.Journal of Food Engineer-ing,1998,35(1):23-4218 Baltsavias A.Mechanical properties of short doughs andtheir corresponding biscuits.Wageningen University andResearch,199619 Denn MM.Extrusion instabilities and wall slip.Annual Re-view of Fluid Mechanics,2001,33(1):265-28720 Devi A,Khatkar BS.Physicochemical,rheological andfunctional properties of fats and oils in relation to cookiequality:a review.Journal of Food Science and Technology,2016,53(10):3633-3641