卓越拓扑之路----锥体的结构化网格 I

说起来搞有限元已经有段时间了，深知空间离散对于仿真的重要性，在这方面下了不少功夫，不断的追求更高的网格品质。有些几何体看起来是那么简单，但是画起网格来却总是让我崩溃。锥体就是个典型！当然，对于简单的集合体，非结构化的网格不在考虑范围之内。

今天有些时间，干脆就研究研究三棱锥。













常规做法，先做CAD，然后建立一个BLOCK包围离散集合体。



















对于三棱锥而言，底面只有三个节点，所以需要对BLOCK其中一面上的节点进行合并。为了保证网格的连贯性，索性我把相对应的BLOCK上表面的节点合并了。




















把相应的节点跟几何体相关联以后，再用一个简单得Y-BLOCK，差不多就可以设定网格的大小了。

记得我再做渗透分析的时候，模型的每条边上我都定义了很多的网格，以便计算收敛的好些，模拟的更贴近实际，结果看到了千丝万缕的空间，如同万箭穿心一般。为了不至于眼花缭乱，对于这个集合体，我只定义了有限几个单元。

现在问题已经很简单了，锥体上端只有一个节点，只有BLOCK上表面的所有节点都合并起来，问题似乎就应该解决了。

上面的网格，看似漂亮，结构清晰明了，质量很高。但是没用，最后的一步是一直以来困扰我的难题。虽然，我可以把节点合并，却最终改变了最顶端网格的性质，违背了结构化网格的基本定义。冥思苦想一上午，还是没有想到好地解决办法。

难道锥体很难建立高质量的结构化空间离散吗？

未完，待续。。。