塑料挤出成型过程的数值模拟研究
- 格式:pdf
- 大小:320.92 KB
- 文档页数:6
赵艳志等.塑卡辛按出戎型过程的数擅模拟研究
挤出成型过程
的数值摸拟研究
赵艳志,麻向军
(华南理工大学聚合物新型成型装备国家工程研究中心,广东 广卅I 5 l 0640)
摘要:先简单介绍了常用数值计算方法(有限差分法、有限元法、边界元法、有限体积法),并将其运用于塑料挤出成 型过程中。从螺杆挤出过程、口模的设计以及挤出胀大等方面综述了模拟技术在塑料挤出成型加工过程中的研究结果,号召 业内人士主动使用CFD(Computational Fluid Dynamics)软件包进行优化设计,可提高国际竞争力。 关键词:挤出成型;数值模拟;口模;计算流体力学 中图分类号:TQ320.663 文献标识码:B 文章编号:1009—797X(2007)04—0007—06
挤出成型在塑料加工领域中是一种非常重
要的加工手段,塑料的挤出过程是一个连续成
型的过程,物米斗从加米斗斗加入后,在机筒中经
压缩、加热、混合成熔融态流体,被螺杆输送
到挤出流道中,物料在旋转螺杆与料筒之间进
行输送、压缩、熔融塑化。熔体通过塑模(口
模)在压力下成为形状与塑模相似的一个连续
体。塑料在整个挤出机内始终处于运动状态,
具有很强的动态特性。将计算流体力学应用于
塑料成型过程的模拟,可以优化塑料的加工过
程,调控产品的结构和性质,有助于实现加工
的自动化控制,提高劳动生产率。对聚合物流
体流动的研究是复杂流体研究领域中最为活跃
的分支之一。其主要研究内容包括了2个方面:
是如何描述聚合物的流动性质,聚合物的流
动性质是用本构方程描述的,即如何建立其本
构方程;二是当确立其本构方程后,在加工机
械的复杂边界条件下,如何精确求解高度非线
性、非局域的流体力学控制方程。考虑到聚合
物熔体流变性质、螺杆形状及流道几何形状的
复杂性,该数值的计算往往是一件十分艰巨的
工作。 1 常用数值计算方法
常用的数值计算方法包括有限差分方法、
边界元方法、有限体积方法、有限元方法等。这
些方法的基本出发点是~致的,即对整个求解
空间作离散处理,通过划分网格作局部线性化
近似,把求解非线性偏微分方程约化成为求解
每一时刻的线性方程组,之后再依照一定的初
始条件或边界条件进行求解。 1.1 有限差分方法
有限差分法是一种直接将微分问题转变为
代数问题的近似数值解法,曾经是一种在计算
流体中用来求解物体运动时占主导地位的数值
方法。但用此方法求解边界条件复杂,尤其是
椭圆型问题不如有限元方法方便。 1.2 有限元方法
有限元法的基础是极值原理和剖分差值。
它吸收了有限差分法中离散处理的内核,又采
用了变分计算中选择逼近函数并对区域进行积
分的合理方法,具有广泛的适应性,特别适合
作者简介:赵艳志(1968),女,华南理工大学博士研究生 高分子材料加工专业,已发表论文十几篇。 收稿日期:2006—07 1 8
2007 ̄ 第33卷 ・7・
维普资讯 http://www.cqvip.com CHINA RUBBER7PLAs’f ICS TECHNOL0GY AND EQUlPMEN’f ●■■——■■■■——嘲螭豳鞠疆隧鞋耥姚§簿
中的流动作了初步分析和计算,分析了轴向压
力分布对流体流动的影响。同年,BOOY 讨论
了宽浅和窄深两种螺槽中牛顿流体的等温流动,
建立了拖曳流、压力流和啮合区漏流的方程,
分析了螺槽的非充满现象。但聚合物在机筒中
流动受到温度、压力变化的影响,若假设为等
温状态,与实际情况有很大差距。一方面,由
于成型工艺或制品性能的要求,对机筒须加热
或用冷却剂制冷,这使得机筒内各处温度分布
不均匀;另一方面,聚合物熔体在流动过程中
以黏滞性方式形变,所损耗的能量转变为热,
从而使熔体温度上升。而温度的变化会引起聚
合物黏度的显著变化,从而对聚合物的流动产
生影响。实验测定,即使模具壁与芯部温差在
1~2℃,黏度的变化亦十分明显。因此,为更精
确地描述聚合物的流动,聚合物熔体采用非等
温非牛顿模型来模拟。彭炯等" 建立了非牛顿
黏性聚合物熔体在同向旋转双螺杆挤出机中的
三维非等温流动模型,利用流体动力学分析软
件Polyflow,给出了螺杆区内速度场、温度场、
压力场以及黏性耗散热的分布。
(2)随着挤出机的发展,螺杆挤出的数值
模拟逐渐由简单的挤出传递现象模拟转向双螺
杆复杂挤出流动模拟;几何模型由二维发展为
三维流动的模拟。
对双螺杆挤出机所进行的数值模拟研究,
国外一直处于领先地位。1998年,PotenteH 对
幂律流体在啮合同向旋转双螺杆挤出机的流动
行为进行了三维有限元分析,并通过计算所得
的流线来确定停留时间 Rio sI51利用边界元法
(BEM),比较了双螺杆的不同几何构造对混合
效率的影响。Cheng 使用FIDAP模拟了ZSK一
5 3同向旋转双螺杆挤出机内的三维等温流动,
并以长度拉伸比、面积拉伸比以及应变分布等
作为分布混合程度的衡量指标,使用粒子示踪
技术对其分布性和混合情况进行了数值研究。
YaoI l1吏用FIDAP讨论了正反向螺纹输送元件的
不同组合对物料在LCMAX一40机中的分散性混
合情况。2000年,B1.avo[ 】在准静态假设下,对
啮合同向旋转双螺杆挤出机进行了三维等温数 于几何、物理条件复杂的问题,而且便于程序
的标准化,在黏弹流体流动模拟中应用广泛。
1.3 边界元方法
边界元法是将全域的求解化为区域边界上
的求解,简化了计算;有限元法的近似范围在
整个区域上,而边界元法的近似范围仅在区域
边界处,因此边界元法的精度高于有限元法。
边界元法的缺点在于,其方程组的系数矩阵不
对称且为满阵,有时是近似的奇异矩阵,因此
对复杂形状的流场,计算效率较低,而且计算
日寸需要采用解析函数的基本解,所以目前只适
用于线性问题及基本解已知的问题。
1.4 有限体积方法
有限体积法的最大优点是即使在粗网格下
也能得到准确的积分守恒,从而在实际计算中
将大大节省内存的占用,提高计算效率。
总之,以上各种数值解法各有特点,在聚
合物流体数值模拟中都有应用。具体选用何种
数值方法需根据问题的复杂程度决定,如求解
线性问题可以考虑选择简便快捷的有限差分法
和边界元法;对于非线性问题选用有限元方法
和有限体积法,一般来说,对椭圆型问题使用
有限元法,对抛物线问题和双曲线问题则使用
有限体积法。
2 挤出机中的数值模拟
塑料在螺杆挤出机中一般经历固体输送、
压缩、熔融、熔体输送。螺杆是挤出机关键部
件,它影响到产品的产量、挤出熔体出口温度、
熔体质量、混合物混合均匀性等。利用数值模
拟可以预测固体床的分布、沿螺杆物料的温度
和压力、质量流率,以及局部速度、温度、应
力、剪切速率和热传导等。通过对塑料在挤出
机中的数值模拟可以优化螺杆设计,检测挤出
机的螺杆转速、料筒温度分布、材料性能等。挤
出过程的数值模拟发展历程如下:
(1)聚合物模型从等温牛顿流体模型到非
等温非牛顿流体模型。
l980年,Denson等…忽略啮合区的影响,用
有限元法对牛顿流体在不同几何形状螺纹元件
・8・
第33卷 第4期 维普资讯 http://www.cqvip.com 赵铤志等,塑拱按出戎氆过程的数馥模拟研究
值模拟,分析了捏合块的构造以及进出口边界
条件对压力和速度分布的影响。Take Shi1 对
同向旋转双螺杆挤出机TEx 30的捏合盘区域和
异向旋转非啮合连续混合机LCM l 00G的熔体输
送段进行了三维非等温流动的模拟研究,考察
了压力分布和温度分布随螺杆转速及流量的变
化规律。Berzint¨ 考虑重均分子量的变化对反应
体系黏度的影响,模拟了聚丙烯在啮合同向旋
转双螺杆挤出机中由过氧化物引发的可控热降
解。Fukuokat 利用流动分析网络法(FAN),
对聚乙烯接枝乙烯硅在双螺杆挤出机中的接枝
反应过程进行了二维非等温数值模拟研究。
Yo shinagat ]等用有限元方法对紧密啮合同向双
螺杆挤出机啮合段三维流动进行数值模拟,对
双螺杆的分散及分布性混合进行了量化。
M OY S eY等 】用三维离散单元法(Di S C r et e
E1ement Method)模拟了单螺杆固体输送段的
粒子运动情况。在运用数值模拟手段对双螺杆
进行研究的这一领域,与国外相比,国内起步
较晚。在近几年里,北京化工大学的研究者对
双螺杆挤出机在等温情况下的流场做了一些研
究工作【l6 ”。李鹏、耿孝正【16 使用ANSYS对
同向啮合双螺杆的捏合块和螺纹元件流道进行
了三维等温流动的有限元分析。马秀清、耿孝
正等Il9 o 使用ANSYS分别对非啮合异向旋转双
螺杆及啮合同向旋转双螺杆的轴向循环流道进
行了三维等温模拟分析。刘志芳等【2 使用
AN s Y s对新型双锥形双螺杆元件——腰鼓形和
三角形混合块的熔体输送特性进行了计算模拟。
彭炯等使用POlYflOW对同向旋转双螺杆挤出机
在非等温情况下的流场进行了一系列的数值模
划研究【3,22—2 31。
3 挤出口模中的数值模拟
挤出口模的主要作用是将塑料熔体分布于
流道中,使物料以均匀的速度从口模内挤出。
塑料熔体实际分布取决于聚合物的流动性能、
流道几何形状、通过口模的流率以及口模中的
温度分布状况。塑料熔体在挤出口模中的流动
通常是稳态和非等温的,但随着新型挤出机的 出现,挤出产量越来越大,塑料熔体在口模中
的稳定流动易被破坏,塑米斗的临界剪切速率也
易被超越,其结果往往造成型材制品的表面出
现毛糙、螺旋形和波浪形等不规则现象(即鲨
鱼皮症和熔体破碎)。由于挤出量及挤出速度的
提高,型材型坯在定形模和水箱中的停留时间
缩短,因此,在模具设计中,必须精确分析和
计算塑料熔体和型坯在模具中流动、形变和温
度的关系,以及流道的几何形状、操作条件、流
变和热力学性能等。国外近年来已开始大量运
用各种专用的CFD(COmPUtatiOnal Fl LIid
DYnamiC S)软件来对口模进行流动分析和优化
设计。l 985年,BU sh[24】首次将边界元法应用于
牛顿和广义牛顿流体的数值模拟中。近来,以
高斯离散理论为基础的有限体积法在黏弹流体
流动模拟中有一些应用,Sasmall 对K—BKZ流
体的4:1突缩流动进行了数值模拟。但有限体
积法中边界条件的加入较有限元法复杂,因此
仅限于简单几何流场中的黏弹流动模拟。如前
所述,有限元法在黏弹流体流动模拟中也有广
泛的应用。Feng等f26】采用Oldroyd—B本构方程,
用有限元方法开展了聚合物熔体在口模内的二
维动态等温模拟,研究了模壁一粒子、粒子一
粒子相互作用的流动机制。HUY nh[2 7 对牛顿流
体通过环形口模的非等温挤出过程进行了有限
元模拟分析。l 999年,Kihara等【28 对口模内三
维黏弹性流体的流动进行了数值模拟。
Matallah等【29】考虑不同的流道几何形状,考虑
聚合物的弹性,模拟熔体在口模内的流动。
A rda等【30 研究了模具出口形状、模具表面质量
和含氟聚合物添加剂对制品表面质量(主要针
对鲨鱼皮现象)的影响。在有限元求解方法中,