环件轧制三维有限元模拟中质量缩放方法的运用.PDFVIP

环件轧制三维有限元模拟中质量缩放方法的运用 作者：钱东升华林左治江袁银良 摘要：阐述了环件轧制三维有限 元模拟中运用质量缩放方法的意 义及其理论依据，并对显式动力 学有限元模拟中各种质量缩放方 法进行了分类，给出了环件轧制 模拟中质量缩放方法的选取原 则。最后基于 ABAQUS/EXPLICIT 操作平台， 建立了铅环件的环件轧制三维有 例 限元模型，通过具体模拟及分析比较计算结果对选取原则进行了验证，从而总结出了适合于 环件轧制模拟的有效的质量缩放方法。 关键词：环件轧制 有限元模拟质量缩放 案 中图分类号：TG E A 1 前言 C 大变形动态非线性问题（nonlinear problem of dynamic large deformation ）[1]的计算和分 析可以说代表目前有限元分析的最高水平，也是目前还在处于迅速发展的领域,其包括的方 模 面主要有：板料成型工艺的数值模拟（汽车覆盖件），汽车安全性碰撞的全过程计算，材料 的体积成形工艺（挤压、锻造、轧制）的全过程数值模拟等，环件轧制的三维有限元模拟就 属于其中的一种材料体积成形工艺的全过程数值模拟。对于此类问题中存在的几何非线性和 凯 材料非线性问题，用于计算求解的方法主要有两种：静力隐式算法和动力显式算法。静力隐 O 式算法的特征是迭代计算，在每一时间步内需反复迭代，迭代收敛性会受许多因素的影响， E 在计算中须调整迭代以满足收敛，而且其计算时间随模型单元数量呈指数增长，所以其计算 M 时间长。而动力显式算法是递推计算，无需迭代，就时间步长直接进行递推计算，计算时间 随单元数量呈线性增长，计算时间短，而且可以通过质量缩放来缩短计算时间，所以对于大 A 规模计算和高度非线性问题[2]，动力显式算法总体效率要比静力隐式算法高，优势明显。 C 近年来，随着在板料成型及碰撞冲击领域得到日益广泛的应用，对于材料的体积成形工艺领 域，人们也都开始采用动力显式有限元法对其进行过程模拟研究[3]，而且取得了较好的成 效。 动力显式算法在有限元分析中通常被用来解决两类典型的问题：结构动力学响应分析和准静 态分析（其中包括复杂的非线性效应和接触条件），而环件轧制的模拟则是其中的一种准静 态分析。由于动力显式算法采用的中心差分法是通过对时间求导进行计算，引入了质量矩阵 [4]，所以在准静态分析中用于平衡方程式中的离散质量矩阵对计算效率和精度都起着至关 重要的作用。从而在这类问题的模拟中，如何正确适当的运用质量缩放方法，使其能够起到 在保持计算精度和稳定结果的同时提高计算效率的作用，也成为了一个至关重要的问题。目 前，针对这个问题，人们也开始对其进行了相应的分析研究，东北大学在利用显式动力学有 限元法模拟平板轧制过程中，通过计算实例研究了采用提高模拟轧制速度和质量缩放技术来 缩短计算时间等方法的有效性，在忽略惯性惯性力影响的前提下得出了与实验结果相符合的 计算结果[5]。燕山大学在用显式动力学有限元法分析板带轧制压力分布时提出了为了加快 计算速度而采用的虚拟速度和虚拟质量带来的虚拟惯性效应会影响计算的结果[6]。而在环 件轧制成形模拟方面，并没有人通过具体的实例分析来验证质量缩放技术在环轧模拟中运用 的有效性和合理性。本文运用Abaqus 软件建立环件轧制的三维有限元模型，并采用动力 显式算法对铅环件的轧制过程模拟进行了实例计算，验证了质量缩放技术在环轧模拟中运用 的有效性，并通过对计算结果的分析比较，给出了适合环件轧制模拟的有效的质量缩放方 法。 2 质量缩放在环轧模拟中的意义及理论依据 在Abaqus 中利用Explict[7]来进行准静态模拟分析的求解，为了在较快的时间内获得可以 库 接受的静态解