不同螺槽深度下异向双螺杆挤出机流道仿真分析.pdf

基金项目：国家自然科学基金面上项目（编号：）作者简介：张一明，男，北京工商大学在读硕士研究生。通信作者：黄志刚（），男，北京工商大学教授，博士。：收稿日期：改回日期：犇 犗 犐：犼 狊 狆 犼 狓 文章编号 （）不同螺槽深度下异向双螺杆挤出机流道仿真分析 张一明，犣犎犃犖犌犢 犻 犿 犻 狀 犵，黄志刚，犎犝犃犖犌犣 犺 犻 犵 犪 狀 犵，杨亚楠，犢 犃犖犌犢 犪 狀 犪 狀，（北京工商大学人工智能学院，北京 ；塑料卫生与安全质量评价技术北京市重点实验室，北京 ）（犛 犮 犺 狅 狅 犾 狅 犳犃 狉 狋 犻 犳 犻 犮 犻 犪 犾 犐 狀 狋 犲 犾 犾 犻 犵 犲 狀 犮 犲，犅 犲 犻 狀 犵犜 犲 犮 犺 狀 狅 犾 狅 犵 狔犪 狀 犱犅 狌 狊 犻 狀 犲 狊 狊犝 狀 犻 狏 犲 狉 狊 犻 狋 狔，犅 犲 犻 犼 犻 狀 犵 ，犆 犺 犻 狀 犪；犘 犾 犪 狊 狋 犻 犮犅 犲 犻 犼 犻 狀 犵犕 狌 狀 犻 犮 犻 狆 犪 犾犓 犲 狔犔 犪 犫 狅 狉 犪 狋 狅 狉 狔狅 犳犎 犲 犪 犾 狋 犺犪 狀 犱犛 犪 犳 犲 狋 狔犙 狌 犪 犾 犻 狋 狔犈 狏 犪 犾 狌 犪 狋 犻 狅 狀犜 犲 犮 犺 狀 狅 犾 狅 犵 狔，犅 犲 犻 犼 犻 狀 犵 ，犆 犺 犻 狀 犪）摘要：目的：讨论聚乳酸在异向双螺杆挤出机不同螺槽深度下的流道分布规律。方法：使用 软件建立种不同螺槽深度的螺杆元件及其对应的流道模型，运用 流体仿真软件进行数值模拟，借助 软件对后处理结果进行数值分析；并对种不同螺槽深度状态下的异向双螺杆挤出机的聚乳酸流体的压力场、剪切速率场和黏度场，以及停留时间分布曲线 进行对比分析。结果：通过改变螺杆元件的齿根圆直径及其对应的流道直径可实现螺槽深度的改变，在其他条件相同的情况下，螺槽深度越浅，压力波动越大，剪切速率波动越大，但黏度波动越小；螺槽深度越深，轴向混合性能越高。结论：螺槽深度的变化对流场混合性能存在影响并具有一定的规律性，可根据试验选择混合性能全面的双螺杆挤出机进行实际生产。关键词：螺杆挤出机；异向啮合；螺槽深度；流场混合性；数值模拟；流场分析犃 犫 狊 狋 狉 犪 犮 狋：犗 犫 犼 犲 犮 狋 犻 狏 犲：（）犕 犲 狋 犺 狅 犱 狊：，（）（）犚 犲 狊 狌 犾 狋 狊：，犆 狅 狀 犮 犾 狌 狊 犻 狅 狀：犓 犲 狔 狑 狅 狉 犱 狊：；聚乳酸（）是一种环境友好型材料，同时也是一种具有完全可降解性、高生物相容性和机械性能优异的热塑性材料，可作为石油基塑料的替代品，已被广泛应用于包装、农业、纤维纺织、医用和工程塑料等领域。双螺杆挤出机由于具有良好的加料性能、混炼塑化性能、排气性能、挤出稳定性能等特点，已被广泛应用于挤出制品的成型加工。双螺杆挤出机分为同向双螺杆挤出机和异向双螺杆挤出机两大类。异向双螺杆由于特殊的构型，形成了型室，物料被封闭在彼此隔开的型室中，封闭在型室的物料不随螺杆转动，只沿螺杆轴线方向做正位移动，故其被广泛应用于挤出成型和配料造粒及其他方面。范涛等研究发现，减小螺槽深度会使熔体热量的传递效果提高、有利于熔体的平稳传递。张雪晨等发现，螺杆元件构型的变化对流道的剪切作用和建压能力 第 卷第期 总第 期 年月有优化作用，且螺槽越深，轴向混合能力越强。等研究发现，增加捏合块错列角会提高材料的温度，有利于提高转化率，同时还可提高挤出料的拉伸效率，延长材料在挤出 机中 的停 留 时 间，有 利 于 强 化 混 合 及 反 应。彭涛发现，开槽螺纹元件的建压输送能力较低，分散混合性能较弱，但其分布混合性能优于常规螺纹元件。综上，机械结构的设计对物料加工起到了决定性影响。例如单螺杆和双螺杆在进行轴向熔融输送时，单螺杆依靠摩擦力，而双螺杆依靠啮合螺纹的推力，这就导致了物料的内部分子键破坏方式不同，得到的产物性质发生改变。目前，有关螺槽深度对螺杆挤出机混合效果的影响研究集中在单螺杆和同向双螺杆方面，有关异向双螺杆螺纹元件螺槽深度对聚乳酸流道内流场影响的研究尚未见报道。研究拟针对聚乳酸在啮合异向双螺杆挤出机中的流道进行数值模拟，应用 软件建立不同螺槽深度的螺杆元件三维模型，利用 流体仿真软件进行数值模拟，并采用粒子示踪法模拟粒子在流道内的运动轨迹，研究不同螺槽深度对聚乳酸流场中的压力场、剪切速率场和黏度场，以及轴向停留时间分布的影响，选择出合适的异向双螺杆构型来提高聚乳酸的性能，降低损耗成本。理论模型的建立建立几何模型端面曲线的建立根据双螺杆几何学的基本原理，通过相对运动法在 上推导出全啮合异向双螺杆的理论端面曲线，得到理论端面曲线如图所示。推导出的全啮合双头异向双螺杆的理论端面曲线的段圆弧方程分别为：、段：狓犚犫 狑 狋狔犚犫 狑 狋狑 狋 狑 狋（），（）、段：狓犚狊 狑 狋狔犚狊 狑 狋狑 狋 狑 狋（），（）段：狓犆犔 狑 狋犚狊 狑 狋（）狔犆犔 狑 狋犚狊 狑 狋（）烅烄烆 犆犔犚狊狑 狋（），（）段：狓犆犔 狑 狋犚狊 狑 狋（）狔犆犔 狑 狋犚狊 狑 狋（）烅烄烆 狑 狋 犆犔犚狊（），（）段：狓犆犔 狑 狋犚狊 狑 狋（）狔犆犔 狑 狋犚狊 狑 狋（）烅烄烆 犆犔犚狊狑 狋（），（）段：狓犆犔 狑 狋犚狊 狑 狋（）狔犆犔 狑 狋犚狊 狑 狋（）烅烄烆 狑 狋 犆犔犚狊（）。（）图理论端面曲线 端面曲线的修正通过 软件螺旋扫描端面曲线即可得到全啮合异向双螺杆三维实体。实际应用中，如果依照理论曲线加工，螺杆之间完全没有间隙，若在螺杆加工中实现无间隙不仅对加工刀具和工艺要求非常高，而且由于热胀冷缩两螺杆很容易干涉，存在很大的安全隐患，同时在螺杆根部容易形成积料和死角。因此，实际应用中通常对理论曲线进行修正（图），使两螺杆螺棱侧面之间有一定侧间隙，修正后螺槽角 ，螺棱角 。图修正端面曲线 食品装备与智能制造 总第 期 年月几何模型的建立通过表、表在 软件中绘制如图所示的异向双螺杆三维模型，将左螺杆的中心点设定为坐标原点，犡轴方向向右，犢轴方向向上，轴的方向和旋转方向根据右手定则确定：左螺杆逆时针旋转、右螺杆顺时针旋转、犣轴方向与挤出方向相同。由于螺槽深度变化不大，因此 模型代表 和 组合得到的模型，以此类推，模型代表 和 组合得到的模型；模型代表 和 组合得到的模型。表螺杆元件建模参数 组合长度螺槽深度螺杆外径螺杆内径 表流道区域建模参数 组合长度中心距流道外径流道内径 图异向双螺杆三维模型 数学模型的建立基本假设在进行异向双螺杆仿真前，考虑到实际情况和环境等众多影响因素，作出以下假设，并在此基础上进行仿真分析：（）熔融体是不可压缩的且充满整个流道。（）熔融体在流场中的各个位置温度都是相同的。（）机筒壁面是无滑移的。（）重力与惯性力忽略不计，且远远小于黏性力的体积力。（）熔融体为雷诺系数较小的层流流动。基本方程熔体流动必须满足流体力学的三大基本方程：质量守恒方程（也称连续性方程）、动量守恒方程（也称运动方程）和能量守恒方程。因为在等温条件下进行，所以能量守恒方程不必求解。基于上述条件，基本方程为：连续性方程：?狏，（）运动方程：犘?犜，（）式中：?哈密尔顿算子；狏 速度矢量，；犘 流体静压力，；犜 应力张量，。其中，选用 本构方程 进行模拟数值计算，探究聚乳酸剪切速率和黏度之间的关系：（）（）（狀），（）式中：零剪切黏度，；剪切速率，；松弛时间，；无穷剪切黏度，；狀 流动指数。模拟过程设定为等温条件，温度设为 ；聚乳酸材料在 下的零剪切黏度为 、松弛时间为 、流体指数为 、无穷剪切黏度为 。有限元模型网格划分使用 软件将建立好的组螺杆元件和对应的流体区域进行网格划分，由于螺杆元件几何形状相对复杂，因此使用四面体网格直接进行划分；流体区域使用六面体规则形网格进行划分。划分后的螺杆元件和流体区域如图所示。图螺杆元件及流道的有限元模型 边界条件的设定流体区域入口流体和出口流体为自由流动，因此法向力和切向力设为；左右内孔属于滑移 边 界，没 有 流 体 能 贯 穿 边 界，因 此 切 向 力 设 为，法向速度设为；流道内壁面无滑移，因此法向速度和切向速度设为，螺杆转速设为 。，张一明等：不同螺槽深度下异向双螺杆挤出机流道仿真分析仿真结果分析压力场如图所示，因为异向双螺杆结构复杂性，所以啮合区出 现 局 部 高 压 情 况，这 种 情 况 被 称 为 压 力 突 变。种螺槽深度不同的螺杆压力分布普遍不均匀，螺棱处压力较高，螺槽处压力较低。为了方便分析种螺杆对应的流道的压力分布情况，在螺杆啮合区沿挤出方向建立压力参考轴线，得出轴向距离和压力的关系图如图所示。随着螺槽深度的变浅，压力波动越大，对物料的延压有利，同时提高了双螺杆挤出机的建压能力。剪切速率场由图可知，种螺槽深度不同的螺杆剪切速率分布不均匀，在螺杆啮合区和机筒内壁的间隙处出现最高剪切速率。由于剪切速率梯度的存在，聚乳酸熔体可以充图不同螺槽深度下的压力云图 图轴向压力折线图 分地混合。为了方便分析种螺杆对应的流道的剪切速率分布情况，在螺杆啮合区沿挤出方向建立剪切速率参考轴线，得出轴向距离和剪切速率的关系图如图所示。剪切速率的峰值点出现在螺棱位置，随着螺槽深度的变浅，剪切速率波动变大，剪切作用更强，分散混合效果更好。模型在螺棱处出现突变值，说明此处产生漏流现象。黏度场如图所示，根据“剪切变稀”现象，即在流动速度较大处剪切作用较大而黏度较小；在流动速度较小处剪切作用较小而黏度较大。通过对比图和图，发现符合此现象。为了方便分析种螺杆对应的流道黏度分布情况，在螺杆啮合区沿挤出方向建立黏度参考轴线，得出轴向距离和黏度的关系图如图 所示。图不同螺槽深度下的剪切速率云图 图轴向剪切速率折线图 食品装备与智能制造 总第 期 年月图不同螺槽深度下的黏度云图 图 轴向黏度折线图 停留时间分布采用示踪粒子轨迹跟踪法，在流道入口布置 个示踪粒子，利用 模块计算分析示踪粒子在流道内的运动轨迹，累积停留时间分布曲线如图 所示。种模型的累积停留时间分布曲线趋势大体一致，当概率密度为时，随着螺槽深度的变浅，停留时间逐渐变短。停留时间分布曲线如图 所示，螺槽越深，拐点之间的距离越大，曲线总长越长，轴向混合能力越好，但自清洁能力减低。结果与分析试验条件异向双螺杆挤出机：选择，种转速，进料速度为 。图 累积停留时间分布曲线 图 停留时间分布曲线 注塑机：融融温度 ，注塑温度 。旋转流变仪：平板模式，频率为 ，温度为 ，应变为。试验结果分析复合黏度由图 可知，随着螺杆转速的升高，复合黏度随频率的增大逐渐减小，符合“剪切变稀”现象，即转动速度较大时剪切作用较大而黏度较小；转动速度较小时剪切作用较小而黏度较大。聚合物熔体的流动性图 不同螺杆转速下 材料的复合黏度 ，张一明等：不同螺槽深度下异向双螺杆挤出机流道仿真分析取决于分子链段的运动能力，螺杆转速的增加改变 的分子结构，使其分子链断裂，降低其黏度。储能模量与损耗模量由图、图 可知，随着螺杆转速的升高，的犌与犌均随之巨幅减小。这是由于螺杆转速增加，中分子链相互断裂，进一步提升了分子链段的自由运动，表现为犌与犌明显减小。图 不同螺杆转速下 材料的储能模量 图 不同螺杆转速下 材料的损耗模量 结论研究运用 流体仿真软件进行数值模拟，借助 软件对后处理结果进行数值分析；对种不同螺槽深度状态下的异向双螺杆挤出机的聚乳酸流体的压力场、剪切速率场和黏度场，以及停留时间分布曲线 进行了对比研究。结果表明：通过改变螺杆元件的齿根圆直径及其对应的流道直径来实现螺槽深度的改变，在其他条件相同的情况下，螺槽深度越浅，压力波动越大，剪切速率波动越大，但黏度波动越小；螺槽深度越深，轴向混合性能越高。仿真挤出过程中，随着螺槽深度的变浅，停留时间逐渐变短。螺槽越深，拐点之间的距离越大，曲线总长越长，轴向混合能力越好，但自清洁能力降低。随着螺杆转速的增加，聚乳酸材料的复合黏度逐渐降低，储能模量与损耗模量均呈现“剪切变稀”现象，符合仿真规律。后续应改善螺杆模型，对比与分析仿真结果，同时通过大量试验设计出更适于聚乳酸加工及改性的螺杆构型。参考文献?刘玲?任国庆?潘鑫?等?聚乳酸生产技术及市场分析?化学工业?（?）?（?）?耿孝正?双螺杆挤出机及其应用?中国塑料?（?）?（?）?范涛?陈作炳?杜佳佳?等?螺槽深度对单螺杆挤出机熔融输送段熔融过程的影响?塑料工业?（?）?（?）?张雪晨?张丽梅?黄志刚?等?啮合同向双螺杆挤出机不同螺杆构型对其流场的影响?食品与机械?（?）?（?）?（?）?彭涛?贾朝阳?同向双螺杆不同螺纹元件混炼效果的数值研究?中国塑料?（?）?（?）?许国智?粒子示踪法研究胶体悬浮液在复合外场下的动态?广州?华南理工大学?耿孝正?双螺杆挤出机及其应用?北京?中国轻工业出版社?陈龙?董力群?张亚军?等?基于?的异向双螺杆的参数化设计?中国塑料?（?）?（?）?（下转第 页）食品装备与智能制造 总第 期 年月?（?）?刘魁立?非物质文化遗产及其保护的整体性原则?广西师范学院学报?（?）?（?）?张玉霞?钧瓷系与禹州钧窑?中原文物?（?）?（?）?阎夫立?阎飞?王双华?中国钧瓷?郑州?河南科学技术出版社?王文章?第二批国家级非物质文化遗产名录图典?北京?文化艺术出版社?杨国政?杨永超?中国钧瓷工艺?郑州?黄河水利出版社?商亚敏?钧瓷的审美特征与创新性设计研究?陶瓷 学 报?（?）?（?）?董?钧瓷?北京?北京工艺美术出版社?邱族周?钱惠?刘云?基于闽东传统民居特色的非遗食品包装设计?食品与机械?（?）?（?）?褚兴彪?雍自高?唐田甜?等?地理茶叶品牌标志设计中字形应用分析?湖南包装?（?）?（?）?左文?民俗文化要素在地方特色食品包装设计中的表达?食品与机械?（?）?（?）?张玉山?裴金秀?基于“绿色设计”理念下的陶瓷食品包装容器设计要素探究?湖南包装?（?）?（?）?刘春媛?王春鹏?徐丽平?意境在现代家居产品设计中的魅力?包装工程?（?）?（?）?周雅琴?穆政臣?基于互动型设计的食品包装创新研究?包装工程?（?）?（?）?（上接第 页）?王少峰?张媛?刘美莲?等?不同温度状态下啮合异向双螺杆挤出机的流道分布规律?食品与机械?（?）?（?）?韩海川?刘谨?李开国?等?双螺杆啮合元件与反螺纹元件对流场影响的仿真?工程塑料应用?（?）?（?）?蒋卫鑫?黄志刚?赵玉莲?等?基于?田字形网格食品托盘挤出成型分析?食品与机械?（?）?（?）?黄志刚?蒋卫鑫?李鑫?等?啮合异向双螺杆挤出机仿真研究?包装学报?（?）?（?）?（?）?杨明山?赵明?塑料成型加工工艺设备?北京?印刷工业出版社?雷雁洲?王少伟?吕秦牛?等?碳纤维?聚乳酸复合材料的结晶性能和流变特性?复合材料学报?（?）?（?）?包装与设计 总第 期 年月