基于无网格Galerkin法模拟金属塑性成形的研究.pdfVIP

基于无网格Galerkin 法模拟金属塑性成形的研究1 刘新，谢桂兰，彭建新 湘潭大学机械工程学院，湖南湘潭（411105) E-mail ：xieguilan@ 摘 要：本文利用无网格Galerkin 法和弹塑性小应变理论，建立了金属塑性成形的EFG 求 解模型，其中采用变刚度法进行增量迭代计算，并利用坐标变换方法使之满足本质边界条件。 编程对压印和斜扎两个塑性成形过程进行了模拟研究，讨论了离散节点数目和迭代步数对计 算精度的影响。模拟结果表明所提出的方法在模拟金属塑性成形问题上是可行的，为基于 EFG 法在大变形问题的研究奠定了基础。 关键词：无网格Galerkin 法，小应变弹塑性理论，金属塑性成形，坐标变换方法，牛顿-拉 夫森迭代法 中图分类号：TG3 1 前言 金属塑性成形过程具有物理非线性和几何非线性等特点，对该过程的模拟，弹塑性有限 元法一直是研究和应用的主流。有限元对网格具有依赖性，容易出现网格畸变现象[1,2,3]，而 无网格法采用基于点的近似，可以彻底或部分地消除网格，不需要网格的初始划分和重新划 分，不仅可以保证计算的精度，也可以减小计算难度，为塑性成形的数值模拟技术带来了新 [4,5] 的希望 。 1994 年，Belytschko 等[6]在DEM 基础上进行了两点改进，提出了EFG(Element-free Galerkin Method) 。2001 年，Belytschko 和李光耀首次将其应用到平面应变镦粗和反挤压的 [7] 模拟中得到了理想的结果 。从此，EFG 法在金属塑性成形模拟中的应用引起了国内外学者 的浓厚兴趣。 目前，随着无网格理论的不断完善，EFG 模拟金属成形取得了很大成果，尤其在体积 成形领域[8,9] 。本文基于EFG 法，首先推导了金属塑性成形过程中的材料本构关系及模拟计 算公式；其次对整个工艺过程，通过施加位移载荷，编写了模拟程序，取得模拟数据；最后 对比有限元软件(ANSYS)的模拟结果，考虑影响结果精度的节点离散数目和迭代步数，得出 相关结论，为无网格模拟金属塑性成形奠定了基础。 2 理论分析 2.1 无网格伽辽金法 EFG 方法采用移动最小二乘法(MLS)来构造近似函数[8] h u (x ) φ(x )u (1) h T 式中u (x )代表近似函u [u u u ] 为节点参数向量。φ(x ) 为形函数矩阵，可表示为 1 2 n φ(x ) [φ(x )=,φ (x )φ (x )] 1 2 n (2) P T (x )A-1 (x )B (x ) T 其中 A P W (x )P (3) T B P W (x ) (4) 1 本课题得到国家自然科学基金 、湖南省优秀青年项目(06B096)、湘潭大学博士启动基金 (06QDZ16)项目的资助。 - 1 -