基于混合元的三维注塑充填数值模拟-材料加工工程专业论文.docxVIP

Ill Ill l l Tl l l ll ll l ll I I IIl {Y1 928535 原创性声明 本人郑重声明：所呈交的学位论文，是本人在导师的指导下，独立进行研究 所取得的成果。除文中已经注明引用的内容外，本论文不包含任何其他个人或集 体已经发表或撰写过的科研成果。对本文的研究作出重要贡献的个人和集体，均 已在文中以明确方式标明。本声明的法律责任由本人承担。 学位论文作者： R期： 年 月 同 使用授权声明 本人在导师指导下完成的论文及相关的职务作品，知识产权归属郑州大学。 根据郑少I,I大学有关保留、使用学位论文的规定，同意学校保留或向国家有关部门 或机构送交论文的复印件和电子版，允许论文被查阅和借阅；本人授权郑州大学 可以将本学位论文的全部或部分编入有关数据库进行检索，可以采用影印、缩印 或者其他复制手段保存论文和汇编本学位论文。本人离校后发表、使用学位论文 或与该学位论文直接相关的学术论文或成果时，第一署名单位仍然为郑州大学。 保密论文在解密后应遵：了：此规定。 学位论文作者： 日期： 年 月 同 v 摘 摘 要 摘 要 充填过程对高分子聚合物材料注塑产品的最终质量具有重要影响。在过去 的几十年中，以有限元为代表的计算机数值技术在注塑成型领域的应用研究已 经经历了从一维到二维再到三维的过程。目自i『三维注翅充填模拟面临的主要障 碍是解的稳定性和求解速度问题。本文针对该充填过程，研究并选择了一套基 于三维实体模型的求解方案，独立编写了一套具有完全自主知识产权的计算机 软件，对该充填过程的速度、压力以及温度场进行了模拟。其中所做的主要工 作有： 1、分析了采用P1／PO单元(速度线性、压力常数)离散求解流场时，由经 典Galerkin法得到的有限元方程无法求解的原因在于不可压缩条件造成的过约 束问题。进而提出了一套基于P1／PO宏元的求解方案：对二维和三维速度场采用 线性三角形和四面体进行离散，而对压力场的离散则打破原有的单元界限，由 若干相邻单元组成一个个所谓的宏元，在每一个宏元内部压力取常数。经过这 样处理可以适当地放开过约束条件，从而可以使得原来基于P1／PO元的无法求解 的有限元方程变得稳定可解。该方案可以降低待解方程的自由度数量，从而能 够提高计算机的求解速度。 2、分析了在采用有限元／有限差分法求解注塑充填瞬态温度场时，对于固定 的网格，若时I’日j步长选择的不合理，容易使得方程的稳定性要求得不到满足而 导致求解失败。采用了分步法，将温度场方程分解成为一个对流方程和一个传 导方程，然后采用不同的时间步长独立进行求解，从根本上解决了由对流项导 致的数值问题。对瞬念项时间域的离散采用向后差分格式，利用集中质量热容 矩阵代替协调质量热容矩阵，较好地抑制了由于时间步长过小而导致的瞬态热 传导方程求解过程中出现的数值问题。该温度场求解方案的物理意义明确，得 到的结果合理。 3、对于含有不同几何尺度的注塑制品，特别是含有微特征结构的注塑制品， 采用统一的数值方法进行模拟时，计算机求解速度和计算精度的矛盾比较突出。 鉴于此，本文对小尺度部分采用P2／P1单元离散，以获得较高的计算精度，而对 大尺度部分则采用P1／PO宏元离散，以获得较快的整体求解速度，在两种单元之 摘要间采用过渡元进行连接，以满足速度场的连续性要求。该模拟方案兼顾了求解 摘要 间采用过渡元进行连接，以满足速度场的连续性要求。该模拟方案兼顾了求解 速度和计算精度，实现了两者的有机统一。 4、应用本文的求解方案计算了变厚度制品以及含有加强筋制品沿截面方向 的速度场和压力场分布；模拟了在熔接线形成过程中两股熔体沿截面的熔接行 为；模拟了使用“H’’型平衡流道的一模四腔制品充填过程中的不平衡充填现象； 采用多种混合元模拟了含有微结构注塑制品的充填过程；分别采用轴对称单元 和口【{面体单元模拟了三维气辅注塑成型中的聚合物充填和高压气体的一次穿透 行为。 关键词：注塑有限元P1／PO宏元三维流动模拟 「 「 AbstractAbstract Abstract Abstract Filling stage of the polymer injection molding plays an important role on the quality of the final product．In the last several decades，the numerical analysis technology,especially finite element method(FEM)，has been applied to the simulation of the plastic injection molding．The researc