型竺！！些！叁型兰！！旦圣竺！-计算物理.pdfVIP

塾鬟!量!!型!!! 型竺!!些!叁型兰!!旦圣竺!塑!型竺!丝!!!!!竺!!竺竺!!!!旦竺塾些!垫型竺 羔篓』：：鋈 [文章编号】 1001．246x(2002)04—03“一05 注塑模充模过程动态分析的有限元／控制体积法 刘春太， 申长雨， 陈静波， 王利霞 (郑州大学(北校区)檬塑模具国家工程研究中心，河南郑州450002) [摘要]翅料熔体在三维薄壁型腔内的流动属于带有运动边界的粘性不可压流体的流动，基于粘性广义H。l} shaw流动理论，利用控制体积法建市了求解压力场的有限元方程，耦合利用有限元，控制体积法自动跟踪熔体的运 动边界．实现了充模过程的动态模拟． [关键词]注塑成型；充模流动；数值分析；有限元，控制体积法 [中圈分类号]TQ320．5 【文献标识码]A 的扩大而扩大．移动网格法虽能较准确地确定流动 O引言 前沿位置，但在实施过程中必须对时间增量进行特 注塑成型充模流动过程是一个相当复杂的物理 殊处理，以保证计算出的前沿点不落在型腔的边界 过程，非牛顿的高温塑料熔体在压力的驱动下通过 之外，有时甚至需要人工干预．更重要的是，这种方 流道、浇口向较低温度的利腔充填，一方面熔体由于 法很难用于三维制件的流动分析，因为在两个面的 模其传热而快速冷去|J，另一方面因高速剪切产生热 交界处或曲面上，速度方向发生变化，需要许多附加 量，同时伴随有熔体固化、体积收缩、分子取向及可 的几何信息和网格节点信息才能确定前沿位置．固 and 能的结晶过程．随着塑料熔体压力、温度、剪切速率 定网格法主要包括MAc(Markercell)法和FAN (now 等物理量在充模过程中的发展变化，塑料熔体的流 AnalysisNelwork)法1“，其基本算法是，先将整 动行为将惠接影响蠼件的最终形状、内部结构、取向 个型腔划分成网格，再形成相应于各节点的体积单 和残余应力． 元，流出或流人体积单元的节点流量可由节点压力 塑料熔体在二维薄壁型腔内的流动属于带有运 计算出来，最后根据体积单元的充填状况近似确定 动边界的粘性不可压流体的流动，数值分析充填过 流动前沿位置Tadmor最早采用控制体积法概念，利 程将面临3个主要问题：①熔体运动前沿的自动跟 用有限差分法来追踪流动前沿，该方法沿用了FAN 踪；②复杂的几何形状；③流动区域的多连通性．对 方法的基本思想． 本文将结合有限元和控制体积法两种思想的优 运动边界的处理依赖于数值算法(FEM、FDM、BEM) 和阀格系统(固定网格和移动网格)的选择．移动网 点，即利用控制体积法实现对熔体运动前沿的追踪， 格法主要包括网格重划法和网格扩展法，其基本思 而利用有限元处理复杂几何形状的数值分析． 想足根据当前时刻的流动前沿位置、速度分布以及 时间增量，确定下一时刻的流动前沿位置，再对流动 1 充模过程的数学模型 前沿的局部区域划分网格，或者网格随着充填区域 对充模过程，实际的特征值和特征变量见表l 表l充模过程中的特征变量和特征值 Tlblel Va由M怖a耐ch触ristIcⅦll瞄|n咖h喧pr毗蜘 收稿日期】2000一12—26；[修回日期】2001—04—26 基金