1. 前期准备:
导入.vtk模型文件,内置一系列三维点。
先在CSV文件生成的点云数据,然后使用像ParaView或带有VTK库的Python软件包将其转换为VTK格式并保存为一个.vtk文件。
2. python代码读取模型,为gmsh创建新模型。
2.1 代码初始化新的Gmsh模型,名称为 "ConcreteBar"。
model_name="ConcreteBar"
gmsh.initialize(sys.argv)
gmsh.model.add(model_name)
gmsh.initialize 函数初始化 Gmsh API,并将传递给脚本的命令行参数作为一个参数 (在本例中是 sys.argv) 。 gmsh.model.add 函数用给定的名字创建一个新的模型。所有对 Gmsh API 的后续调用都将适用于这个模型,直到有新的模型被添加或当前的模型被删除。 因此,运行这段代码后,使用Gmsh函数为 "ConcreteBar "模型定义几何体、网格和物理组。
2.2 读取.vtk文件,初始化颗粒数据集。
fname="/home/lalith/Downloads/GMSH/Concrete/PoissonSampling/Poisson_sample.vtk"
m = meshio.Mesh.read(fname, "vtk") # same arguments as meshio.read
num_spheres=len(m.points) #reset the number of spheres
print("number of random pints in fname ",num_spheres)
SPH_info=np.zeros((num_spheres, 4))
#for pt in m.points:
# print(pt[0])
for i in range (0,num_spheres):
SPH_info[i] = [m.points[i][0], m.points[i][1], m.points[i][2], 0.5*RS]
#//
m = meshio.Mesh.read(fname, "vtk")
meshio被导入,读取函数被调用,有两个参数:文件名 "mesh.vtk "和字符串 "vtk",它指定了要使用的文件格式。运行这段代码后,m将是一个Mesh对象,它包含了从VTK文件中读取的网格的节点和元素。可以使用m.points和m.cell属性来分别访问Mesh节点和元素。
num_spheres=len(m.points): 将num_spheres的值重置为等于Mesh中的点数。
SPH_info=np.zeros((num_spheres, 4))创建了一个名为SPH_info的二维NumPy数组,