非结构网格自适应细化的实现与应用.pdfVIP

学兔兔 第 1期(总第 176期) 机 械 工 程 与 自动 化 No．1 2013年 2月 MECHANICAL ENGINEERING AUTOMATION Feb． 文章编号 ：1672—6413(2013)01—0001—03 非结构网格 自适应细化的实现与应用来 李于锋 (中国工程物理研究院计算机应用研究所 ，四川 绵阳 621900) 摘要 ：基于非结构网格 的有限元法在 CAE中有着广泛 的应用。网格 自适应细化方法利用后验误差估计 自动 决定局部细化，用较低 的计算代价获得较好的计算精度。非结构网格的 自适应细化实现相对复杂 ，同时还要 保证 网格的质量 。基于z—Z后验误差估计对三角形和 四面体 网格进行网格 自适应细化，给出了具体 的实现算 法，开发了自适应计算模块，达到了提高计算效率的目的，通过算例表明自适应实现的正确性和可靠性。 关键词 ：白适应网格细化；非结构网格 ；有限元法 中图分类号 ：TP391．7：TB115 文献标识码 ：A 0 引言 e-~-- 一 。 … … … … … … … … … … … … … (1) 在工程数值模拟领域 ，针对复杂结构的有限元计 考察常见 的线性微分方程 ： 算应用十分广泛。在机械零件制造 中如何高效得到结 L “+ 一0 。 …… …… …… ……… ……… … (2) 构关键点的应力，以及在流场中如何快速准确地捕捉 其中：L为线性微分算子，P为已知函数 。 激波等问题 ，若单纯依靠不断全局细化 的网格增加计 定义能量范数为： 算量，是不经济的。网格 自适应技术兼顾了计算精确 lj一 [IeTLedl~]。一[I(“一五)L(u一 性和计算效率，能够使用较少的网格计算代价获得较 为准确的计算结果uJ。 )dn] 。 ……………………………………… (3) 自适应网格技术在 国外发展十分迅速，在成熟的 对于弹性 问题 ，这个能量范数可写为 ： 商业软件中都有功能体现，如LS—Dyna的重划分网 lIglI一[I(＆)DSedg2]。。…………… (4) 格的 自适应实现等 ；在美 国Sandia实验室开发的多物 其 中：D为弹性矩阵；S为应变关于位移的算子。由算 理耦合计算框架 sIERRA2『中，也实现了H 自适应 的 子 S定义的应变和应力分别为 ： 策略。在 国内，北京应用物理与计算数学研究所开发 ￡= Su， —S五 。 … … … … … … … … … … (5) 的JASMIN[3结构化网格 自适应框架已经在多个应用 =Ds， 一D 。 … … … … … … … … … … (6) 程序中得到实际应用。中科院科学与工程计算国家重 将式(1)代入式(4)，联立式(5)、式(6)，能量范数 点实验室开发 的PHG_4能够对 四面体 网格进行 自适 也可以写为： 应细化。然而涵盖二维和三维非结构化 网格 自适应的 软件模块还