随着计算机技术和数值计算方法的不断发展,刚塑性有限元法在金属体积成形过程中得到了广泛的应用。但是,有限元法在处理金属塑性大变形问题时遇到了困难。有限元法是基于网格的计算方法,当工件变形到一定程度时,有限元网格将发生畸变,此时必须对畸变网格进行网格重划分。网格重划分不仅耗时,而且降低了计算精度。特别地,在三维刚塑性有限元法中实现可靠、高质量的三维网格的自动划分及重划分仍然是一个世界性的难题,成为制约金属成形三维有限元数值模拟技术发展的瓶颈。无网格法是20世纪90年代发展起来的一种新的数值计算方法。由于该方法仅仅采用基于点的近似,不需要节点的连接信息,因此避免了有限元中繁琐的单元网格生成,消除了网格畸变。作为最主要的无网格法之一,再生核质点法(RKPM)不仅具有一般无网格法的特点,还具有其他无网格方法所不具备的优点(如变时-频特性、多分辨率特性等),因而得到国内外学者的广泛研究。本文采用刚(粘)塑性不可压缩材料模型,对再生核质点法(RKPM)在金属二维和三维体积成形过程数值模拟中的应用进行研究,发展基于刚(粘)塑性流动理论的体积成形过程无网格RKPM法数值模拟的基础理论、关键算法和系统开发,扩展无网格RKPM法在金属塑性领域里的应用范围并推动无网格法在体积成形过程模拟分析中的通用化、自动化和实用化。主要研究成果如下:(1)将二维刚(粘)塑性流动理论与无网格RKPM法相结合,提出了二维刚(粘)塑性无网格RKPM法。采用再生核近似构造速度场,分别采用边界奇异核法和修正的罚函数法处理本质边界条件和体积不可压缩条件,采用反正切摩擦模型描述摩擦边界条件,详细推导了金属二维体积成形过程无网格RKPM法数值模拟的刚度方程,实现了金属二维体积成形过程的无网格RKPM方法数值模拟。(2)对实现二维刚(粘)塑无网格RKPM法数值模拟通用化和自动化的一系列关键技术进行了研究,并提出了有效的算法和处理方法。积分过程利用四边形背景网格;采用矩形影响域和张量积核函数;利用直接迭代法建立迭代过程初始速度场来求解非线性方程组;由直线和圆弧的组合来描述二维任意模具几何形状,将动态接触边界自动识别技术引入无网格RKPM法。提出在刚塑性无网格法中采用修正的罚函数方法来处理体积不可压缩条件和体积闭锁,并从理论上推导证明了修正的罚函数法与体积应变率映射法向常数场进行映射时二者是一致的,修正的罚函数法可以有效解决包括最严重过度约束情形在内的体积闭锁问题以及避免压力震荡现象的产生。提出了一种基于背景积分网格的自适应动态影响域大小确定方法。其优点是:影响域大小随节点疏密程度自适应确定;对线性基,影响域内参与拟合计算的节点数目大约在4~8个,保证了在影响域内包含合适数量的节点和计算精度;基于局部节点搜索而不必搜索全局,提高了计算效率。(3)在二维问题的基础上,将三维刚(粘)塑性流动理论与无网格RKPM法相结合,提出了三维刚(粘)塑性无网格RKPM法。对实现三维刚(粘)塑无网格RKPM法数值模拟通用化、自动化的关键技术进行了研究:积分过程利用六面体背景网格;采用砖形影响域和三维张量积核函数;利用本文提出的基于背景积分网格的自适应动态影响域大小确定方法,对于三维线性基函数取砖形影响域内参与拟合计算的节点数目大约在8~12个,可以达到较好的计算精度和效率;采用网格法来描述任意形状的三维模具型腔曲面,给出了三维无网格法数值模拟系统中动态接触边界自动处理的非线性接触算法,实现了非稳态任意模具形状三维体积成形过程的无网格法分析。(4)自主开发了金属二维、三维体积成形过程无网格数值模拟分析系统-2D-MForming和3D-MForming。分别对二维镦粗、金属圆环压缩和挤压等二维典型工艺以及镦粗、锻造、反向挤压、复合挤压等三维典型工艺过程进行了分析,通过与刚塑性有限元计算结果和实验数据的对比,验证了本文建立的理论与方法的正确性。(5)无网格法的计算效率低下已经成为其推广到实际工程计算中的瓶颈,在三维问题中尤其明显。本文分析了影响无网格法计算效率的主要因素并提出了刚塑性无网格法里系数矩阵半带宽的计算公式,提出并实现了在三维刚塑性无网格法里采用半带宽方法进行分析。通过调整影响域内节点数目大小和适当编排影响域内节点编号(即进行无网格半带宽优化),使得三维无网格法的计算效率得到较大幅度地提高。三维实例分析表明:在相同计算条件和采用相同半带宽的情况下,对于节点数较少(节点数少于1000)的简单三维问题, FEM的计算效率明显高于无网格;对于节点数较多的复杂三维工程计算问题,无网格法与有限元法的计算消耗时间比例可以达到(1.5~2.0):1。