有限元法处理金属塑性成型过程的接触问题.pdfVIP

第 7 卷 第 4 期 7　 4 塑性工程学报 V o l N o 2000 年 12 月 JOU RN A L O F PLA ST IC IT Y EN G IN EER IN G D ec　 2000 有限元法处理金属塑性成型过程的接触问题 (北京科技大学　100083) 　苏　岚　王先进　唐　荻　孙吉先　田荣彬 摘　要: 用有限元法模拟金属塑性成型过程的难点之一是接触边界条件的处理, 本文着重介绍了处理金属塑性成型 过程中接触问题的一般方法, 包括接触边界的搜索、接触力的计算及接触热传导的计算三方面内容。运用AN SY S 提 供的接触单元, 模拟了Y 型三辊轧制过程及钢轨矫直过程的接触状态, 结果与实际相符。 关键词: 接触边界条件; 有限元法; 金属成型 中图分类号: T G 302 何时发生接触。为了搜索接触边界, 首先要对可能发 1　引　言 生接触的边界进行有限元离散, 即在这些边界之间建 立一层虚拟的接触单元。大多数情况下, 在有限元模 近十几年来, 随着计算机硬件能力及软件水平的 型中不可能处处形成精确的节点对节点的接触, 因此 提高, 有限元法在金属压力加工行业中得到进一步的 接触单元一般表示为点对面或面对面接触。图 1 所表 应用。在对金属塑性成型过程进行有限元模拟时, 如 示的就是一种二维的点对面接触单元示意图, 其几何 何处理工具及工件之间的接触将对模拟结果产生重 形状为三角形, 底边为其中一个接触边界 ( 目标面) 要影响。 ( ) 上单元的边, 顶点为另一个接触边界 接触面 上单 接触问题的处理是用有限元法模拟金属塑性成 元的节点。这种接触单元的每个节点都具有X 和 Y 型过程的一个难点。首先, 在金属塑性成型过程中, ( 两个方向的位移自由度 若要考虑接触热传导, 则还 接触边界状态随变形过程的发展不断变化, 事先难以 ) 需增加一个温度 自由度 。另外值得注意的是, 由于 确定。其次, 接触面之间存在着摩擦和相对滑动, 而 随着变形的发展, 接触面上的每一个节点将会有可能 摩擦的机理仍然是尚未解决的难题, 目前一般采用简 与 目标面上所有的单元发生接触, 所以在离散接触边 化的理想模型来进行模拟计算, 如滑动库仑摩擦模 界时要保证没有遗漏, 否则会导致结果错误。 型、粘着摩擦模型。另外, 接触问题大多属于单侧接 触, 而且是不可逆的, 在数学上可用单边约束不定变 分描述, 但接触和摩擦构成的泛函一般不可微, 有时 非凸, 在数学处理上比较困难。正是由于这些原因, 使得用有限元处理这类问题非常困难。 本文将简要介绍有限元程序中处理接触问题的 一般方法, 着重研究了运用商业有限元分析软件 AN SY S 所提供的接触单元处理 Y 型三辊轧制过程 和