金属塑性成形过程无网格伽辽金方法及其关键技术研究-材料加工工程专业论文.docxVIP

山东大学博士学位论文摘要 山东大学博士学位论文 摘要 金属塑性成形技术在金属零件的制造过程中起着十分重要的作用。它不仅具 有生产效率高、产品质量稳定、原材科消耗少的优点，而且还可以有效地改善工 件的组织性能。随着计算机技术的发展和数值计算方法的同益完善，有限元方法 在工程实际中得到了广泛的应用。但是，当工件变形到一定程度时，有限元网格 将产生畸变现象，此时，这种以单元作为基本概念的有限元方法面临着一些难以 处理的问题。无网格方法基于离散节点的近似，避免了有限元方法对于网格的依 赖，在涉及到网格畸变的大变形问题分析中具有一定的优势，并且在数据准备和 后处理方面也比有限元方法灵活简单。 无网格方法作为一种较为新颖的数值方法，经过十余年的发展，已经逐渐应 用于金属塑性成形过程的模拟，并且取得了一定的成果。但是对于非稳态大变形 金属塑性成形问题，由于成形过程的复杂性，许多关键应用技术还有待于研究。 尤其当金属流动比较剧烈或变形比较复杂时由于节点分布的非均匀性增强和变形 域形状更加复杂，速度场的近似精度和变形域的积分精度会有所下降，进而导致 数值模拟出现错误。因此，本文针对非稳态大变形成形过程，主要研究了基于刚 (粘)塑性材料假设的无网格伽辽金方法、关键处理技术及其在金属塑性成形过 程分析中的应用。 基于无网格近似方案，结合伽辽金离散方法，提出了刚(粘)塑性无网格伽 辽金方法。速度场采用经过变换法修正处理的最小二乘近似和再生核质点近似或 者直接采用径向基函数插值进行近似。采用反正切摩擦模型描述摩擦接触边界条 件。对于模具形状任意的塑性成形过程，在局部坐标系下施加摩擦力边界条件， 给出了局部坐标系和整体坐标系的变换矩阵，解决了模具形状任意的二维塑性成 形过程分析中的摩擦接触边界条件施加问题。在此基础上推导建立了金属塑性成 形过程无网格方法分析刚度矩阵方程。进而采用直接迭代法获得初始速度场，利 用Newton。Raphson迭代方法求解刚度方程。最后给出了模拟等温塑性成形问题的 分析步骤。 针对非稳态金属塑性成形过程的复杂性、数学处理上的困难性和在成形过程 模拟分析模型建立以及分析程序的通用化方面存在的不足，研究了刚(粘)塑性 无网格伽辽金方法应用于金属塑性成形过程分析的关键处理技术。对于任意边界 形状的二维金属塑性成形过程的无网格伽辽金方法分析，建立了诸如任意形状模 具描述方法、迭代收敛判据、接触脱离判断等问题的处理方法，实现了非稳态任 摘要意模具形状塑性成形过程的无网格伽辽金方法分析，提高了分析程序的通用化程 摘要 意模具形状塑性成形过程的无网格伽辽金方法分析，提高了分析程序的通用化程 度。重点研究了应用无网格方法分析金属塑性成形问题时所面临的体积闭锁问题 的缓解算法，采用体积应变率映射方法，对能量速率泛函中的体积应变率进行修 正处理，通过将速度场计算的体积应变率映射到低维空间，降低独立约束方程数 目，从而建立了体积闭锁现象的缓解算法。 当金属流动比较剧烈或变形比较复杂时，近似精度和积分精度会有所下降， 为了保证分析精度和提高效率需要进行一定的人为的干涉，所以自动化程度比较 低。为此，本文提出了自适应刚(粘)塑性无网格伽辽金方法。无网格伽辽金方 法基于紧支试函数加权残量法，所以权函数影响域的确定会影响近似的精度，本 文提出了修正的权函数影响域的确定方法以提高近似的精度。无网格伽辽金方法 的分析精度依赖于变形域的积分精度。积分时，背景网格的疏密和网格内的高斯 积分规则对计算精度有一定的影响，尤其在轴对称非稳态金属塑性成形过程模拟 中，如果以固定的背景网格来获得高斯点时，可能会造成高斯点减少，导致精度 损失。文中提出了自适应的积分背景网格的动态划分方法，通过控制高斯积分节 点数目，使其保持在一定高斯积分节点数量，以保证积分的精度。在大变形的成 形过程分析中，一部分摩擦接触边界上的边界节点密度会随着变形的累积而逐渐 变小，从而边界条件约束能力下降，导致误差的产生，累积到一定程度会出现内 部节点穿透边界的现象。因此通过动态控制边界节点分布密度来增强边界条件的 约束能力，以保证分析的精度。 对于中高温条件下的塑性成形过程，推导了刚(粘)塑性无网格伽辽金方法 热力耦合分析公式，给出了热力耦合分析步骤，拓展了无网格方法在塑性成形过 程分析中的应用范围。在热力耦合分析时，变形初始阶段有时由于热载荷变化速 度过快而产生温度场震荡的现象，文中采用集中热容法修正热容矩阵解决了温度 场震荡的问题。此外，对于Mn含量低于1％的碳锰低碳钢，采用Y耐a模型描述 微观组织演变过程，实现了动态再结晶发展演变过程的无网格方法模拟。完善了 无网格伽辽金方法在塑性成形过程无网格伽辽金方法分析中的应用。 基于上述剐(粘)塑性无网格伽辽金方法理论研究与关键技术的开发，采用 c++