0

我正在尝试在 FiPy 中导入使用 Gmsh 生成的 3D 网格。使用 2D 网格进行的测试效果很好。如果随后使用 Gmsh3D 挤压和导入模型，我会收到一条错误消息。

GmshException: Gmsh 没有产生任何细胞！检查您的 Gmsh 代码。

我正在使用 Python 3.7.3、Fipy 3.1.3 和 Gmsh 3.0.6（如推荐）在 Win10 上工作。

test2D.geo 测试文件：

SetFactory("OpenCASCADE");
cl = 0.5;
bs = 2.;
Point(1) = {0, 0, 0, cl};
Point(2) = {0, bs, 0, cl};
Point(4) = { bs,  0, 0, cl};
Point(3) = {bs,  bs, 0, cl};
Line(5) = {1, 2};
Line(6) = {2, 3};
Line(7) = {3, 4};
Line(8) = {4, 1};

Line Loop(10) = {6, 7, 8, 5};
Plane Surface(1) = {10};
Extrude {0, 0, 1} {
  Surface{1}; 
}

和：

from fipy import *
mesh = Gmsh3D("test2D.msh")

错误信息是：GmshException: Gmsh 没有产生任何细胞！检查您的 Gmsh 代码。

我没有看到我的错误，希望有人可以在这里帮助我。提前致谢

为 Gmsh 输出编辑：

Gmsh output:
Info    : Running 'gmsh C:\Users\Tinka\AppData\Local\Temp\tmpj4zr8g_c.geo -3 -nopopup -format msh -o C:\Users\Tinka\AppData\Local\Temp\tmpnz1bp4vu.msh' [Gmsh 3.0.6, 1 node, max. 1 thread]
Info    : Started on Tue May 28 19:50:42 2019
Info    : Reading 'C:\Users\Tinka\AppData\Local\Temp\tmpj4zr8g_c.geo'...
Info    : Done reading 'C:\Users\Tinka\AppData\Local\Temp\tmpj4zr8g_c.geo'
Info    : Finalized high order topology of periodic connections
Info    : Meshing 1D...
Info    : Done meshing 1D (0 s)
Info    : Meshing 2D...
Info    : Done meshing 2D (0 s)
Info    : Meshing 3D...
Info    : Done meshing 3D (0 s)
Info    : 0 vertices 0 elements
Info    : Writing 'C:\Users\Tinka\AppData\Local\Temp\tmpnz1bp4vu.msh'...
Info    : Done writing 'C:\Users\Tinka\AppData\Local\Temp\tmpnz1bp4vu.msh'
Info    : Stopped on Tue May 28 19:50:42 2019
4

2 回答 2

0

我将参数的名称更改Gmsh3D为test2D.geo并从地理文件中删除了第一行，一切似乎都在工作。

>>> from fipy import Gmsh3D
>>> mesh = Gmsh("test2D.geo")
>>> print(mesh.cellCenters)
[[1.34821429 1.24404762 1.34821429 ...
...

我不确定第一行是做什么的，但是我得到Error : Gmsh requires OpenCASCADE to add vertex并且如果包含它，则不会生成任何顶点或单元格，但生成网格不是必需的。

我认为 FiPy Gmsh 类同时采用 geo 和 msh 格式的文件，但文件名确实需要引用磁盘上的实际文件。

我正在使用 FiPy 版本 3.2+2.gccec299e 和 Gmsh 版本 3.0.6。

于 2019-05-28T15:43:11.673 回答
0

gmsh 和 spyder 的这个问题已在今天早些时候发布的 FiPy 3.3 中得到修复；谢谢你的报告。

您在聊天中报告的另一个问题是不同的。正如Gmsh2D 所记录的，但不是 Gmsh3D：

... // attention: if you use any "Physical" labels, you *must* label
... // all elements that correspond to FiPy Cells (Physical Surface in 2D
... // and Physical Volume in 3D) or Gmsh will not include them and FiPy
... // will not be able to include them in the Mesh.
...
... // note: if you do not use any labels, all Cells will be included.

添加Physical Volumes("cells") = {1};到您的.geo脚本将解决该问题。

于 2019-07-01T23:36:09.120 回答