将自然单元法与刚（粘）塑性流动理论相结合，对自然单元法在金属塑性成形过程数值模拟中的应用进行了研究。采用基于Voronoi图和Delaunay三角化结构的Non-Sibsonian插值方法构造近似速度场向量，实现无网格方法中速度边界条件的直接精确施加，提出了基于刚（粘）塑性流动理论的无网格自然单元法。运用不完全广义变分原理，采用罚函数法实现体积不变条件，推导出基于刚（粘）塑性流动理论的无网格自然单元法的离散控制方程，并给出了基于剐（粘）塑性流动理论的自然单元法及其关键算法，拓展了自然单元法的应用范围。典型算例的数值计算结果表明了该方法的可行性和有效性。
刚(粘)塑性流动理论 自然法单元法 Non—sibsonian插值 金属塑性成形 数值模拟