网格剖分之人头.doc

网格剖分之人头——ICEM 先来看看头像的原始几何模型吧，原始文件是影像扫描技术成型的自由曲面蒙皮模型，当然先要把曲面转换成实体格式（这一步对专业的网格平台不是必需），还好，手头先进的三维CAD程序正好发挥用途，一个简单的操作，便完成了曲面向实体的转变，很显然，这是一个对称体，为了减少工作量，当然，对称化的处理是必须的（见图1、2），对不住了，我只有残忍地把这家伙一劈为二了。 图1 整体头像 图2 对称头像好了，完成了基本头像几何的边界，接下来的工作应该是几何拓补规划，把脑筋转起来，OK，很快就发现了，整体头像的几何拓补非常简单，几乎就是一个方块+圆柱+切削圆球的组合体（图3、4），麻烦的制造者在于这家伙脑袋上长了眼睛和一个嘴巴，嘴巴是一个洞穴，眼睛也是一个洞穴，显然，最终拓补构造的时候，这个地方可能得多挖两锄头，至少得搞两个坑出来，把多余的“赘肉”给割掉（鼻子很简单，不过就是广阔的平原上冒出一个小山丘，而且和立脸颊面部连接非常光滑，没有凹凸拐点，因此，应该不影响咱们的网格分布构造）。OK，到目前为止，美容手术之前的准备工作已经结束，可以躺到手术床上了。 图3 整体拓补 图4 面部拓补 上了手术台（导入ICEM）,当然首先要做的工作便是一系列“宽衣解带”工作，ICEM里面这个过程称之为几何清理和修复，当然，这一切过程都基本上是自动化的，也有需要手工清理的时候，那种情况的前提是“相貌”长得实在太丑（几何模型太“烂”），不手工打底，实在没有办法。这个头像的几何模型还可以，基本上没有什么多余操作，自动清理就是了，这个过程实在没有啰嗦的必要，大多数网格程序都差不多。接下来，先回顾一下几种典型的几何拓补关系的分块方法（图5），彻底理解了这些分块思路，这个模型就简单了，图6为分块组合方式（几乎是高质量结构化六面体网格的必须方式）——从左到右，依次为圆球挖孔、1/8圆球+1/4圆柱、1/4圆球+1/2圆柱、圆柱、圆球、1/2圆柱+方体，当然，圆球也可异化为椭球体，方体可以异化为梯形体、楔形体、1/8圆球体可以异化为三棱锥体，这些异化体的几何拓补以及分块组合方式是完全相同的，因此图5的几何组合元素具有相当的普适性。图6、7是对应的分块拓补和网格，值得一提的是，最后的网格质量是相当地高。 图5 典型的几何拓补组合形式 图6 典型的几何分块拓补 图7 网格示意好了，言归正传，基于上述分块原则和方法，下面给出本案例的分块示意图（图8、9、10） (a)底座分块(b)底座+颈子 （c）整体分块 （d）头部分块局部 （e）头部局部分块图8 案例分块block好了 ，经过分块block，可以进行网格试划分了，经过初步调整vertices和edge,得到最低网格质量（2x2x2行列式检验标准）为0.24的网格，再次经过网格优化（全自动过程），网格质量达到令人满意的0.31，美容手术顺利结束。 （a） (b) (c) 图9 最终网格示意图美容手术顺利结束，必要的经验总结是不可缺少的，这个案例和F6飞机流场模型据说代表了“网格”剖分的典型技术难点（本案例属于腔体域内网格，F6则属于腔体域外绕流网格，个人认为F6短舱绕流模型应该更复杂一些），因此，本案例应该代表了复杂结构化网格划分的绝大部分方法与技巧，经验总结如下：1、合理的几何分块规划和识别能力是成功进行结构化网格的前提和必要能力，这个谁也教不了，得靠自己学习和悟性，最根本的分块拓补形式为方体，所有的其它分块形式只是其变异处理或者组合处理的结果，在高级网格工程师眼里，所有的几何形状均都应该看作是方块的堆砌，水平的差异主要体现在对方块的“发现能力”。 2、Y型block、C型block和O型block是分块组合常用的形式，目的一个：提高网格交接处（如x\y两向相交位置）几何边界夹角过小或者夹角圆弧型过渡造成的网格角度过小的问题，对于网格楔率过大或者扭曲、翘曲严重，也是有效解决手段。 3、对于个别位置的复杂分块，与其用从高到低的“砍块”的方式进行分块划 ...