金属体积成形过程的无网格RKPM分析.pdfVIP

第二届中国CAE工程分析技术年会 金属体积成形过程的无网格RKPM分析 刘永辉陈军 (L海交通大学模具CAD国家工程研究中心上海20t)030) 摘要：将刚(粘)塑性流动理论与RKPM相结合．对金属体积成形过程无弼格RKPM数值模拟关键技术进行 了研究．推导了金属体积成形过程无网格RKPM数值模拟的刚度方程．利用RKPM分析了二维和三维金属体 积成形过程，并将数值计算结果与刚塑性有限元及实验数据进行了比较，RKPM与有限元和实验结果均符合 良好，表明了本文方法的正确性．同传统的刚塑性有限元法相比，无网格法由于在构造场函数时只需要节点 信息而不需要划分网格，消除了网格畸变，因而在金属体积成彤领蛾具有广阔的发展与应用前景。 关键词：无网格；再生核质点法；金属体积成形；刚(粘)塑性；不可压缩性；修正曲罚函数浩 中圈分类号1G302 文t标讽碍：A 1引言 金属体积成形技术在金属零件的制造过程中起着十分重要的作用。应用传统刚塑性有限元 法模拟金属体积成形过程时，会在变形剧烈区域产生网格畸变，导致计算停止。无网格法是近 年来发展起来的一种新的数值模拟方法，它将连续体离散为有限数目的节点，在构造场函数时 只需节点信息而不需要节点的连接信息，避免了有限元中繁琐的单元阐格生成，消除了网格畸 变。有代表性的无网格法有光滑质点流体动力学方法SPH、扩散单元法DEM、无网格伽辽金 以SPH、EFGM和RKPM在金属体积成彤中的研究与应用最为广泛。 Chen等I”首先将无网格法用于解决金属塑性成形问题，对弹塑性材料的圆环压缩、轴对 称镦粗、二维板料成形等工艺过程进行了拉格朗ElRKPM数值模拟。Xiong等[L31采用RKPM 对刚塑性微压缩材料的轧制、二维镦粗和反向挤压等工艺进行了数值分析。在国内，李光耀等 l|1对三维弹塑性金属成形问题进行了RKPM数值模拟；崔青玲等91对矩形坯料平面无摩擦镦粗 过程进行了RKPM分析，工件假定为刚塑性可压缩材料。 在传统的金属体积成形工艺(如镦粗、拉拔、挤压等)的数值模拟分析时，由于塑性变形 远人于弹性变形，弹性变形可忽略，往往将材料模型简化为刚(粘)塑性不可压缩材料模型161。 本文将刚(粘)塑性流动理论与RKPM相结合，对金属体积成形过程无网格RKPM数值模拟 关键技术进行了研究，推导了金属体积成形过程无网格RKPM数值模拟的刚度方程。最后利 用RKPM分析了二维和三维金属体积成形过程．数值计算结果表明了本文方法的正确性。 2无网格RK蹦基本原理 Liu等171在积分变换的基础上对核函数进行修正，提出再生核质点法，即 “”(x)=』c(弼X--S)中。(X--S)“o)ds (I) n 式中：UR(x)为速度场函数Ⅳ(x)的再生函数{o。O一5)为核函数；C(罨z—J)为修正函数： s为积分哑元。 修正函数C(x；x—s1定义如F c(x；x—S)=H7(O)M一1(工)H0一s) 式中：H(x-s)为单项式基函数矢量，即 H 7(工一J)=[1，‘一Ⅳl，x2一s2，x3一墨，(_一s1)2，L，(x—J)“】 (3) M(x)=￡H(x～s)H7(x-s)O。(x-s)ds (4) 由于方程(1)完全复制N次单项式，因此这种方法满足N次一致性条件。式(1)的离散形式 为 — 1| S mX S △ Ⅳ “ ) 口 ( ) 咋 5 ilh(z)=∑C(x；x 卉 M∑H 点，相关的值。与有限元法不同，无网格法的速度场函数一般不具有插值性质，却 妒，(x．)≠巧“；需要采用相应的处理技术来施加本质边界条件。本文采用边界奇异核方法吲能 够在无网格法里实现本质边界条件的直接、精确旌加。 3 lll【蹦用于刚(粘)塑性不可压缩性材料 3．1等效、体积应变速率矩