基于Fluent管件短射法水辅助注塑数值模拟.docVIP

基于Fluent管件短射法水辅助注塑数值模拟 　　【摘 要】水辅助成型是基于气体辅助成型技术发展起来的，因其充填介质不同，使得水辅助成型相比气体辅助成型有很多的优势。水辅助注塑工艺的主要方法有短射法和溢流法，其主要应用对象是管类制件。鉴于目前短射法水辅助注塑工艺理论研究的极度缺乏，本文对管件的短射法水辅助注塑工艺进行数值模拟研究。 　　【关键词】管件；短射法；水辅助注射成型；数值模拟 　　1.引言 　　水辅助成型是基于气体辅助成型技术发展起来的，因其充填介质不同，使得水辅助成型相比气体辅助成型有很多的优势[1]。水辅助注塑工艺的主要方法有短射法和溢流法，其主要应用对象是管类制件。鉴于目前短射法水辅助注塑工艺理论研究的极度缺乏，本文对管件的短射法水辅助注塑工艺进行数值模拟研究，借助Fluent软件实现水辅助注塑成型流动和冷却过程的数值模拟。 　　2.水辅助注塑成型充填过程的数学模型 　　水辅助注塑成型的充填过程首先是将熔体注入模具型腔，部分填充，而后在高压水的推动下进行完全填充的一个非等温、非稳态、纯粘度的复杂流动问题，精确地描述这一流动过程十分的困难，必须从流体力学的基础理论出发，结合水辅助成型工艺的实际情况进行合理的简化与假设处理，才能建立起较为准确地描述该过程，容易实现的数学模型。 　　2.1 水辅助注塑成型充模流动的理论基础[2] 　　流体流动控制方程是物理学中守恒定律的数学描述。守恒定律包括质量守恒定律、动量守恒定律和能量守恒定律。相对应的流体流动控制方程为连续性方程、动量守恒方程和能量守恒方程。 　　（1）质量守恒方程（连续性方程） 　　（1） 　　式中：表示流体密度；div表示矢量微分算子；u表示流体速度。 　　在笛卡尔坐标系下其展开式为： 　　（2） 　　在柱坐标下其展开式为： 　　（3） 　　（2）动量守恒方程 　　（4） 　　式中：为应力张量；为重力矢量。 　　根据柯西应力张量的描述，其可以分为两部分： 　　（5） 　　其中：p为流体净水压力，为偏应力张量。 　　根据广义牛顿本构方程，偏应力张量表示为： 　　（6） 　　其中：为粘度系数，与温度和压力有关，为应变速率张量。 　　在笛卡尔坐标系下动量方程展开为： 　　（7） 　　其中i=x、y、z 　　在柱坐标系下动量方程展开为： 　　方向： 　　（8） 　　方向： 　　（9） 　　方向： 　　（10） 　　（3）能量守恒方程 　　能量守恒定律表述为：流体微元的能量变化率等于外界对流体微元的加热率加上外界作用力对流体微元所做的功。 　　（11） 　　在笛卡尔坐标系下的展开式为： 　　（12） 　　在柱坐标系下的展开式为： 　　（13） 　　2.2 粘度模型[3] 　　本文采用非牛顿流体幂律模型对聚合物熔体的粘度进行描述。 　　非牛顿流体幂律模型 　　非牛顿流体幂律粘度的形式为： 　　（14） 　　FLUENT允许为幂律函数设置上下限： 　　（15） 　　式中k是流体平均粘度（稠度指数）的量度，n是幂律指数，T0为参数温度，和分别为幂率的上下限。 　　2.3 边界条件 　　水辅助充填过程中，熔体前沿基本上与模腔内气体相接触，其熔体前沿处的压力为零；熔体和介质水在浇口处的边界条件为：熔体和介质水的在浇口处的速度与设定的注塑速度值相等，熔体和介质水的温度在浇口处的温度为设定值。 　　2.4 多相流中界面追踪技术 　　体积函数法是当今移动界面追踪技术中最为有效的方法之一。体积函数法的基本思想是在整个流场中定义一个流体体积函数F，每个网格中该函数定义为一种流体的体积与网格单位体积的比值。体积函数方程可以表示为： 　　（16） 　　对于流体体积函数F，当F=1时表示流体完全充满单元体积网格。当F=0表示流体未填充单元体积，称为空单元。当0F1表示部分填充，称为边界单元。有限体积方法通过求解流体体积函数F，从而实现对移动界面的追踪。 　　3.数值模拟 　　3.1 条件设置 　　（1）模型建立与网格划分 　　本文以一圆管为对象来进行水辅助注射成型的数值模拟研究。该圆管的尺寸为80-8.5（mm）。因为制品为轴对称形状，简化成二维轴对称模型来进行模拟，相对三维模型而言，可大大节省计算量，且模拟结果与三维模拟结果相同。二维轴对称几何模型如图1所示，X轴为对称轴。如图2所示，采用四边形结构化网格对模型进行划分，单元边长0.25mm，单元面17000个。 　　图1 二维模型 　　图2 网格划分 　　（2）材料参数设置 　　材料属性的设定主要是指材料密度、比热、热传导率、粘度及参考温度。水和空气的材料属性在Fluent中可以直接选取，PP材料没有，需要自定义，模