三维刚塑性有限元中网格重分技术的研究.pdfVIP

维普资讯 第21卷 第3期 太 原 重 型 机 械 学 院 学 报 Vd．21 №3 2000年 9月 JOURNALOFTAIY~ANHEAVYMACHINERYINSTITUTE Sep．2000 文章编号：1000—159X(2000)03 0221 04 ≥＼一 ．≥委 三维刚塑性有限元中网格 技术的研究 周 一 刘建生，郭晓霞 院’太原 。。。 1 Dl 摘 要：本文在深入研装究了5-维网格重分技术的基础上，针对轴娄锻件的网格重分提 n 生 出了一种简洁、有效的算法——投影纽合法，并用等参逆变换法对重分前后的场量进行转一r ut- 换。模拟实倒袁明，该算法生成的网格质量好、数据转换精度高。 l 宦 降 ．1耽 一 一 一 0 引言 1 网格重分 j蟊7G 在金属成形的三维有限元数值模拟过程 中，制 投影组合法的基本思路如下：根据锻件的外部 约着三维有限元应用的一个关键 问题是网格重分 形状进行分段。对每段向平行于轴径的平面投影 ， 技术。塑性有限元法模拟金属大变形时，由于摩擦 投影区域经过布尔运算后形成独立的子区域。对 的作用而产生变形的不均匀性，随着变形的增加容 子区域进行二维网格划分，对不同的段调用相应的 易导致网格畸变而逐渐失去其描述变形规律的精 投影区域，形成该段投影区域的单元拓扑关系。沿 度。当畸变单元达到一定数 目时，就需对网格重 纵向剖切组合成六面体单元 此算法中的主要步 分。目前 ，网格重分的算法很多，如映射法、填充 骤如下 ： 法、拓扑分解法、八叉树法等等 这些方法在二维 1．1 形成二维子区域 问题中已经 比较成熟，但在形状较复杂的三维体积 成形方面仍然面临较大的难度。三维方面比较成 不同的段向垂直于轴径的平面投影，各 自的投 熟的TK．Ham…的分区划分法、D Y．Yang[刘的模 影区域将出现包含、相交、相离等不同情况。此时 块法、Delaunay[三角化方法等等。这些方法划分 需对投影区域进行布尔运算，形成待划分的二维子 四面体单元比较成熟，但在划分六面体单元时各有 区域 如图1所示，A、B分别表示不同段的投影区 缺陷，仅在一定范围内适用。本文在综合分析现有 域，进行布尔操作后分成三个子区域。 的三维网格重分技术的基础上，提出了一种简洁、 1．2 二维网格划分 有效的算法．投影组合法，并采用等参逆变换法进行 场量的转换，开发了三维网格重分系统。 采用层进法 【一对二维区域划分，此算法对边界 离散后．自动生成优化后的网格。离散过程 中可根 收稿 日期 ：1999—03—11 作者简介：周晓章(1964 )，男，太原重型机械学院科技外事处工程师，研究方向为有限元数值模拟研究。