液态金属充型过程三维流动场数值模拟.pdfVIP

第 32卷 第5期 华 侨 大 学 学 报 (自然 科 学 版 ) Vol_32 No．5 2011年 9月 JournalofHuaqiaoUniversity(NaturalScience) Sep．2011 文章编号： 1000-5013(2011)05—0481—04 液态金属充型过程三维流动场数值模拟 刘晶峰，李洪友，江开勇 (华侨大学机电及 自动化学院，福建 泉州 362021) 摘要：基于有限差分法建立液态金属充型过程流动计算的数学模型，使用SOIA—VOF数值模拟技术开发 液态金属充型过程三维流动场数值模拟分析软件．然后 ，利用该软件计算标准实验铸件充型过程的三维流动 场，结果与经典算例近似，表明该软件关于流场的计算精度较高． 关键词： 充型；流动场；有限差分法；SOLA-VOF算法 ；数值模拟 中图分类号： TG244；TP391．9 文献标志码： A 液态金属成形的充型过程对铸件质量的影响很大，各种缺陷如冷隔、浇不足等都是在液态金属充型 不利的情况下产生的．因此，研究液态金属的充型过程对获得健全铸件有着重要意义．液态金属的充 型过程是在不透明铸型中进行的，且具有速度快、温度高的特点，研究液态金属的流动规律大都借助于 间接的方法．在充型过程数值模拟研究中，主要采用 SIMPLE算法、MAC—SMAC算法及 SOLA—VOF 算法等几种方法 ．其中，SOLA—VOF算法是 目前应用最为广泛的方法，国外一些著名的铸造过程模拟 软件，如Magmasoft，Procast等均采用SOLA-VOF算法 ．应用这种数值计算方法，可以对极其复杂的 液态金属成形过程进行定量描述，达到改进浇注系统、优化铸造工艺的目的．本文介绍利用 自主开发的 金属液态成形工艺分析系统，进行液态金属充型过程三维流动场数值分析的全过程． 1 数理模型 铸件充型过程中，液态金属的流动为粘性不可压缩流体带有 自由表面的非稳态流动，而且是三维 的，它的运动状态可用质量守恒方程 (连续性方程)和动量守恒方程 (Navier—Stokes方程)来表示 ．在三 维直角坐标系中，其表达式分别为 + +_dW = 0， (1) ： dY d u + )，．02u_十 O2u_十 O2u． Ou+w O 1OPq _ + 一 一 3x - v． ， d d2 … 瓦十 十d d2一一 … 十十 ]J +十 +十 +十叫叫笔一一一一吉一Oy+十gy+十y．82v__十十82v_十十82v．，，} ((22)) +十 +十舅+十碧瓦一一一一1篑瓦+十gz：+十yyc磐+十—Oy—2+十磐．jJ 式 (1)，(2)中：P为流体的密度；7为流体的运动粘度；P为流场中(，，z)点的压力； ，，W为(，Y，) 点的流速在3个坐标轴方向的速度分量 gx，g，g为重力加速度在3个坐标轴方向的分量． 2 计算方法 采用SOLA—VOF方法叫计算流动域及其速度、压力分布规律．速度、压力按交错网格方式离散，进 收稿 日期：2010—11—18 通信作者 ：刘晶峰(1964一)，男，副教授，主要从事材料成形CAD／CAE及计算机图形学的研究．E-mail：ULo27@163． com ． 基金项 目： 国家自然科学基金资助项目(50675