python + gmsh生成颗粒状复合材料盒装内插杆模型的网格

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