汽车覆盖件工具曲面有限元网格自适应转换方法.pdfVIP

第23卷第3期 计算力学学报 V01．23。NO．3 of ChineseJournal Mechanics June2006 2006年6月 Computational 汽车覆盖件工具曲面有限元网格自适应转换方法 龚科家， 杨 光， 胡 平’ (吉林大学车身与模具工程研究所汽车动态模拟国家重点实验室，长春130025) 摘要：提出一种适用于汽车覆盖件曲面有限元网格转化和在单元水平上提高模拟精度的方法，将平面下通过合 并三角形单元成四边形单元的有限元网格转换方法的应用范围扩展到曲面，并且降低了对初始网格形状的要求。 算法的关键在于增加了对曲面相邻单元不同夹角情况下的处理和优化规则，以使其能够更好地拟合原始CAD 曲面。 关键词：工具曲面；三角形单元；四边形单元；有限元网格 中图分类号：TP391 文献标识码：A 自适应转换方法可以为已有三角形网格的汽车覆 1 引 言 盖件曲面生成四边形网格提供一条有效途径。另 在板金属成形领域，由于所涉及的数学和力学 外，在CAE分析时应用该算法可以由三角形网格 问题的复杂性，弹塑性大变形有限元数值仿真方法 快速生成四边形网格，避免了返回前处理重新划分 成为该领域目前广泛采用的分析手段。众所周知： 四边形网格的冗杂步骤。 从有限元单元理论本身而言，四边形单元的模拟精 本文针对已有的平面上三角形网格转换的方 度通常要比三角形单元高；而另一方面，将十分复 法[8以嵋推广到三维曲面网格生成时效果并不是十 杂的汽车覆盖件三维曲面转化为可以进行有限元 分理想的情况，提出了一种适用于曲面的有限元网 分析的网格曲面，通常转化为三角形单元比较简 格自适应转换算法，在平面算法用口和．8值来衡量 便，相对容易实现。目前，国际上许多一步逆成 三角形和四边形质量的基础上，通过增加对曲面相 形[13有限元商业化软件，如AutoForm／One-step， 邻单元夹角的角度量因子0的处理规则，以使其能 对工具曲面均采用简便的网格转换方法。如何将 够更好地拟合原始曲面。该算法对曲面的初始网格 两个方面的优势发挥出来，即采用三角形网格离散 要求并不严格，可以是全三角形网格或者三角形与 工具曲面，而采用四边形网格进行有限元计算，是 四边形的混合网格，同时采用分级合并实现了与平 本文工作的主要目的。 面三角形转换算法相同的O(N)时间复杂度，显著 有限元网格划分在有限元分析中的作用非常 提高了曲面网格的质量。文中着重阐述了角度量因 关键，也是有限元前处理的重点和难点。与成熟、 子0的求法、数据结构的建立和分级合并的处理方 简单的全自动三角形网格生成算法[2‘4]相比，目前 法等关键技术，最后用Laplace平均法n3结合基于 只有少数方法能生成曲面四边形网格，如铺路法、 顶点和对角线优化的方法[1妇进行网格光顺，改善 映射法和模板法等[4。]，而且这些算法往往比较复 网格质量。 杂或受到一定的限制。为了减少四边形网格生成 的难度，通常采用间接算法。对于许多复杂的汽车 2 平面上网格转换方法概述 覆盖件来说，由于大量覆盖件比较容易得到或者已 2．1 三角形的口值 经有三角形网格，本文提出的曲面上三角形网格的 三角形ABC的质量可以用式(1)的口值来衡