注塑成型是目前塑料成型的主要方法之一，生产的制品如磁带盒、计算机键盘等。该工艺的流程为，首先将熔融的高聚物熔体以较大的流速注射入模具型腔中，直至充满整个腔体，然后待塑料熔体充分固化后开启模具脱出塑料制品，最后再循环进入下一个制品的生产。 为提高模具的生产效率和水平，人们迫切需要发展针对该工艺过程的高效的数值分析工具来取代传统的费时费力的试模修模的过程。因而，在过去的数十年中，注塑成型过程的数值模拟引起了越来越广泛的关注。尽管如此，仍然还有很多问题有待进一步的研究，如自由面的准确追踪以及稳定、高效、精确地求解控制方程等。本论文的研究工作正是以解决这些问题为目标而展开的。 移动自由面的准确追踪对于注塑成型过程的数值模拟具有十分重要的意义。根据建立方程所采用参考系的不同，目前已有的处理方法大体上可归类为拉格朗日方法、欧拉方法以及任意拉格朗日欧拉方法(ALE)。在ALE方法中，计算网格的运动独立于物质构型的运动，结合了拉格朗日方法和欧拉方法的优点。本文第三章将给出一个基于ALE方法的自由面追踪和网格生成．重生成模型，可在准确确定移动自由面空间位置的同时尽量避免发生网格扭曲。在所发展的模型中，充填域变质量体系的实时网格生成被简化为移动自由面附近区域的多边形三角化过程，大大减少了网格处理所需的计算量。此外，还引入了局部Laplacian光顺方案来提高计算网格的质量，并发展了多种机制来处理自由面与模具边界的接触问题。 实现注塑成型模拟的另外一个关键问题是须发展能稳定求解控制方程初边值问题的数值方法。众所周知，采用标准Galerkin方法模拟不可压缩流动问题通常会导致两种类型的虚假数值振荡。其一来源于控制方程组的可约化混合特性，该特性限制了速度．压力(u-p)插值空间的自由选取。不满足所谓LBB条件的非协调的插值模式，如u-p的等低阶插值模式，将导致压力场的虚假空间振荡。第二种类型的数值振荡是由控制方程的对流特性引起的，尤其易发生在对流占优的情况下，这是由于标准的空间 Galerkin离散方法只对自伴随的算子方程有效。 对于第一个问题，本文通过在有限增量微积分 (FIC)的理论框架下重建质量守恒方程，并理性地引入一个辅助变量来避免空间高阶导数的计算，提出了压力稳定型分步算法(PS-FSA)，有效消除了压力场的虚假数值振荡。与经典的分步算法(C-FSA) 相比，该算法具有更好的压力稳定性。此外，本文还采用这种压力稳定机制改善了经典的成型充填过程中非等温非牛顿粘性流动的ALE有限元与无网格自适应耦合模拟特征线基分裂算法(CBS)的压力稳定性。PS-FSA算法具体的推导过程以及它的数值验证将在第四章中给出。对于第二个问题，本文发展了基于ALE描述的广义特征线Galerkin(CG)方法，有效处理了对流占优问题控制方程的空间离散。而且，与经典的CG方法相比，该方法能使用更大的时间步长。实际上，经典的CG方法以及经典的 Crank-Nicolson (CN)方法能分别作为两个特例而被纳入到本文所提出的广义CG方法的框架中。此外，本文方法还将经典CG方法由欧拉描述下推广到ALE描述下，以便于和前述ALE自由面追踪技术相结合。应指出的是，若将参考系固定于空间不动，ALE描述可退化为欧拉描述，即相对于经典CG方法来说，本文方法更具一般性。将所发展的广义CG方法和分步算法以及第四章给出的压力稳定技术相结合，本文还提出了一个迭代型压力稳定的广义CBS算法，并将其应用于非等温非牛顿流的数值模拟中。该算法具体的推导过程以及它的数值验证将在第六章中给出。 此外，本文还提出了一个有限元与无网格离散区域剖分的自适应算法，可自动地将整个计算域划分成分别应用有限元法(FE)、无网格法(MF)以及它们的耦合插值进行空间离散的三个子域。基于该自适应算法以及连续掺混法(CBM)，本文提出了有限元与无网格的自适应耦合空间离散方法，并将其应用于注塑成型的数值模拟中。所发展的自适应耦合方法可在充分结合有限元法与无网格法各自优点的同时避免它们各自的缺点。该方法具体的实现细节以及展现它相对于传统有限元法和无网格法优势的数值验证将在第五章中给出。 为了论文的完整性，本文第二章还简单介绍了伽辽金型无网格法的基本概念和它所面临的问题以及相应的处理方法。实现本文算法的计算机程序及其数据结构的说明放在第七章。第八章总结全文并展望进一步的研究内容。