通过Gmsh3D在fipy中加载(.msh)文件时出现问题

2024-09-30 22:11:48 发布

您现在位置:Python中文网/ 问答频道 /正文

我正在尝试将网格文件(.msh)从Gmsh 4.7.0导入FiPy==3.4.2.1(Windows 10)(Python 2.7.9)。 该网格是使用开源FEM软件(Cast3M)创建的,然后作为.MED文件导出到Gmsh。Gmsh在相同数量的元素和节点下,可以完美地查看我的网格,没有任何问题。 然后,我继续从Gmsh将网格导出为.msh文件。(文件----导出----msh文件----保存)。 我在做:

mesh=fp.Gmsh3D('cylinder.msh') #After importing fipy

所以我得到了这个错误:

Traceback (most recent call last):
  File "<itasca-code-2-53a6529af1f3>", line 8, in <module>
  File "C:\Users\Youssef Fawaz\AppData\Roaming\Python\Python27\site-packages\fipy\meshes\gmshMesh.py", line 1993, in __init__
    self._orderedCellVertexIDs_data) = self.mshFile.read()
  File "C:\Users\Youssef Fawaz\AppData\Roaming\Python\Python27\site-packages\fipy\meshes\gmshMesh.py", line 813, in read
    raise GmshException(errStr)
fipy.meshes.gmshMesh.GmshException: Gmsh hasn't produced any cells! Check your Gmsh code.

Gmsh output:
Info    : Running 'gmsh c:\users\yousse~1\appdata\local\temp\tmpoqqytw.geo -3 -nopopup -format msh2 -o c:\users\yousse~1\appdata\local\temp\tmpygqpw3.msh' [Gmsh 4.7.0, 1 node, max. 1 thread]
Info    : Started on Fri Nov 13 10:31:12 2020
Info    : Reading 'c:\users\yousse~1\appdata\local\temp\tmpoqqytw.geo'...
Info    : Done reading 'c:\users\yousse~1\appdata\local\temp\tmpoqqytw.geo'
Info    : Meshing 1D...
Info    : Done meshing 1D (Wall 0s, CPU 0s)
Info    : Meshing 2D...
Info    : Done meshing 2D (Wall 0s, CPU 0s)
Info    : Meshing 3D...
Info    : Done meshing 3D (Wall 0s, CPU 0s)
Info    : 0 nodes 0 elements
Info    : Writing 'c:\users\yousse~1\appdata\local\temp\tmpygqpw3.msh'...
Info    : Done writing 'c:\users\yousse~1\appdata\local\temp\tmpygqpw3.msh'
Info    : Stopped on Fri Nov 13 10:31:12 2020 (From start: Wall 0.00300002s, CPU 0s)

需要注意的是,我正在DEM代码(由Itasca开发的PFC3D)中与fipy耦合实现此网格。例如,问题是导入我的网格(不是我的代码问题)。 以下是我的文件(.msh、.MED)和Gmsh3D中网格图像的链接 files

多谢各位


Tags: 文件info网格localcpuusersappdatatemp