双螺杆多材料3D打印挤出机构设计与仿真分析.pdf

针对目前 3D 打印单一材料不能完全满足应用性能要求的现状,设计出一种双螺杆挤出式 3D 打印多材料的混合打印机构,用于高固含量陶瓷浆料的挤出成型。通过对双螺杆的重要参数进行设计和计算,实现了双螺杆机构的小流率稳定挤出。使用 ANSYS 的 POLYFLOW 模块对螺杆转速 1 15 r/min 下流体的流变特性进行仿真分析,得到整个流道压力场、速度场和剪切速率场的分布规律,仿真结果表明出口流率范围为 0郾 66 6郾 48 mm3/s,打印速度为 2 23 mm/s,适合陶瓷浆料挤出成型。同时对打印机构的混合性能进行评估,得出不同螺杆转速下的材料组分分布,并采用示踪粒子法得到粒子的停留时间分布函数和累积停留时间分布函数;结果表明,螺杆转速越高,其轴向混合能力越弱,自洁能力越强。所设计双螺杆挤出机构可以满足低转速小流率的高精度打印,实现多种材料变组分的实时混合和精准打印。关键词:多材料;3D 打印;双螺杆;陶瓷浆料中图分类号:TP23摇 摇 DOI:10.13543/j.bhxbzr.2023.03.011收稿日期:2022-07-12基金项目:天津市科技计划(18ZXRHGX0020);河北省中央引导地方科技发展资金项目(216Z1804G)第一作者:女,1983 年生,讲师E鄄mail:zhoujing313 引摇 言3D 打印技术又称为增材制造技术,是 20 世纪80 年代后期发展起来的一种新型制造技术1,它是一种集成了机械、控制和计算机技术于一体的复杂的、一体化的成型系统2。随着增材制造技术的快速发展,其在航空航天、交通运输、生物医疗、建筑等方面逐渐得到广泛应用3。近年来,3D 打印技术逐渐向打印多材料零件的方向发展,能用于多材料的3D 打印成型工艺主要有材料挤出、材料喷射和粉末床熔融,而对于材料挤出式的 3D 打印成型工艺,具有良好混合性能的螺杆挤出方式为其提供了一种新的研究路径4。目前,国内外学者对螺杆式 3D 打印机进行了大量研究。Valkenaers 等5的研究表明采用螺杆挤出进料方式可以显著改善材料因长时间受热而发生降解、进料不均匀而挤出不稳定等问题。Li等6研制了螺杆动态混合器,通过实验对各组分的氧化铝和氧化锆梯度元件进行了样品打印、烧结和测试,结果表明挤出成形技术能够打印出功能分级的陶瓷零件。任礼等7分析了传统桌面级熔融沉积制造(fused deposition modeling,FDM)型3D 打印机与螺杆挤出式 3D 打印机的挤出原理,设计了一种新型螺杆挤出式 3D 打印机,可实现小流率高精度打印。丁承君等8改进了陶瓷 3D 打印机螺杆挤出装置的喷头结构,可实现连续供料。陈磊9结合螺杆挤压技术与 FDM 型 Delta 并联式打印机研制了一套小型螺杆挤出装置,利用 STAR鄄CCM+仿真软件对物料在挤压系统中的粒料输送与熔体挤压传动过程进行了数值模拟研究,最终得到了合理的螺杆挤压工艺参数范围。Silveira等10设计了一种基于双螺杆挤出的创新型三维打印头,该打印头专用于允许在过程中混合多种材料并直接沉积产品的 3D 打印机。针对多材料混合打印,现阶段国内外学者对单螺杆挤出的研究较多,可以在一定程度上解决材料应用范围有限、喷头易堵塞的问题,但在打印精度方面仍存在混合效果不佳、挤出连续性差等问题。针对当前所存在问题,本文设计了一种基于双螺杆多材料的 3D 打印挤出机构,不仅可以弥补单螺杆挤压输送能力弱等不足,而且其优良的混合性能、剪切作用等可为解决单螺杆应用受限问题提供有效的解决方案,对在打印过程中需要保持挤出速度稳定,高效率打印出表面平整、纹理清晰、成型精度较高的零件具有较为实际的意义。本文分析了双螺杆挤出机构的运动原理,对螺杆进行了结构设计,最后通过软件仿真对其流场进行了分析以及混合性能的评估。1摇 3D 打印机挤出机构设计1郾 1摇 整体结构挤出机构是整个 3D 打印机的关键部件。本文设计的双螺杆多材料 3D 打印机挤出机构如图 1所示。图 1摇 挤出机构结构Fig.1摇 Structure of the extrusion mechanism摇为使双螺杆能实现同步转动,以步进电机驱动齿轮转动,从而带动双螺杆实现同步旋转。双料筒的设计用于存储两种材料以实现多材料打印,料筒上方的接口与气泵相连,气泵产生的气压能将陶瓷浆料通过输料导管顺利流入料腔中。导管接口与导管连接处易于拆卸,紧固性和密封性都能达到挤出要求。料腔与螺杆之间使用的密封轴承能保证螺杆旋转的稳定性,同时也能防止在料腔浆料压力过大时从料腔上方溢出。浆料进入料腔后,由于螺杆的挤压输送,浆料会在料腔中不断地进行挤压、剪切、拉伸和混合,最后将混合好的浆料通过挤出喷嘴稳定挤出。该 3D 打印挤出机构具有正常的挤出功能,且结构紧凑、易拆卸,对于多材料的陶瓷浆料挤出打印具有兼容性并降低了其打印成本。1郾 2摇 双螺杆结构双螺杆螺纹元件的型线从本质上决定了挤出机构的剪切、混合、输送等性能,螺杆端面型线的设计对提高双螺杆挤出机构的整体性能具有重要意义。本文基于相对运动法设计啮合型同向双螺杆螺纹元件端面型线,其原理图如图 2 所示,其中 O1、O2为圆心,B1 B4为圆上的点,A 为两圆盘中心距。图 2摇 相对运动法原理图Fig.2摇 Schematic diagram of the relative motion method摇根据微流挤出 3D 打印机的特性,螺杆不宜太大,在满足挤出流量的条件下,应合理选用螺杆结构参数。本文初步确定螺纹头数 Z=2,中心距A=5 mm,螺杆外径 Do=6 mm,则螺杆顶径 Ro=Do/2=3 mm,螺杆根径 Ri=A-Ro=2 mm,则螺杆内径Di=4 mm,根据经验公式在螺纹头数为 2 时,其中心距比Cr需满足式(1),计算得 Cr=A/2Ro=0郾9逸 2/2,满足条件。再根据式(2)(5)11分别计算出螺顶角琢=22郾 8毅,螺槽角 茁=22郾 8毅,螺侧角 渍=33郾 6毅,半啮合角 准=33郾 6毅。摇 摇 Cr=A/2Ro逸cos(90毅/Z)(1)摇 摇 琢=仔/Z-2准(2)摇 摇 茁=琢(3)摇 摇 渍=cos-1(A/2Ro)(4)摇 摇 准=渍(5)摇 摇 根据同向啮合型双螺杆相对运动的形成过程,可以得出螺纹元件的端面型线方程(式(6)。39第 3 期摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 周摇 婧等:双螺杆多材料 3D 打印挤出机构设计与仿真分析摇 摇籽(兹)=Ri,-琢2臆兹臆琢2或 仔-琢2臆兹臆仔+琢2-Ro(cos兹-琢)2+A2-R2o sin(2兹-琢)2,琢2臆兹臆仔2-茁2Ro,仔2-茁2臆兹臆仔2+茁2或3仔2-茁2臆兹臆3仔2+茁2-Ro(cos仔-兹-琢)2+A2-R2o sin(2仔兹-琢)2,仔2+茁2臆兹臆仔-琢2-Ro(cos兹-仔-琢)2+A2-R2o sin(2兹-仔-琢)2,仔+琢2臆兹臆3仔2-茁2-Ro(cos 2仔-兹-琢)2+A2-R2o sin(22仔-兹-琢)2,3仔2+茁2臆兹臆2仔-琢2(6)摇 摇 将式(2)(5)计算得出的数值代入式(6)中,可以得到双螺杆的端面型线图如图 3 所示,之后再通过 SolidWorks 三维建模软件建立三维模型。螺杆设计为三段式,包括输送段(L1)、混合段(L2)和计量段(L3):输送段螺距较大,输送能力强,将刚进入料腔的陶瓷浆料沿着挤出方向输送;混合段螺距最小,提升剪切拉伸作用,使得陶瓷浆料得到充分混合;计量段螺距适中,能有效稳定浆料从打印喷嘴挤出。为了防止螺杆与料腔内壁的磨损和碰撞,保持螺杆与螺杆之间、螺杆与料腔内壁之间的间隙(分别为 啄1、啄2)非常重要。双螺杆和料腔总体结构如图 4 所示,螺杆各参数及数值如表 1 所示。图 3摇 双螺杆端面型线Fig.3摇 Profile diagram of the twin鄄screw图 4摇 双螺杆与料腔结构Fig.4摇 Structure of the twin鄄screw and material cavity表 1摇 螺杆各参数及数值Table 1摇 Screw parameters and dimensions参数数值螺纹头数 Z2螺杆内径 Di/mm4螺杆外径 Do/mm6中心距 A/mm5螺顶角 琢、螺槽角 茁/(毅)16郾 26螺侧角 渍、半啮合角 准/(毅)36郾 87输送段长度 L1/mm36混合段长度 L2/mm15计量段长度 L3/mm18螺杆与螺杆间隙 啄1/mm0郾 4螺杆与料腔内壁间隙 啄2/mm0郾 2喷嘴直径 De/mm0郾 62摇 双螺杆区域流体仿真分析2郾 1摇 模型参数与边界条件的设定本文选用可打印的高固含量陶瓷浆料作为研究对象。高固含量陶瓷浆料是一种高黏度的非牛顿流体,浆料黏度由式(7)计算得出。摇 摇 滋=子酌(7)式中,滋 为黏度,Pa s,子 为剪切应力,Pa,酌为剪切速率,s-1。不同流体的剪切应力随剪切速率的变化如图 5 所示。高固含量陶瓷浆料属于屈服假塑性流体,由图可看出浆料黏度随剪切速率的提高而降低,呈现剪切变稀的现象,其本构方程遵循 Herschel鄄Bulkley 模型,表达式为1249北京化工大学学报(自然科学版)摇 摇 摇 摇 摇 摇摇摇 摇 摇 摇 摇 摇 摇 摇 2023 年摇 摇 滋=子0酌+(K酌酌)cn-1(8)式中,子0为屈服应力,Pa,K 为黏度系数,Pasn,n 为流动系数,酌c为临界剪切速率,s-1。本文选择比较适合打印的固含量分别为 54%和 58%的氧化铝及氧化锆陶瓷浆料作为 3D 打印材料,两种浆料的物性参数如表 2 所示13。图 5摇 流体剪切应力随剪切速率的变化Fig.5摇 Fluid shear stress as a function of shear rate摇表 2摇 陶瓷浆料物性参数Table 2摇 Physical parameters of the ceramic slurry浆料子0/PaK/(Pa sn)n酌c/s-154%Al2O316417郾 40郾 58158%ZrO226718郾 60郾 681摇 摇 考虑流体的性质以及双螺杆运动网格划分难度较大,选用 ANSYS 中的 POLYFLOW 模块进行流体仿真分析。POLYFLOW 是一个主要用于黏性流体的动力学仿真 CFD 模块,拥有较多的流体计算模型,其独特的网格重叠技术(mesh superposition tech鄄nology,MST)有效避免了重新划分网格的过程,并具备较高的精度14。流场分析过程中,根据具体情况给出如下假设:1)浆料充满整个流道且不可压缩;2)不计重力、惯性力等远小于黏性力的体积力;3)料腔内壁和螺杆表面无滑移;4)流场中浆料温度不变。同时设定流体域的边界条件:每个入口压力为5 000 Pa;左右两个内孔为无滑移边界,无流体能够穿透该边界流动;外壁面为无滑移边界;出口为outflow;设定两个螺杆为运动部件,螺杆转速为1 15 r/min。2郾 2摇 流体区域网格模型的划分双螺杆挤出过程中料腔与螺杆之间的间隙对浆料的剪切有很大影响。为了描述间隙内的物料流动,需要对间隙局部进行合理的网格划分。在流体外壁设定 4 层膨胀层,每层宽度为 0郾 05 mm,总宽度与料腔和螺杆之间的间隙一致。同时在浆料入口、混合段区域和喷嘴流体位置施加局部网格细化命令。其余流体域与螺杆的网格单元尺寸均设定为0郾 3 mm,且均用四面体网格划分,最终得到网格数2 498 338、节点数 608 915 的网格模型。单元网格质量评估平均值为 0郾 775 81,评估值大于 0郾 6,满足仿真计算的要求。划分好的网格模型如图 6 所示。图 6摇 流体域与螺杆网格划分Fig.6摇 Fluid domain and screw meshing摇2郾 3摇 压力场分析5 r/min 螺杆转速下的压力云图和各转速下压力沿挤出方向的变化趋势如图 7 所示。由图 7 可以看出沿着浆料挤出方向,输送段整体压力呈现增大的趋势,这是由于输送段到混合段螺杆的螺距是逐渐减小的,螺距的减小使得输送能力减弱,浆料在挤出过程中不断向结构密集的流道汇集,造成压力增大;混合段压力又呈现出降低的趋势,这是由于混合段螺距增大,其流道空间逐渐变大,起到一定的释压效果;计量段压力呈现小幅度上升,这是由于在螺杆尾部区域流道逐渐变小,浆料在此处不断被挤压,使得压力上升;最后当浆料脱离有螺杆存在的流道后一直到从喷嘴挤出,其压力逐渐减小到 0,这是由于随着挤出头的不断挤出,流道压力得到了释放。整个流道压力的大体趋势是升-降-升-降。同时从图中也可以看出,在有螺杆存在的区域,每隔一段距离都会有一个较小的先上升后下降的现象,这是浆料在每一个螺距内所呈现的压力变化,越靠近螺棱处压力越大,螺槽处压力最小。另外从图 7 也可以看出随着螺杆转速的提升,流道整体的压力也是逐渐上升的,59第 3 期摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 周摇 婧等:双螺杆多材料 3D 打印挤出机构设计与仿真分析说明螺杆转速上升对螺杆建压能力有着正向作用。在螺杆各转速下,最大压力比入口压力提高3 4 倍,表明本文所设计的三段式双螺杆挤出结构建压能力强,能够将浆料稳定地挤出进行打印。图 7摇 螺杆转速 5 r/min 下的压力云图与各转速下的压力变化对比Fig.7摇 The pressure cloud diagram and the comparison ofpressure change at each speed when thescrew speed was 5 r/min摇图 8摇 螺杆转速为 5 r/min 时的流场速度Fig.8摇 Flow field speed for a screw rotationspeed of 5 r/min2郾 4摇 速度场分析图 8 为螺杆转速 5 r/min 下的流体在料腔中形成的迹线、无螺杆区域速度云图及出口截面速度云图。从迹线图中可以看出浆料在随着螺杆转动的同时也被螺杆挤压沿着挤出方向运动,且越靠近出口速度越大。从出口截面速度云图可以看出,越靠近出口浆料挤出的中心速度越大,在出口边缘处速度很小,这是由浆料黏度高、与流道内壁面间的阻力系数大造成的。对出口截面速度进行数据提取,计算出口挤出速度的平均值,得出当螺杆转速为5 r/min时,双螺杆挤出机构的浆料挤出速度为8郾 53 mm/s。图 9 为不同螺杆转速下的出口挤出速度与出口流率曲线,其中出口流率由流体仿真结果得出,出口理论 速 度 通 过 出 口 流 率 除 以 出 口 截 面 面 积(0郾 28 mm2)计算得出,出口平均速度通过提取出口截面速度数据再求平均值得出。从图中可以看出出口理论速度、出口平均速度两条曲线基本吻合。对不同螺杆转速进行速度仿真分析,结果表明出口流率、出口速度与螺杆转速均成正比,出口挤出速度范围为 2 23 mm/s,对比市面上陶瓷 3D 打印机的挤出速度(1 50 mm/s),其能够满足陶瓷浆料挤出3D 打印机的要求。图 9摇 螺杆转速为 1 15 r/min 下的出口挤出速度与流率Fig.9摇 The exit extrusion speed and flow rate when thescrew rotation speed was 1-15 r/min摇2郾 5摇 剪切速率场分析图 10 给出了不同螺杆转速下料腔内剪切速率的分布状态。从图中可以看出螺纹尖端的剪切速率大,螺槽处的剪切速率小,说明螺纹尖端的剪切作用比螺槽处强;在螺纹尖端陶瓷浆料受到更大的拉伸作用,特别是在啮合区域,浆料会受到向里与向外两个剪切力,所以此处的剪切速率最大,更有利于浆料的混合。以螺杆转速 5 r/min 下的仿真结果为例,其螺纹尖端的剪切速率由式(9)计算得出,其中 N 为螺杆转速,计算结果为7郾 85 s-1,螺纹尖端处最大剪切速率仿真结果为 7郾 88 s-1,计算结果与仿真结果相近。通过仿真实验发现,不同转速的螺纹尖端剪切速率最大值均出现在Z=69北京化工大学学报(自然科学版)摇 摇 摇 摇 摇 摇摇摇 摇 摇 摇 摇 摇 摇 摇 2023 年图 10摇 螺杆转速为 5 r/min 时的剪切速率云图与各转速下的剪切速率变化对比Fig.10摇 Shear rate cloud diagram at 5 r/min and comparisonof shear rate variation for different screw rotationspeeds摇13 mm 处,这是由于此时浆料从入口随着螺杆输送到此处的填充率是最高的。螺杆转速提高,螺纹尖端处的剪切速率也随之增加,且沿挤出方向剪切速率呈以升-降-升为周期的变化趋势,两峰值之间的距离与螺距相对应,表明螺棱处剪切速率最大,螺槽底部剪切速率最小,本文所设计的双螺杆挤出机构剪切速率大小适宜。摇 摇 酌=2仔DoN2啄2(9)3摇 混合性能分析3郾 1摇 组分分布由于双螺杆挤出机构啮合区域存在高剪切速率,所以在双螺杆挤出过程中,输送和混合是同步进行的。当料腔流道趋于稳定时,不同螺杆转速下两种陶瓷浆料在流道中的组分分布情况如图 11 所示,两种陶瓷浆料同时从两个入口进入,进入速率比为1颐 1。随着螺杆的转动,左右螺杆通道的浆料会通过啮合区域进行材料交换,组分云图中显示 Al2O3陶瓷浆料的分布情况,红色代表含量为 100%,蓝色代表含量为 0,绿色部分代表浆料组分占比 50%,认为此时浆料混合完成。从图中可以看出,随着螺杆转速的提高,浆料混合完成所需要的时间缩短。在螺杆转速为1 r/min 时,螺杆输送能力与混合能力均较差,在螺杆的计量段浆料未完成完全混合,但是从图中也可以看出当浆料输送到混合段时,浆料在高剪切速率的作用下,其分散频率明显增加。在螺杆转速为 5 r/min 时,浆料通过混合段的高频分散之后,浆料混合基本完成。对比 Z=70 mm 处 3 种转速的组分分布情况,发现在螺杆转速 10 r/min 时浆料混合效果最好,5 r/min 时基本混合完成,1 r/min 时混合未完成,说明转速太低不适用于混合打印。组分云图在一定程度上反映出浆料在流道中的分布情况,为进行实物实验提供了一定的数据支撑。螺杆混合段的设计在低转速下起着非常大的分散剪切作用,在高转速下也对浆料的混合效率发挥着积极作用。图 11摇 不同转速下流道浆料组分分布Fig.11摇 Distribution of slurry components in theflow passage at different rotation speeds摇3郾 2摇 混合性能评估流体在流道内混合情况的优劣,不仅受物料承受的剪切拉伸应力影响,还与物料在流道内停留的时间有关。因此在评价混合能力时,必须考虑流体在流道中停留的时间。采用示踪粒子分析法(parti鄄cal tracking analysis,PTA)计算浆料在螺杆流道中的停留时间。在入口处设置进入 2 000 个粒子,运动时长为 1 200 s,得出粒子的停留时间分布函数(RTD)和累积停留时间分布函数(CRTD),分别由概率密度函数和概率函数表示15。停留时间分布反映了轴向混合的能力16,其不仅可以反映混合性能及输送行为,还可以结合局部微观分析为研究挤出材料性能提供依据17。停留时间分布的函数表达式为79第 3 期摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 周摇 婧等:双螺杆多材料 3D 打印挤出机构设计与仿真分析摇 摇 E(t)=c(t)乙肄0c(t)dt艿c(t)移肄0c(t)驻t(10)式中,c(t)为示踪粒子的浓度。累积停留时间分布为 0 到 t 时间内对停留时间分布的积分,其表达式为摇 摇 F(t)=乙t0E(t)dt艿移t0c(t)驻t移肄0c(t)驻t(11)根据仿真结果统计,不同螺杆转速下的停留时间分布函数和累积停留时间分布函数如图 12、13 所示。示踪粒子最快流出流道的时间为图像与水平轴的第 一 个 交 点,用 t0表 示,t10则 表 示 转 速 为1 r/min 时的 流 出 时 间,得 出 t10=670郾 6 s、t50=357郾 9 s、t100=173郾 7 s、t150=114郾 4 s,可以发现随着转速的增大,t0是逐渐减小的。通过仿真结果计算得出各转速下的平均停留时间,分别为 t11=956郾 3 s、t51=580郾 2 s、t101=292郾 5 s、t151=206郾 4 s,t1越大代表轴向的混合性能越强,以上结果表明随着转速升高,轴向混合能力呈下降的趋势。图 12摇 不同螺杆转速下的停留时间分布函数Fig.12摇 Residence time distribution functions atdifferent screw speeds摇累积停留时间分布函数可以代表在某一时刻有多少粒子已经流出流道。用 t2表示有 90%的粒子已经流出通道的时间,不同转速下有 t12=1 035郾 8 s、t52=760郾 8 s、t102=402郾 5 s、t152=299郾 6 s,t2越小代表自洁作用越强,以上数据表明转速越高,双螺杆挤出机构的自洁效果越好。螺杆转速越高,粒子流出流道的时间越短,轴向混合能力越弱,自洁效果越好。仿真数据为实际的打印实验提供一定的数据支撑,然而在实际打印零件过程中要达到最佳的打印效果,还需在后续的打印实验中进一步探究。图 13摇 不同螺杆转速下的累积停留时间分布函数Fig.13摇 Cumulative residence time distributionfunctions at different screw speeds摇4摇 结论(1)对双螺杆挤出式多材料 3D 打印挤出机构进行了结构设计,并分析其运动原理;重点对双螺杆的端面型线进行了设计和计算,确定了三段式变螺距螺杆的功能段划分及螺杆与螺杆的间隙(啄1=0郾 4 mm)、螺杆与料腔内壁的间隙(啄2=0郾 2 mm)等关键参数。(2)针对高固含量陶瓷浆料,对所设计的挤出机构进行流体仿真分析。通过对压力场、速度场和剪切速率场进行分析,得出螺杆转速1 15 r/min 下的流道压力、速度和剪切速率沿轴向的分布趋势,结果表明出口流率范围 0郾 66 6郾 48 mm3/s,打印速度范围 2 23 mm/s 为适宜的挤出成型参数。(3)对混合性能进行评估,得出不同螺杆转速下的组分分布,并采用示踪粒子法得到粒子的停留时间分布函数和累积停留时间分布函数,结果表明转速越高,粒子流出流道的时间越短,螺杆的轴向混合能力越弱,自洁效果越好,可为多材料实时混合打印实验提供有效的数据参考。参考文献:1摇 卢秉恒.增材制造技术-现状与未来J.中国机械工程,2020,31(1):19-23.LU B H.Additive manufacturing鄄current situation and fu鄄tureJ.China Mechanical Engineering,2020,31(1):19-23.(in Chinese)2摇 KRAKHMATOVA V Y,ZAKHAROV A I,ANDREEV DV,et al.Methods of additive technologies for the manu鄄facture of ceramic products J.Glass and Ceramics,2019,75(11-12):479-484.89北京化工大学学报(自然科学版)摇 摇 摇 摇 摇 摇摇摇 摇 摇 摇 摇 摇 摇 摇 2023 年3摇 NGO T D,KASHANI A,IMBALZANO G,et al.Addi鄄tive manufacturing(3D printing):a review of materials,methods,applications and challenges J.CompositesPart B:Engineering,2018,143:172-196.4摇 田静,薛平,孙华,等.微型螺杆挤出式3D 打印机的研究进展J.塑料,2018,47(1):46-51.TIAN J,XUE P,SUN H,et al.Research developmentof micro鄄screw extrusion 3D printerJ.Plastics,2018,47(1):46-51.(in Chinese)5摇 VALKENAERS H,VOGELER F,FERRARIS E,et al.A novel approach to additive manufacturing:screw extru鄄sion 3D鄄printingC 椅Proceedings of the 10th Interna鄄tional Conference on Multi鄄material Micro Manufacture.San Sebasti佗n,2013:235-238.6摇 LI W B,ARMANI A,MARTIN A,et al.Extrusion鄄based additive manufacturing of functionally graded ce鄄ramicsJ.Journal of the European Ceramic Society,2021,41(3):2049-2057.7摇 任礼,白海清,鲍骏,等.螺杆挤出式3D 打印机结构设计与仿真分析J.中国塑料,2021,35(4):98-105.REN L,BAI H Q,BAO J,et al.Structural design andsimulation analysis of screw extrusion 3D printer J.China Plastics,2021,35(4):98-105.(in Chinese)8摇 丁承君,吴畏,朱智辉.陶瓷 3D 打印机喷头结构设计与仿真分析J.天津工业大学学报,2017,36(4):84-88.DING C J,WU W,ZHU Z H.Structure design and sim鄄ulation analysis of ceramic 3D printer nozzleJ.Journalof Tianjin Polytechnic University,2017,36(4):84-88.(in Chinese)9摇 陈磊.基于数值模拟的塑料颗粒 3D 打印机关键技术研究D.哈尔滨:哈尔滨工业大学,2016.CHEN L.The research of key technologies of plastic par鄄ticles 3D printer based on simulationD.Harbin:Har鄄bin Institute of Technology,2016.(in Chinese)10 JUSTINO NETTO J M,SILVEIRA Z D C.Design of aninnovative three鄄dimensional print head based on twin鄄screw extrusion J.Journal of Mechanical Design,2018,140(12):125002.11 魏静,张光辉.双螺杆捏合机螺纹元件型线设计M.北京:科学出版社,2018.WEI J,ZHANG G H.Design theory and method of pro鄄file of screw element for twin鄄screw kneaderM.Bei鄄jing:Science Press,2018.(in Chinese)12 李耀刚,叶晓濛,纪宏超,等.先驱体陶瓷材料3D 打印机螺杆挤出装置的设计与优化J.北京工业大学学报,2019,45(12):1173-1180.LI Y G,YE X M,JI H C,et al.Design and optimizationof the screw extrusion device for a precursor ceramic ma鄄terial 3D printer J.Journal of Beijing University ofTechnology,2019,45(12):1173-1180.(in Chinese)13 夏晓光,段国林.基于溶剂挥发的陶瓷直写成型工艺J.硅酸盐学报,2020,48(12):1880-1885.XIA X G,DUAN G L.Ceramic direct ink writing processbased on solvent evaporationJ.Journal of the ChineseCeramic Society,2020,48(12):1880-1885.(in Chi鄄nese)14 陈晋南,胡冬冬,彭炯.Polyflow 软件包在聚合物挤出成型中的应用 J.世界科技研究与发展,2002,24(1):28-34.CHEN J N,HU D D,PENG J.Applications of Polyflowin extrusion processJ.World Sci鄄Tech R&D,2002,24(1):28-34.(in Chinese)15 秦望.螺杆挤出式 3D 打印机结构设计及其数值模拟研究D.汉中:陕西理工大学,2021.QIN W.Structure design and numerical simulation ofscrew extrusion 3D printerD.Hanzhong:Shaanxi Uni鄄versity of Technology,2021.(in Chinese)16 张先明.挤出过程停留时间分布的实验研究和数值模拟D.杭州:浙江大学,2008.ZHANG X M.Studies on residence time distribution intwin鄄screw extruders using experimental method and nu鄄merical simulation D.Hangzhou:Zhejiang University,2008.(in Chinese)17 ZIEGLER G R,AGUILAR C A.Residence time distribu鄄tion in a co鄄rotating,twin鄄screw continuous mixer by thestep change methodJ.Journal of Food Engineering,2003,59(2-3):161-167.99第 3 期摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 周摇 婧等:双螺杆多材料 3D 打印挤出机构设计与仿真分析Design and simulation analysis of the extrusion mechanism ofa twin鄄screw multi鄄material 3D printerZHOU Jing1摇 HE HouXiang1摇 LI Ya1摇 DUAN GuoLin2(1.School of Mechanical Engineering,Tianjin University of Science and Technology,Tianjin 300222;2.School of Mechanical Engineering,Hebei University of Technology,Tianjin 300401,China)Abstract:In view of the fact that current 3D printing single materials cannot fully meet the performance require鄄ments,a twin鄄screw extrusion 3D printing multi鄄material hybrid printing mechanism has been designed for the ex鄄trusion molding of a high鄄solid content ceramic slurry.A small flow rate and stable extrusion of the twin鄄screwmechanism was achieved by designing and calculating the important parameters of the twin鄄screw.The POLYFLOWmodule from ANSYS was used to simulate the rheological characteristics of the fluid at a screw speed of 1-15 r/min,and to obtain the distribution of the pressure,velocity and shear rate fields throughout the flow channel.The simu鄄lation results yielded an exit flow rate range of 0郾 66-6郾 48 mm3/s and a print speed of 2-23 mm/s,which is suit鄄able for ceramic paste extrusion.At the same time,the mixing performance was evaluated,and the distribution ofmaterial components at different screw speeds was obtained.The tracer particle method was used to obtain the parti鄄cle residence time distribution function and the cumulative residence time distribution function.The results showthat the higher the screw speed,the weaker the axial mixing capacity and the greater the self鄄cleaning capacity.The twin鄄screw extrusion mechanism can meet the requi