论文摘要
高分子材料、木材、金属和硅酸盐并称世界四大材料体系,是国民经济和国防建设中重要的生产材料。围绕高分子材料,目前已形成了门类齐全的高分子材料加工工业体系并在国民经济中占有重要地位。挤出是高分子材料成型加工中一类重要工艺,通过挤出模具能够模塑所有的热塑性材料和某些热固性材料,可以生产管材、板材、棒材、网材、单丝、薄膜、异型材、发泡型材、多种材料复合制品以及线缆等带包覆层的工业制品。挤出成型过程中,受温度、压强、应力以及作用时间等变化的影响,体系中聚合物熔体的聚集态结构和化学结构会发生变化。挤出工艺条件在很大程度上决定着材料的结构和性能,并最终影响成型制品的外观和质量。由于挤出成型过程影响因素较多,借助常规实验方法,成本高且耗时费力。目前,围绕聚合物挤出流动过程所开展的实验研究大多在实验室中完成,其实验目的主要是从研究材料自身流变性出发,以简单剪切和拉伸流动为测试模型。由于先进测试方法难以直接引入挤出工艺现场,普通测试方法又不可避免会对挤出加工流场中聚合物熔体流动行为造成影响并导致大量系统误差,因此,实验方法只能定性但难以定量描述聚合物加工中的实际流变行为。数值模拟技术经过近几十年的发展,其对于大规模复杂物理和工程问题的优秀解算能力吸引了科研人员的目光,并逐步在力学、传热学、材料学等诸多领域得到广泛应用,加快了现代科技发展步伐。本文探讨了数值模拟技术在聚合物流变学中的应用,研究了聚合物挤出成型过程的模型化理论和数值计算方法,构造了能够描述挤出过程中聚合物流变行为特点的物理和数学模型。通过数值模拟技术,成功预测挤出加工流场中聚合物熔体的流动速度、应力和温度等重要场变量分布特点及其变化规律,讨论并分析聚合物的复杂流变行为及其成型机理。将成型过程数值模拟技术与现代优化设计理论相结合,建立并运用相应的优化模型和算法,实现成型工艺与制品质量的优化控制。在聚合物流变学基础上,结合流体动力学理论,建立了基于Euler描述的非牛顿流体三维非等温流动数学模型。采用基于低阶插值的罚函数有限元方法,成功避免了混合有限元方法中为稳定压力项求解所采用的高阶插值,从而有效利用运算资源,将求解空间扩展至三维。充分考虑了剪切速率和温度变化对材料流动性的影响,采用非线性粘度模型,实现了流动与传热的耦合求解。通过线性化交替迭代算法,在迭代过程中实现非线性项的线性化,减小了初始变量分布对计算收敛性影响。采用流线迎风Petrov-Galerkin（SUPG）方法,通过构造非对称权函数来增大流场中来流上游效应,克服了标准伽辽金（Galerkin）方法在处理对流扩散方程时因对流占优所引起的数值振荡问题。通过理论分析,推导了幂率流体圆管泊肃叶（Poiseuille）流动中的各场变量分布函数,将模拟与理论结果进行比较,以验证该非牛顿流动数学模型与数值计算方法的可靠性。本文所建立的非牛顿流体三维非等温流动数值模拟技术不仅可用于分析一般非牛顿流动问题,在处理目前流动模拟技术难以预测的复杂工艺问题时也表现出较强的适应性。基于该方法,本文首次针对一类新型复合共挤出工艺——异型材钢塑共挤工艺过程开展了数值建模及其工艺分析工作。根据其不同于常规挤出工艺的特点,建立了该工艺过程中聚合物成型的数学模型,通过模拟计算,得到聚合物熔体由测向导入至复合成型整个流动过程的速度、温度、应力和压力等各场变量的分布,掌握了聚合物熔体的流动特点与成型规律,讨论了体积流量和钢衬移动速度变化对各场变量分布的影响。定义流速分布相对均差作为判断挤出流动平衡性的依据,模拟并得到不同导入角角度和分流段长度对挤出流动平衡的影响,分析结论可为实际异型材钢塑共挤工艺设计提供理论指导与支持。挤出加工流场中的聚合物熔体除了具有非牛顿流动特性外,还表现出较强的弹性流动特点。针对该问题,本文建立了能够真实反映成型过程中聚合物熔体粘弹流变特性的三维粘弹非等温流动数学模型,构造了稳定的有限元数值求解方法并将其成功应用于聚合物挤出工艺模拟与分析中。采用PTT（Phan-Thien—Tanner）本构模型描述聚合物熔体的粘弹流变行为,在反映聚合物熔体剪切流动特点的同时,能较真实地反映其拉伸流动特点。考虑粘弹介质特有的能量耗散模式,在热力学第一定律基础上,根据非平衡不可逆热力学理论首次推导了该粘弹介质的能量守恒方程。将粘弹性附加应力张量作为有限元基本解,采用解耦方法实现了三维空间中速度场、温度场和流动应力场的多物理场稳定求解。将应力张量作为拟体力项处理后,动量方程会失去椭圆性并导致计算结果发散。通过分离粘弹分裂（DEVSS）方法,引入稳定化因子对动量方程进行椭圆化处理,提高了速度场求解的稳定性。采用非协调流线迎风（SU）方法克服了本构方程在对流占优时的数值振荡问题,实现了应力场的稳定求解。通过对粘弹流体4:4:1收缩流动模拟结果与polyflow软件模拟结果的比较,证明本文所建立的粘弹流动模型和数值计算方法的合理与可靠性。将该模型与方法成功应用于异型材挤出工艺过程模拟,讨论了网格密度、罚数和能量分割系数等计算控制参数对模拟结果的影响。分析了中空异型材挤出过程中聚合物熔体的流动速度、温度和应力分布,讨论了加工流场中聚合物熔体的粘弹流动特点,获得了工艺条件和口模结构参数变化对聚合物流变行为的影响规律。挤出胀大是聚合物挤出成型工艺中无法回避的一个问题。由于加工中聚合物熔体的粘弹流变特性,熔体离开口模时,形变回复等会导致聚合物熔体的挤出胀大,表现为挤出物截面形状和尺寸发生变化,对挤出制品的尺寸和精度造成影响。本文在粘弹流动数值模拟技术基础上,建立了聚合物熔体三维挤出胀大的数学模型和数值求解方法并编制了相应的有限元模拟程序。针对一种工业用低密度聚乙烯（LDPE）的挤出胀大问题开展了实验及其数值模拟研究。通过控制应变流变仪分别得到小幅振荡剪切流动中储能和耗能模量的分布以及稳态剪切流动中剪切粘度和第一法向应力差的分布。采用非线性回归方法拟合流变测量实验结果,得到以PTT本构模型表征的材料线性和非线性粘弹流变参数。通过间接测量,得到不同螺杆转速时LDPE通过圆形口模时的出口挤出胀大比。采用本文所建立的聚合物挤出胀大数学模型及其数值求解方法,对实验条件下LDPE的挤出胀大过程进行模拟,比较了挤出胀大比的实验和模拟结果。通过挤出胀大数值模拟,进一步讨论了LDPE通过圆环口模时的出口挤出胀大问题,得到实验中难以测得的流动速度和应力等场变量分布,可对挤出胀大特点及其形成机理进行定性与定量分析和预测。聚合物挤出成型过程数值模拟是被动式的,实际应用中需依靠专业工程技术人员的智力、知识和经验,对计算结果进行分析、评价,然后修改设计。本文将挤出成型过程数值模拟技术与优化设计理论相结合以实现成型过程与制品质量的优化控制。在成型过程模拟技术基础上,提出了一种基于数值模拟、前馈神经网络和遗传算法的聚合物挤出工艺与模具优化设计方法。根据挤出流动平衡原则,建立了以出口流动均匀性为目标,以工艺和模具结构参数为设计变量的优化模型。通过成型过程数值模拟获得目标函数值以建立训练神经网络模型的样本库,采用反向误差传播算法进行网络学习,建立用于预测隐目标函数的神经网络模型,从而有效减小有限元模拟计算量。通过遗传算法与神经网络的交互运算,得到优化结果,使设计建立在科学分析的基础上从而提高挤出加工工艺设计水平。探讨了该优化设计方法各模块计算实施中的关键技术问题,编制了聚合物挤出工艺与模具优化设计程序。分别针对异型材钢塑共挤工艺过程和片材挤出工艺过程进行优化设计并达到相应的优化目标。数值模拟技术在聚合物加工工程领域的应用已成为计算流变学的重要分支之一。近年来,尽管流体力学数值方法取得了较大进展,但在处理复杂流动问题时,其解算能力仍受到流动区域和计算稳定性的限制。尤其对于复杂工程与工艺问题,其数学建模和数值模拟关键技术研究鲜有报道。本文针对挤出加工流场中聚合物熔体非牛顿粘性和弹性流动特点所建立的数学模型以及所构造的稳定数值计算方法,对于丰富计算流变学具有一定的理论意义。针对挤出工艺过程中聚合物熔体的复杂流变行为及其成型机理所作的分析与研究,以及基于此所开展的成型工艺模拟与优化工作具有较大的工程应用价值。