棒材通过锥形模静液挤压成形的无网格法分析.pdfVIP

维普资讯 第21卷第6期 计 算 力 学 学 报 Vo1．21．No．6 2004年 12月 ChineseJournalofComputationalM echanics December2004 文章编号：1007—4708(2004)06—0647—06 棒材通过锥形模静液挤压成形的无网格法分析 卿启湘 ， 李光耀 ， 王永和 ， 刘潭玉 (1．湖南大学 机械与汽车工程学院，湖南 长沙 410082；2．中南大学 土木学院，湖南 长沙 410083) 摘 要 ：建立 了一种基于初始构形及有限变形的粘塑性弹性本构关系，并由空间描述 的Galerkin能量弱交分原 理，经一致转换得一整体拉格 朗 日方程描述下的动量平衡方程．同时经线性化处理给 出了显示中心差分法求解 格式，以棒材通过锥形模 的静液挤压成形为例进行 了全面的EFG法数值模拟，从而证 明了有限变形粘塑性EFG 法对实际成形工艺分析、优化及设计的有效性。 关键词 ：无网格法；有限变形；粘塑性；静液挤压 中图分类号 ：TG 335 文献标识码 ：A 1 引 言 2 基本原理 计算力学的发展面f临着许多难 以处理的问题， 2．1 EFG方法[I] 如金属压力加工塑性成形问题、裂纹扩展问题等。 考虑一个连续体，它初始构形为区域。，即时 传统的数值方法如差分法、有 限元法、离散单元法 构形为区域 ‘，它们之间的关系： 等都难 以处理此类问题，其主要原因是网格的存在 X— X( ，f)一 x + U (1) 妨碍了处理与原始网格线不致的不连续和大变形。 这些基于网格的方法 ，在处理不连续和大变形时， 式中 X为即时构形 ‘的场点位置矢量； 为初始 常用的是网格重构，这样不仅计算费用昂贵，而且 构形。 的场点位置矢量 tu为位移矢量。 会使计算精度严重受损。而无网格方法，网格可以 在 EFG法中，域。 内场点 的参数值 (f)可 彻底或部分地消除，从而保证了计算精度。 以由滑动最小二乘法(MLS)拟合出函数U(，f) 构造无网格方法途径很多，目前最常用是采用 移动最小二乘法进行近似。Belytschko等人提出了 Uh(，f)一∑P()口(，f)一pW()口(，f) i= 1 无 网格 的 Galerkin法 (ElementFreeGralerkin (2a) Method)，简称为EFG 法，它具有计算协调性好、 式 中 P ()为完备多项式基，对于二维问题 ，其线 稳定性强 的特点。在本文中，首先给出EFG法， 性基为 Galerkin能量弱变分原理、接触边界处理，然后给 出一种基于初始构形及有限变形粘塑性本构关系 P ()一 [1 X y]，m一3 (2b) 以及整体的拉格朗 日EFG法离散格式 ，并 以显示 平方基为 中心差分法求解 ，最后以棒材通过锥形模静液挤压 P ()一 [1 x y x xy y]，m一6 成形为例阐明有限变形粘塑性EFG法在指导金属 (2c) 塑性成形工艺上的有效性。