Gmesh的网格剖分模块

第6章 网格剖分模块

Gmsh的网格剖分模块集中了几种1维，2维，3维的网格剖分算法，他们都能产生适合于有限元的网格：

2维的非结构化算法产生三角单元或者三角单元与四方形（当使用拼接曲面时）。3维的非结构化算法产生四面体单元。

3维的结构化算法（超限或拉伸）默认产生三角单元，但是通过recombine命令，能得到四方形。3维结构化算法产生四面体，六面体，三菱柱、椎体，取决于曲面网格的类型。

通过使用Mesh.SubdivisionAlgorihm选项，所有的网格都可再细分为完全的四方形或六面体单元。但是，要注意到是，非结构化方法产生的初始网格，经过细化后，单元的质量可能非常差。

6.1 选择正确的非结构化算法

Gmsh有不同的2维的非结构化算法和不同3维的非结构化算法可供选择。每种算法有它的优点和缺点。对于所有的2维非结构化算法，会首先使用分治法产生初始的Delaunay网格，包含了所有的点。使用边交换的方法，调整不符合delaunay特性的边。之后，可使用3个不同的算法来产生最终的网格：

1、网格自适应算法，基于局部网格调整。这个技巧使用边交换、分割、叠合：长的边被分割，短的边被叠合，如果交换边后得到更好的网格就交换边。

2、Delaunay算法。依次对具有最大外接圆的单元，插入新的点在其外接圆心上。然后重新连接网格，使用各向异质的Delaunay准则。

3、波前法。

各个算法特性列表如下：

自动化程度 性能 单元质量</