一种基于能量理论的板料成形反向模拟初始域计算方法.pdfVIP

一种基于能量理论的板料成形反向模拟初始域计 算方法 (上海交通大学塑性成形工程系上海200030)赵震，唐炳涛，阮雪榆 摘要：针对目前在板料成形反向模拟法中常用到的初始域计算方法的缺点，提出采用一种基于弹簧系统能量理论 的空闻网格映射法。该算法克服了垂直投影法因零投映面积而不能迭代计算的问题，解决了几何映射法需要寻找 对称线、断面线展开法无法解决局部网格畸变的问题。NuⅦsHEET2002一标准考题中，通过与增量模拟软件 Dynaform的比较，证明了采用基于弹簧系统能量理论的调和映射思想来求解反向模拟法初始域具有较好的适用 性，同时验证了反向模拟法较为可靠的计算结果及较高的计算效率。 关键词：板料成形初始域反向模拟调和映射 中图分类号：TG386 量较大。 1前言 因为反向模拟法所要解决的问题之一正是零件 的初始构形．并且由于冲压件形状的复杂性和多样 目前，常用的板料成形数值模拟方法主要有两 性，很难找到一种适用于任何构形的初始域计算方 种：增量方法和反向模拟方法…。板料成形数值模 法。由于几何方法不能精确考虑到材料的性能和变 拟中的增量方法是基于流动理论。采用有限变形弹 形的特点，本文基于大量的分析比较，提出了一种 塑性有限元方法，并采用静态隐式、静态显示或动 基于弹簧系统能量理论的调和映射法来解决板料成 态显示算法求解。基于形变理论的板料成形反向模 形初始域计算问题。实例证明该方法可以很好地解 拟法(也称一步模拟法)具有信息要求比较低、计 决垂直壁问题、多凸起特征及负角零件问题等，具 算速度快等优点。利用反向模拟法软件可以预测成 有良好的通用性。 形过程中可能出现的缺陷(例如破裂和起皱)和毛 坯形状。和优化方法相结合，可以对压边力、工艺 2反向模拟法的基本理论 补充面、拉延筋等进行快速优化设计。 反向模拟法实质上相当于有限元增量法的一个 反向模拟法假设成形过程比例加载，仅仅考虑 加载步，除了不用作板料和模具的接触搜索和判断 初始的毛坯和变形终了时的状态，不考虑变形过程 外，和增量法的实现比较类似。如果采用 中的中间状态(图1)。基本思想是1161：从最终工 Newton-Raphson方法来求解非线性方程组，则箕牧 件C出发，将其作为变彤终了时工件的中面．对其 敛性要求给出的初始域和精确解足够接近才篚保 离散，通过有限元方法确定在满足一定的边界条件 证。垂直映射法较简单，但如果单元的法向和水平 面平行时，其中水平面上的投影面积为零，无法进 下工件中各个节点P在初始毛坯Co上的位置Po。 行迭代。几何映射法口’考虑了由于几何形状特征导 致的变形特点，但它要求能找到工件上的对称线， 即已知条件为最终工件的几何形状和初始状态的厚 而实际上在很多复杂的汽车覆盖件上很难找到对称 度：未知量是初始状态的形状和最终构形的厚度。 线。射线展开法(也称截面线展开法)J。吸收了 基于几何展开方法的优点，同时克服了几何映射法 必须寻找对称线的困难，使得反向模拟法更具有实 用性。但当构件上有比较剧烈的局部胀形或者比较 狭长时，会导致相邻的两条剖面线长度差别很大， 严重时可能使初始域单元重叠。也会使计算失败。 基于有限元模型的几何影射法p1采用由工件中心向 外周逐层映射的思想，可以解决垂妻壁的问题，但 是该方法的计算成功与否与零件的复杂程度密切相 图l反向模拟法示意图 关．在网格畸变较为严重区域难以凑效，并且计算 schematicrelationsofinversemethod