基于背景网格变形的动态网格移动方法 Fast dynamic mesh moving based on background grid morphing.pdf

第29卷第1期 计 算 力 学 学 报 VoI．29，No．1 2012年 2月 ChineseJournalofComputationalMechanics February2012 文章编号 ：1OO7—47O8(2O12)O1—0105—06 基于背景网格变形的动态网格移动方法 林天军， 关振群 (大连理工大学 工程力学系 工业装备结构分析国家重点实验室，大连 116024) 摘 要：基于Delaunay背景网格插值技术的动态网格生成方法无需迭代计算，效率较高。但对复杂构形大幅运动 的动边界 问题 ，尤其当边界大幅转动时，背景网格极易交叉重叠 。重新生成背景网格和重新定位 网格节点信息不 仅费时而且会导致网格质量的严重下降。本 文提出改进的基于背景 网格 的动态网格 变形方法，通过在初始 Delaunay背景网格中添加辅助点，生成一层新的背景网格和新的映射关系；采用ball-vertex弹簧法驱动新背景网 格的变形，进而牵动 目标 网格的变形。算例表明，本文提出的动态网格变形方法对所关心区域的网格具有 良好保 形性，边界可转动更大角度而不会 出现网格交叉重叠问题，总体上提高了动态网格更新的效率和质量。 关键词：动态网格；网格变形；Delaunay背景网格；运动边界；ball-vertex弹簧法 中图分类号：0242．21 文献标志码 ：A 关，从而防止单元各边重叠甚至交错导致网格变形 1 引 言 失败。但由于拉伸弹簧与扭转弹簧针对 的变形状 流固耦合现象越来越被关注 ，可通过数值模拟 况不同，应用于大扭转变形情况时结果仍不够理 研究流固耦合发生时流场与结构的动力学现象。 想。 数值模拟涉及到计算流体力学、计算结构力学和动 最近 ，刘学强等口提出基于 Delaunay背景网 态网格计算三个研究领域。随时间变化的网格称 格插值技术的快速动 网格变形方法 ，该方法首先用 为动态网格。动态网格的适应性调整方法有两种 ： 计算域的边界点和 Delaunay三角化方法[8 生成 一 是随时问步推进 网格重新生成[1]，该方法计算耗 背景网格--Delaunay背景 图，再将计算 网格点逐 费很大，且会 引入较大误差；二是 网格的移动变形 个映射到背景图单元 。背景 图随着边界移动 ，计算 方法 ，即初始网格随时间步推进而变形 ，以保持边 网格根据映射关系随着背景图的移动而变形 ，网格 界的一致性 。网格变形方法主要有弹簧近似方法 变形过程无需迭代 ，效率较高和稳定性好 ，随后得 和弹性体方法。弹性体方法的网格变形能力强且 到应用[1。。，但是对较复杂的有多体相对运动或大 网格质量较高，但计算工作量大，弹簧近似方法则 变形 ，尤其是有大幅相对转动的问题，背景 图单元 较为简单高效 。弹簧法的概念 由Batina~23提出，即 极易出现交叉重叠。文献ET]虽然给出了背景图出 把流场网格各节点之间的连线看作弹簧，弹簧刚度 现交叉时的处理策 略，即重新生成 Delaunay图和 与节点间的距离成反 比。弹簧模拟在气动弹性问 重新定位网格点的信息，但这会导致新的问题：(1) 题 中应用较多[3“]，可以处理边界小变形的情况 ， 重新生成背景图和定位网格点要花费更多的计算 但在某些复杂几何变形[5或大变形情况下会失效 。 时间；(2)背景图交叉后 ，重新生成背景图和重新定 Farhat等 认为拉伸弹簧模拟只假设了弹簧的轴 位网格点，相当于丢弃初始 网格质量的良好特性 ， 向伸缩 ，在变形较大情况下，不能阻止弹簧