根据@Paddy Harrison的建议,我使用了
我有一组顶点和三角形面,它们一起形成一个三角形网格:
import numpy as np
verts = [[0.1, 1., 1. ] [1., 1., 0.1] [1., 0.1, 1. ] [1., 1., 1.9] [1., 1.9, 1. ]
[1.9, 1., 1. ] ]
faces = [[ 2, 1, 0] [ 0, 3, 2] [ 1, 4, 0] [ 0, 4, 3] [ 5, 1, 2] [ 3, 5, 2]
[ 5, 4, 1] [ 4, 5, 3]]
我知道从这些创建三角形网格的唯一方法是使用Poly3Dcollection
from mpl_toolkits.mplot3d.art3d import Poly3DCollection
myMesh = Poly3DCollection(verts[faces])
接下来,我想使用pygalmesh模块来创建体积网格。它应该接收一个曲面网格并输出一个体积网格,如图所示。
根据教程,我应该能够使用创建体积网格
import pygalmesh
mesh = pygalmesh.generate_volume_mesh_from_surface_mesh(
"elephant.vtu",
facet_angle=25.0,
facet_size=0.15,
facet_distance=0.008,
cell_radius_edge_ratio=3.0,
verbose=False
)
然而,当我运行时:
mesh = pygalmesh.generate_volume_mesh_from_surface_mesh(myMesh)
我被退回:
TypeError: expected str, bytes or os.PathLike object, not Poly3DCollection
我认为这是因为我创建的网格不正确,或者我打算在使用pygalmesh.generate_volume_mesh_from_surface_mesh
之前以某种形式保存它。我不确定。我已经在github模块中提出了这个问题,但还没有收到好的反馈。
meshio
的write
函数。具体来说,对于我的数据集,它是这样做的:
import meshio
points = np.array(verts)
cells = [("triangle", np.array(faces)]
meshio.write_points_cells('out.vtu',points,cells)