板材冲压成形性模拟CAE 软件系统KMAS.pdf.pdfVIP

板材冲压成形性模拟CAE 软件系统KMAS1 胡 平 柳玉起 王锦程 李运兴 那景新 杨 光 吴志强 郭 威 张向奎 (吉林工业大学 汽车覆盖件成形技术研究所，长春, 130025) 摘要:本文阐述了自行研制的覆盖件成形性仿真CAE 软件系统KMAS 的基础理论方法和关键技术。对一些 典型冲压过程的板材成形缺陷进行模拟预示，并与实验结果相比较，旨在检验 KMAS 软件在金属塑性成 形屈服准则、塑性流动本构理论、以及弹塑性大变形接触问题有限元方法的有效性，为解决金属覆盖件 塑性加工成形与模具制造工业面临的关键力学问题提供有效的数学工具。 关键词 覆盖件成形，板材成形缺陷，CAE 软件系统KMAS 1 前言 板材冲压成形技术已广泛应用于汽车、 机电等机械制造行业， 在塑性加工领域中占有重要的 地位。由于对冲压成形产品的形状、质量及精度的要求日益提高，促使工程技术人员不断地探索克 服板材冲压成形缺陷的有效途径。冲压成形过程中经常发生的起皱、鼓动、制耳、破裂以及成形后 的回弹等主要冲压成形技术缺陷一直是实际工程技术人员十分关心的力学疑难问题，同时也是制约 我国覆盖件冲压模具制造业发展的关键所在。 针对汽车车身及模具制造业目前的状况，以本文作者为骨干的项目组经过十几年的研究工作积 累，独立开发了专门用于复杂冲压件成形性模拟仿真的 CAE 软件系统 KMAS (King Meshing Analysis System) 。本文着重介绍了以笔者提出的塑性本构理论和基于率形式的Updating Lagrangian 弹塑性大变形接触问题有限元方法为基础的KMAS 软件的关键技术以及结构框架；作为理论与方法 可靠性的验证，文中例举了三个具有代表性的成形实例，并与实验结果进行了比较。 2 板材成形弹塑性大变形有限元仿真理论 2.1 有限元算法 目前在国际上常见的弹塑性大变形有限元算法分为两大类。一类为隐式的UL 时间积分算法 [1] （implicit time integration scheme ） 。这类算法的显著优点是能够连续地模拟板料从塑性加载直至 弹性卸载全过程，并允许采用较大的加载时间步长。但由于遇到复杂的接触问题往往迭代计算过程 难以收敛以及求解切线刚度矩阵需要耗费较大的计算机存储和计算时间等原因，另一类显式时间积 [2,3] 分算法（explicit time integration scheme ） 被普遍用于目前的商品化软件中。这类算法采用动力学 的中心差分格式，显著优点是对变形相关接触问题算法稳定，且不需直接求解切线刚度。由于算法 本身是条件稳定的，加载时间步长必须充分小，计算时间直接受到最小单元尺寸的限制。本文采用 率形式的 Updated Lagrangian 有限元列式方法。在工具和成形板料的接触界面采用了修正的 Coulomb 定律。以t 时刻构形为参考构形的虚功率增率原理为[4,5] ∇ S ) ∫V ( i j −σ ε −σ ε +σ v δ v d V i k k j k j k i i k j ,k i , j （1） ∫V P δ v d V + ∫A P δ v d A + ∫A f δ v d V i i i i i i