不同收缩口模中聚合物熔体流变行为的数值模拟.pdf

D O I:1 0.3 9 6 9/j.i s s n.1 0 0 3-0 9 7 2.2 0 2 3.0 3.0 1 2 文章编号:1 0 0 3-0 9 7 2(2 0 2 3)0 3-0 4 1 5-0 5不同收缩口模中聚合物熔体流变行为的数值模拟刘庆升*,刘友琼,柯春海(信阳师范大学 数学与统计学院,河南 信阳 4 6 4 0 0 0)摘 要:挤出成型是聚合物加工工业中的一种重要成型方法。在挤出成型过程中,聚合物熔体在收缩口附近容易出现应力集中现象,从而导致最终制件的表面缺陷。以最终消除或减小表面缺陷为目的,基于数值模拟的方法,主要研究了口模的收缩外形对聚合物熔体流变行为的影响,包括具有不同收缩长度的V-型收缩口以及喇叭型收缩口。采用R o l i e-P o l y模型刻画聚苯乙烯熔体在41收缩流中的应力-应变关系,基于有限体积法及S I MP L E算法模拟了不同收缩长度条件下,聚苯乙烯熔体在V-型收缩口模以及喇叭型收缩口模中的主应力差、第一法向应力差、伸展率和速率。数值结果表明:V-型收缩口、喇叭型收缩口以及收缩长度的增大,都能有效改善收缩流中的应力集中现象,在收缩长度相同的条件下,喇叭型收缩口模具有更好的效果。关键词:R o l i e-P o l y模型;41收缩流;有限体积法中图分类号:O 2 4 2.1;O 3 5 7.1 文献标识码:A开放科学(资源服务)标识码(O S I D):N u m e r i c a l S i m u l a t i o n o f R h e o l o g i c a l B e h a v i o r s o f P o l y m e r M e l t i n D i f f e r e n t C o n t r a c t i o n D i e sL I U Q i n g s h e n g*,L I U Y o u q i o n g,K E C h u n h a i(C o l l e g e o f M a t h e m a t i c s a n d S t a t i s t i c s,X i n y a n g N o r m a l U n i v e r s i t y,X i n y a n g 4 6 4 0 0 0,C h i n a)A b s t r a c t:E x t r u s i o n m o l d i n g i s a n i m p o r t a n t m e t h o d i n p o l y m e r p r o c e s s i n g i n d u s t r y.T h e s t r e s s c o n c e n t r a t i o n o f p o l y m e r m e l t c a n o c c u r e a s i l y a t t h e n e a r b y r e g i o n o f c o n t r a c t i o n e n t r a n c e i n t h e p r o c e s s o f e x t r u s i o n m o l d i n g,w h i c h c a u s e s v o l u m e d e f e c t s o f t h e f i n a l p a r t s.F o r t h e p u r p o s e t o e l i m i n a t e o r r e d u c e t h e v o l u m e d e f e c t s,t h e e f f e c t s o f c o n t r a c t i o n p r o f i l e s o n t h e r h e o l o g i c a l b e h a v i o r s o f p o l y m e r m e l t a r e m a i n l y i n v e s t i g a t e d b a s e d o n n u m e r i c a l m e t h o d,i n c l u d i n g V-s h a p e d a n d t r u m p e t-s h a p e d c o n t r a c t i o n p r o f i l e s a t d i f f e r e n t c o n t r a c t i o n l e n g t h.R o l i e-P o l y m o d e l i s u s e d t o d e s c r i b e t h e s t r e s s-s t r a i n r e l a t i o n s h i p o f p o l y s t y r e n e m e l t i n 41 c o n t r a c t i o n f l o w.N u m e r i c a l s i m u l a t i o n s o f t h e p r i n c i p a l s t r e s s d i f f e r e n c e,t h e f i r s t n o r m a l s t r e s s d i f f e r e n c e,s t r e t c h r a t i o a n d v e l o c i t y i n V-s h a p e d a n d t r u m p e t-s h a p e d c o n t r a c t i o n f l o w s a r e i m p l e m e n t e d b a s e d o n f i n i t e v o l u m e m e t h o d a n d S I MP L E a l g o r i t h m.N u m e r i c a l r e s u l t s s h o w t h a t t h e V-s h a p e d,t r u m p e t-s h a p e d c o n t r a c t i o n p r o f i l e s a n d t h e i n c r e a s e o f c o n t r a c t i o n l e n g t h c a n r e d u c e t h e s t r e s s c o n c e n t r a t i o n,t h e t r u m p e t-s h a p e d c o n t r a c t i o n p r o f i l e s h a v e b e t t e r r e s u l t s u n d e r t h e s a m e c o n t r a c t i o n l e n g t h c o m p a r e d w i t h V-s h a p e d c o n t r a c t i o n p r o f i l e s.K e y w o r d s:R o l i e-P o l y m o d e l;41 c o n t r a c t i o n f l o w;f i n i t e v o l u m e m e t h o d0 引言随着黏弹本构方程和模拟技术的发展,准确刻画聚合物熔体的流变特性成为可能。数值模拟作为一种描述和预测黏弹流体流变性的重要方法,成为聚合物加工工业中的一种重要的设计辅助工具。一直以来,人们都希望能够对复杂流体尤其是聚合物熔体的非线性动力学行为进行准确刻画。为此发展了一系列的黏弹本构方程,包括P T T模型1、W a g n e r模型2、D o i-E d w a r d s模型3、P o m-P o m模型4、R o l i e-P o l y模型5以及相应的改进模型6,并将这些本构模型应用于非牛顿流体力学的研究7-9。上述模型各有优缺点,根据聚合物分子链的特征,R o l i e-P o l y模型适宜描述线型聚合物熔体的流变特征,而P o m-P o m模型适宜描述分支型聚合物熔体的流变特征。聚合物在复杂的加工成型条件下常牵扯到熔体的流动与变形,这与收缩流中复杂流体的流动类 收稿日期:2 0 2 2-0 6-0 6;修回日期:2 0 2 2-0 9-2 5;*.通信联系人,E-m a i l:q i n g s h e n g 4 0 81 6 3.c o m 基金项目:国家自然科学基金项目(1 1 9 0 1 5 0 4);河南省自然科学基金项目(2 0 2 3 0 0 4 1 0 3 4 0)作者简介:刘庆升(1 9 8 2),男,河南驻马店人,讲师,博士,主要从事科学与工程计算、复杂流体的多尺度建模和仿真等方向的研究。514信阳师范学院学报(自然科学版)J o u r n a l o f X i n y a n g N o r m a l U n i v e r s i t y第3 6卷 第3期 2 0 2 3年7月 N a t u r a l S c i e n c e E d i t i o n V o l.3 6 N o.3 J u l.2 0 2 3似。因此,收缩流常作为基准流动用于研究聚合物熔体的流变行为1 0-1 3。HA S S E L L等1 4基于实验研究了不同类型的聚苯乙烯熔体在收缩扩张型腔中的流变行为,详细讨论了收缩口的曲率、熔体的流率对流动不稳定性的影响。C OMB E AUD等1 5通过挤出成型实验研究了口模的收缩角对流动双折射的影响,并得出了制件表面缺陷受到收缩角强烈影响的结论。AGA S S ANT等1 6详细综述了聚合物熔体挤出成型过程的流动不稳定的实验后发现,这些挤出不稳定是影响挤出生产线最大输出量的关键因素;作者还详细总结了影响收缩流中挤出不稳定的可能因素,并提出了消除或尽可能减少这些挤出不稳定的方法。AGA S S ANT等1 7简要综述了聚合物加工在工业、实验研究及模型研究的发展现状,同时指出了精确地预测聚合物熔体的流变行为一直是一个挑战。在挤出成型中,聚合物熔体在收缩口附近往往会产生应力集中现象,该现象会导致熔体破裂以及最终制件出现鲨鱼皮型的表面缺陷。本文的主要目的是基于有限体积法和R o l i e-P o l y本构模型研究聚苯乙烯熔体在收缩口模中的流变行为,探究在不同收缩长度的条件下V-型收缩口模及喇叭型收缩口模对聚苯乙烯熔体流变性的影响,同时给出流场中有关主应力差、剪切力、伸展率等更多更详细的信息,能够更好地理解聚合物熔体进入收缩口时的动力学响应过程。1 数学模型1.1 质量和动量守恒方程基于不可压缩及等温条件,对聚苯乙烯熔体流变性进行研究,同时不考虑表面张力的影响。因此,质量和动量守恒方程可以写成如下的形式:肀u=0,(1)(ut+u肀u)=-肀p+肀+s肀2u,(2)式中:u、分别表示速度矢量和应力张量,、p及s分别表示熔体的密度、压力及溶剂的黏度。1.2 R o l i e-P o l y本构方程R o l i e-P o l y本构方程,常用该模型来描述线型聚合物熔体的应力-应变状态,其形式为ddt=肀u+肀uT-1d(-I)-2(1-(3/t r)R(+*(t r 3)(-I),(3)式中:d和R分别表示聚合物链的取向和松弛时间,I和t r分别表示单位矩阵和张量的迹,*是与聚合 物 链 的 约 束 松 弛 有 关 的 参 数,经 常 取 为-0.5,伸展率定义为=(t r/3)。1.3 动量方程及R o l i e-P o l y本构方程的无量纲形式数值结果将采用无量纲形式来表示,定义U和L分别为速度和长度尺度,则动量方程和R o l i e-P o l y本构方程的无量纲化可通过式(4)得到X*=XL,u*=uU,t*=t UL,*=LpU。(4)动量方程的无量纲化形式为:R e(ut+u肀u)=-肀p+肀+肀2u,(5)式中:R e=U L/表示雷诺数,=s+p为总黏度,=s/(s+p)为溶剂黏度与总黏度的比值,p表示聚合物黏度。R o l i e-P o l y本构方程的无量纲化形式为:W et+u肀-肀u -肀u T =1-肀u+肀u T -21-W eRI1-3/-1+2W eW eR +2W e 3/W eR1-*+6W e*W eR,(6)式中:=W e/1-t r+t r I,W e=dU/L和W eR=RU/L分别为基于链的取向和伸展松弛时间所定义的W e i s s e n b e r g数。1.4 边界条件速度矢量u和应力张量的入口边界条件为u=ui n f l o w,(7)=0。(8)出口边界采用齐次N e u m a n n边界条件:un=0,(9)n=0,(1 0)式中:n表示出口的熔体的流出方向。固壁采用无滑移边界条件,相应的应力张量可以由方程(6)推导得到,可参阅文献1 2,在此不再赘述。1.5 流动诱导双折射聚合物熔体的流动诱导双折射可以通过偏光弹性仪得到,这种现象产生的原因是应力的各向异性引起了光折射率的各向异性。实验观测到的双614第3 6卷 第3期信阳师范学院学报(自然科学版)h t t p:/j o u r n a l.x y n u.e d u.c n2 0 2 3年7月折射与模拟得到的主应力满足应力光学定律,即n=c,(1 1)式中:n表示双折射;表示主应力差;c为常数,又称为应力光学系数。对于二维情况,的形式为=(x x-y y)2+42x y。(1 2)2 数值方法基于同位网格有限体积法对质量守恒方程、动量守恒方程和R o l i e-P o l y本构方程进行离散,并采用S I MP L E算法1 8进行求解。在控制体内对质量守恒方程进行积分,可得到其离散形式(ue-uw)y+(vn-vs)x=0,(1 3)式中:uw和ue分别表示控制单元左侧和右侧界面上x方向的流率,vn和vs分别表示控制单元上侧和下侧界面上y方向的流率,x和y分别表示x和y方向上的网格步长。以同样方式对动量方程进行离散,对流项采用中心差分格式,可得aPuP=aEuE+aWuW+aNuN+aSuS+S(u),(1 4)式中:S为动量方程的源项,系数aP、aE、aW、aN、aS与对流项和扩散项密切相关,其形式为aE=DeAPe +m a x-Fe,0 ,aW=DwAPw +m a x-Fw,0 ,aN=DnA(|Pn|)+m a x(Fn,0),aS=DsA(|Ps|)+m a x(Fs,0),aP=aE+aW+aN+aS+xy/t,式中:De、Dw、Dn、Ds为单元面上的扩散系数,Pe、Pw、Pn、Ps为单元面上的P e c l e t数,Fe、Fw、Fn、Fs为单元面上的通量,对流项的离散采用中心差分格式A(|P|)=1-0.5|P|。上述系数的表达式为Fe=uey,De=y/R e(xE-xP),Pe=Fe/De,Fw=uwy,Dw=y/R e(xP-xW),Pw=Fw/Dw,Fn=vnx,Dn=x/R e(yN-yP),Pn=Fn/Dn,Fs=vsx,Ds=x/R e(yP-yS),Ps=Fs/Ds。本构方程的离散形式可类似得到。3 数值结果所采用的聚合物材料为巴斯夫公司提供的聚苯乙烯(代号为P S 2)。数值模拟基于单模式的R o l i e-P o l y本构模型描述聚苯乙烯熔体在不同收缩口模流动过程中的应力-应变关系,熔体温度设为2 0 0。首先,基于41收缩流验证模型及算法的有效性;接着,讨论在不同的收缩长度L下,V-型收缩口模及喇叭型收缩口模对聚苯乙烯熔体流动的影响。当收缩长度L=0时,此时型腔为41突然收缩型腔,如图1(a)所示。图1 41收缩口模示意图(a)突然收缩口模;(b)V-型口模;(c)喇叭型口模F i g.1 S c h e m a t i c d i a g r a m o f 41 c o n t r a c t i o n d i e s(a)A b r u p t c o n t r a c t i o n;(b)V-s h a p e d c o n t r a c t i o n;(c)T r u m p e t-s h a p e d c o n t r a c t i o n当收缩口为V-型口或喇叭口时,如图1(b)和图1(c)所示,为了保证数值算法的顺利实施,采用区域扩充法1 8将计算区域扩充为图1(a)所示区域。计算区域的直角坐标系如图1所示,计算过程中采用的网格为均匀网格x=y=0.0 2 5。本文算法基于V i s u a l C+语言在3.6 GH z的酷睿i 7-4 7 9 0处理器上实现。3.1 模型及算法验证收缩流常用于挤出成型和注射成型等聚合物加工领域,高剪切和高拉伸的流动分别出现在模壁附近及收缩口附近区域,收缩流能更好地反映聚合物流体的物理特性,而41突然收缩流常作为基准算例。选用文献1 0 中的算例验证模型及算法的有效性以及R o l i e-P o l y本构方程刻画聚合物熔体应力-应变关系的能力。主要目的是基于上述模型和算法研究收缩口的不同外形对聚苯乙烯熔体流变行为的影响,下面仅简要给出验证性结果。图2为不同流率下聚苯乙烯熔体在41突然收缩口模中的双折射实验及模拟结果对比图,其中图2(a)(c)分别对应的流率为1.4 3 3 3、2.8 6 6 7、4.3 mm/s,图的下半部分为实验结果1 0,上半部分为本文的模拟结果。从图2可以看出,当流率为1.4 3 3 3 mm/s时,所预测的条纹数与实验是相同的,而在高流率时,预测出的条纹数比实验测的条纹数略少。众所周知,对黏弹流体流变性的准确预测一直是个难题,为了更准确地描述黏弹流体的流动和变形特性,人们一直在不断地改进已有本构方程及相应算法。总体来说,本文预测的双折射结果与实验结果基本上是一致的,能够反映聚苯乙烯熔体在收缩流714刘庆升,刘友琼,柯春海.不同收缩口模中聚合物熔体流变行为的数值模拟中的流变性,同时说明了本文的模型及相应的求解算法的有效性。图2 不同流率下由实验1 0及模拟所得到的双折射图(a)1.4 3 3 3 mm/s;(b)2.8 6 6 7 mm/s;(c)4.3 mm/sF i g.2 C o m p a r i s o n b e t w e e n e x p e r i m e n t a l a n d p r e d i c t e d b i r e f r i n g e n c e f o r P S 2 m e l t a t d i f f e r e n t f l o w r a t e s(a)1.4 3 3 3 mm/s;(b)2.8 6 6 7 mm/s;(c)4.3 mm/s3.2 V-型收缩口模对聚合物熔体流动的影响本节主要研究V-型收缩口模在不同收缩长度条件下对流动的影响。图3为本文预测的V-型收缩口模在不同收缩长度条件下沿入口和出口中心线上的主应力差、第一法向应力差、伸展率以及速率。图3 本文预测的不同收缩长度的V-型口模沿型腔中心线上的主应力差、第一法向应力差、伸展率和速率F i g.3 T h e p r e s e n t s i m u l a t i o n r e s u l t s o f p r i n c i p a l s t r e s s d i f f e r e n c e,f i r s t n o r m a l s t r e s s,s t r e t c h r a t i o a n d v e l o c i t y a l o n g t h e c e n t e r l i n e o f V-d i e s w i t h d i f f e r e n t c o n t r a c t i o n l e n g t h从图3可以看出,随着收缩长度的变化,这些量在远离收缩口的上游和下游区域并没有发生变化,而在收缩口区域附近发生了明显的变化。主要原因可解释为:在远离收缩口的上游区域,流场的变化趋势是相对平稳的,对聚合物的分子链几乎没有影响,分子链的构型并没有受到扰动,仍处于平衡状态。在离收缩口一定距离的下游区域,流体流动充分发展,聚合物分子链达到了完全松弛状态,相应的物理量不发生变化。而当流体流经收缩口及其附近区域时,分子链受到了较强的剪切和拉伸影响,相应的物理量在此区域变化强烈。随着收缩长度的变大,主应力差、第一法向应力差和伸展率的最大值明显降低,而速率的最大值明显增大,这也是符合理论分析的,众所周知,当流体流经收缩的流道 时,截面积的 减小将会导 致流体速率 的增大。3.3 喇叭型收缩口模对聚合物熔体流动的影响本节主要研究喇叭型收缩口模在不同收缩长度条件下对流动的影响。图4为本文预测的V-型收缩口模在不同收缩长度条件下沿入口和出口中心线上的主应力差、第一法向应力差、伸展率以及速率。图4 本文预测的不同收缩长度的喇叭型口模沿型腔中心线上的主应力差、第一法向应力差、伸展率和速率F i g.4 T h e p r e s e n t s i m u l a t i o n r e s u l t s o f p r i n c i p a l s t r e s s d i f f e r e n c e,f i r s t n o r m a l s t r e s s,s t r e t c h r a t i o a n d v e l o c i t y a l o n g t h e c e n t e r l i n e o f t r u m p e t-s h a p e d d i e s w i t h d i f f e r e n t c o n t r a c t i o n l e n g t h从图4可以看出,主应力差、第一法向应力差、伸展率和速率的变化趋势与V-型收缩口模的情况类似,随着收缩长度的变大,主应力差、第一法向应力差和伸展率的最大值明显降低,而速率的最大值明显增大。从图3和图4容易看出,在相同收缩长度的条件下,喇叭型的收缩口模中聚合物熔体的主应力差、第一法向应力差和伸展率要比V-型口模中的更小,而速率更大。所以相对于突然收缩口,采用V-型口模或喇叭型口模都将会减小收缩口附近的应力集中现象(应力的局部增高),而喇叭型口模的效果更好,能够更大程度地减小应力集中现象,从而能够减小挤出成型中由应力集中导致的制品表面缺陷。814第3 6卷 第3期信阳师范学院学报(自然科学版)h t t p:/j o u r n a l.x y n u.e d u.c n2 0 2 3年7月4 结论主要研究了不同收缩长度下的V-型和喇叭型口模对聚苯乙烯熔体流变性的影响。数值结果表明:(1)本文所采用的数值方法和算法切实有效,R o l i e-P o l y模型能够刻画聚苯乙烯熔体的流变特征,这也是复杂流体可计算建模的一次成功实践;(2)相对于突然收缩流,具有V-型和喇叭型收缩口的型腔能够有效减小流场中的主应力差、第一法向应力差和伸展率,而喇叭型收缩口减小的效果更好;(3)不论是V-型还是喇叭型收缩口模,随着收缩长度的增大,聚苯乙烯熔体在收缩流中的主应力差、第一法向应力差和伸展率都具有明显的下降趋势;(4)在挤出成型中,为尽量减少或避免由应力集中导致的制件表面缺陷,应尽可能选择具有喇叭型外形且具有尽可能长的收缩长度的口模。参考文献:1 TH I E N N P,T ANN E R R I.A n e w c o n s t i t u t i v e e q u a t i o n d e r i v e d f r o m n e t w o r k t h e o r yJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,1 9 7 7,2(4):3 5 3-3 6 5.2 WA GN E R M H,R U B I O P,B A S T I AN H.T h e m o l e c u l a r s t r e s s f u n c t i o n m o d e l f o r p o l y d i s p e r s e p o l y m e r m e l t s w i t h d i s s i p a t i v e c o n v e c t i v e c o n s t r a i n t r e l e a s eJ.J o u r n a l o f R h e o l o g y,2 0 0 1,4 5(6):1 3 8 7-1 4 1 2.3 D O I M,E DWA R D S S F.T h e t h e o r y o f p o l y m e r d y n a m i c sM.N e w Y o r k:O x f o r d U n i v e r s i t y P r e s s,1 9 8 6.4 ME A D D M,L A R S ON R G,D O I M.A m o l e c u l a r t h e o r y f o r f a s t f l o w s o f e n t a n g l e d p o l y m e r sJ.M a c r o m o l e c u l e s,1 9 9 8,3 1(2 2):7 8 9 5-7 9 1 4.5 L I KHTMAN A E,G R AHAM R S.S i m p l e c o n s t i t u t i v e e q u a t i o n f o r l i n e a r p o l y m e r m e l t s d e r i v e d f r o m m o l e c u l a r t h e o r y:R o l i e-P o l y e q u a t i o nJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 0 3,1 1 4(1):1-1 2.6 KA B AN EM I K K,H TU J F.N o n e q u i l i b r i u m s t r e t c h i n g d y n a m i c s o f d i l u t e a n d e n t a n g l e d l i n e a r p o l y m e r s i n e x t e n s i o n a l f l o wJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 0 9,1 6 0(2/3):1 1 3-1 2 1.7 F R AN A H L,O I S H I C M,THOMP S ON R L.N u m e r i c a l i n v e s t i g a t i o n o f s h e a r-t h i n n i n g a n d v i s c o e l a s t i c b i n a r y d r o p l e t c o l l i s i o nJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 2 2,3 0 2:1 0 4 7 5 0.8 G AO P u y a n g,Z HAO Z i t o n g,YAN G Y a n g.T h e n u m e r i c a l m o d e l i n g a n d s t u d y o f g a s e n t r a p m e n t p h e n o m e n o n i n n o n-i s o t h e r m a l p o l y m e r f i l l i n g p r o c e s sJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 2 1,2 9 4:1 0 4 5 7 5.9 A B UG A J G,CH I NYOKA T.B e n c h m a r k s o l u t i o n s o f t h e s t a b i l i z e d c o m p u t a t i o n s o f f l o w s o f f l u i d s g o v e r n e d b y t h e R o l i e-P o l y c o n s t i t u t i v e m o d e lJ.J o u r n a l o f P h y s i c s C o mm u n i c a t i o n s,2 0 2 0,4(1):0 1 5 0 2 4.1 0 T E N CHE V R,G OUGH T,HA R L E N O G,e t a l.T h r e e d i m e n s i o n a l f i n i t e e l e m e n t a n a l y s i s o f t h e f l o w o f p o l y m e r m e l t sJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 1 1,1 6 6(5/6):3 0 7-3 2 0.1 1 O L L E Y P,GOUGH T,S P A R E S R,e t a l.3 D s i m u l a t i o n o f t h e h i e r a r c h i c a l m u l t i-m o d e m o l e c u l a r s t r e s s f u n c t i o n c o n s t i t u t i v e m o d e l i n a n a b r u p t c o n t r a c t i o n f l o wJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 2 2,3 0 4:1 0 4 8 0 6.1 2 L I U Q i n g s h e n g,L I U Y o u q i o n g,J I AN G C h u n t a o,e t a l.N u m e r i c a l s i m u l a t i o n o f v i s c o e l a s t i c f l o w s d u r i n g i n j e c t i o n m o l d f i l l i n g b a s e d o n R o l i e-P o l y m o d e lJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 1 9,2 6 3:1 4 0-1 5 3.1 3 L I U Q i n g s h e n g,L I U G u i x i a n,L I U Y o u q i o n g,e t a l.N u m e r i c a l s t u d y o f t h e f l o w o f p o l y s t y r e n e m e l t s i n c o n t r a c t i o n f l o w u s i n g R o l i e-P o l y m o d e lJ.R h e o l o g i c a A c t a,2 0 2 2,6 1(2):1 3 9-1 6 1.1 4 HA S S E L L D G,EMB E R Y J,MC L E I S H T C B,e t a l.A n e x p e r i m e n t a l e v a l u a t i o n o f t h e f o r m a t i o n o f a n i n s t a b i l i t y i n m o n o d i s p e r s e a n d p o l y d i s p e r s e p o l y s t y r e n e sJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 0 9,1 5 7(1/2):1-1 4.1 5 C OMB E AU D C,D EMAY Y,V E R G N E S B.E x p e r i m e n t a l s t u d y o f t h e v o l u m e d e f e c t s i n p o l y s t y r e n e e x t r u s i o nJ.J o u r n a l o f N o n-N e w t o n i a n F l u i d M e c h a n i c s,2 0 0 4,1 2 1(2/3):1 7 5-1 8 5.1 6 A G A S S AN T J F,A R D A D R,C OMB E AU D C,e t a l.P o l y m e r p r o c e s s i n g e x t r u s i o n i n s t a b i l i t i e s a n d m e t h o d s f o r t h e i r e l i m i n a t i o n o r m i n i m i s a t i o nJ.I n t e r n a t i o n a l P o l y m e r P r o c e s s i n g,2 0 0 6,2 1(3):2 3 9-2 5 5.1 7 A G A S S AN T J F,MA C K L E Y M R.A p e r s o n a l p e r s p e c t i v e o n t h e u s e o f m o d e l l i n g s i m u l a t i o n f o r p o l y m e r m e l t p r o c e s s i n gJ.I n t e r n a t i o n a l P o l y m e r P r o c e s s i n g,2 0 1 5,3 0(1):1 2 1-1 4 0.1 8 陶文铨.数值传热学M.2版.西安:西安交通大学出版社,2 0 0 1.TAO W e n q u a n.N u m e r i c a l h e a t t r a n s f e rM.2 n d e d.X ia n:X ia n J i a o t o n g U n i v e r s i t y P r e s s,2 0 0 1.责任编辑:郭红建914刘庆升,刘友琼,柯春海.不同收缩口模中聚合物熔体流变行为的数值模拟