Gmsh 中空心圆环的 3-D 网格划分问题

Question

我想创建一个空心或厚半环面，为此，我使用以下命令：

SetFactory ("OpenCASCADE");
Torus (1) = {0,0,0, 170,30, Pi};
Torus (2) = {0,0,0, 170,20, Pi};
BooleanDifference (8) = {Volume {1}; Delete; } {Volume {2}; Delete; };

当我尝试创建 3D 网格时，出现以下错误：

PLC Error: A segment and a facet intersect at point
Info: (122,229,106,391, -9.48334).
Info: Segment: [314,311] # -1 (0)
Info: Facet: [7,54,60] # 1
Error: Invalid boundary mesh (segment-facet intersection) on surface 1, intersection (122.229,106.391, -9.48334)
Error: No elements in volume 8

为什么会这样？怎么解决？

Answer 1

您的几何定义没有问题。一切都是合法的。

但是，在 GEO 文件中，您没有为网格指定所需的单元尺寸。而且，在这种特殊情况下，GMSH 无法使用它选择的默认四面体网格创建适当的四面体网格。

以下将允许您在空心半环面上创建适当的四面体网格：

Mesh.CharacteristicLengthMin = 5;
Mesh.CharacteristicLengthMax = 10;


SetFactory ("OpenCASCADE");
Torus (1) = {0,0,0, 170,30, Pi};
Torus (2) = {0,0,0, 170,20, Pi};
BooleanDifference (8) = {Volume {1}; Delete; } {Volume {2}; Delete; };

在这里，我手动指定了最小和最大元素大小。我的选择是任意的，主要是通过获得视觉上吸引人的网格来决定的。

您可以在the corresponding section of GMSH documentation 中阅读有关指定网格元素大小的各种选项（以及关于网格元素大小本身）的更多信息。