应用计算机图形技术划分单元网格.pdfVIP

应用讨‘算机图形技术划分单元网格 王伟 游亚戈 10070) (中国科学院广州能源研究所，广州5 摘要 本文研究了运用计算机图像，图形技术划分计算对象单元网格的方法。猷应用边界元法数 值模拟某一波浪电站的波浪作用过程为例，探讨了利用图像扫描／处理技术与CAD控术，直接生成描 述不规则地形的几何模型、自动划分单元网格以及校验单元划分正确性的方法。实验表明，直接运用 成熟的计算机图像，图形处理工具，可以准确描述计算对象的几何结构，自动划分具有任意不规则形状 的对象以及对所划分单元进行有效的检验。与常规方法相比，这种方法具有划分精度高、工作量小的 特点。 关键词 图像，图形处理，CAD技术，单元弼格划分，计算流体力学 1前言 对几何体的单元网格划分是流体力学数值计算的关键步骤，是数值分析工作者面 临的一个共同问题。不同的计算方法所采用的单元网格划分方法各有不同，但保证划分 单元准确地描述计算对象的几何结构则是一个基本要求。对具有简单边界形状的几何体 而言，这一要求可通过适度编程来实现。但计算对象通常是一个形状复杂的三维几何体， 要在几何上作准确的单元划分确非易事。 一般而言，对一个三维几何体的单元网格划分有两利，方式：其一是三维实体划分， 采用四面体或长方体单元来描述整个几何体(如有限元法、有限差分法等)；其二是三 维边界划分，采用空间四边形或三角形单元来描述几坷体的内外边界曲而(如边界元 法)。三维几何体形状越不规则，单元网格划分的工作量就越大。通常在单元网格划分 时，针对每一个具体计算对象都需用手工准备庞大的几何描述数据文件，再编制专门的 划分程序来划分单元网格。而手工准备数据文件是一件十分艰巨的工作，耗时费力，而 且极易出错；划分单元所采用的划分程序又往往是针对某～具体计算对象编制的，很难 具有通用性。至于所划分的单元网格的正确性校验，目前也没有较成熟、较简便的方法， 而确保划分单元网格的正确性(几何上不重复、无遗漏、单元外法向一致)是整个计算 模拟成败的关键。因此，单元网格的划分工作十分重要，其工作量在整个流体力学计算 工作中占了很大的比重。 现代计算机技术发展十分迅速，硬件水平和软件技术可以蜕是日新月异，特别是 图像／图形处理技术发展更快，对一任意形状的三维几何体作精确描述已十分容易。以 流行的大众化矢量图形处理平台AutoCAD为例，其数据文件中包含了足以精确描述任 意形状三维几何体结构的几何数据，并且可以很方便地从文件中提取所关心n勺数据。这 样就使我们有可能直接运用这些优秀的通用图形处理软件来解决单元网格划分中的难 题，而避免或减少编制仅仅适合于某些特定形状几何体划分的专用划分程序，位数值分 析工作者从繁重的数据准备工作中解放出来，达到事半功倍的效果。 本文以某～边界元计算为例，探讨了运用图像扫描／处理技术与CAD技术实现计 ·】8· 算对象边界的单元网格划分以及校验单元划分正确性的方法。该算例的计算对象为一个 具有复杂地形的岸式波浪电站，对象的原始几何资料为工程用手工描绘的等高线图。从 此图出发，借助于AutoCAD平台，运用某些特殊的图形处理工具，实现了复杂形状对 象的边界单元网格的划分，以及所划分单元的正确性校验。与传统方法(人工描述加程 序处理)相比，该方法不需手工准备几何描述数据文件，不需复杂编程，直观可视，失 误率小，划分精度和划分效率都有很大的提高。 2基本方法 运用边界元法研究波浪的作用过程与地形影响，要求将选定计算域的水体(计算 对象)的外边界作空间四边形或三角形单元网格划分。因此，首先需要将计算对绿输入 计算机，以供进一步处理。 图1计算对象的原始资料 图l所示为某一岸式波浪电站的原始几何资料。图中数据为水位值，曲线段为实 测等高线。 通常做法以手工方法为主，先用控制点与折线段摘述等高线，节点输入计算机后， 再用手工或程序形成所需要的网格单元。 本文所采用的方法以计算机图像／图形处理技术为主，其主要步骤为图纸扫描与识 别、网格单元划分与输出以及划分单元的检验与改正等。 3计算对象的几何描述 3．1图像的预修整