如何使用gmsh或pygmsh加载并检查*.msh网格文件？

你完全不需要依赖meshio，用Gmsh官方的Python绑定或者pygmsh就能实现读取*.msh文件并提取节点、单元、边、面信息的需求，具体方法如下：

Gmsh自带的Python绑定完全支持读取网格文件，且能直接获取各类网格实体的详细信息，步骤如下：

1. 初始化Gmsh并读取文件

import gmsh

# 初始化Gmsh环境
gmsh.initialize()
# 读取目标msh文件
gmsh.open("your_mesh_file.msh")

2. 获取节点数量

# 获取所有节点的标签、坐标（第三个返回值是节点所属物理组，这里用不到）
node_tags, _, _ = gmsh.model.mesh.getNodes()
node_count = len(node_tags)
print(f"节点总数: {node_count}")

3. 统计单元、边、面的数量
   Gmsh中不同维度的实体对应不同的单元类型：1D单元是边，2D是面，3D是体单元。我们可以通过单元类型的维度属性来筛选统计：

# 获取所有存在的单元类型
element_types = gmsh.model.mesh.getElementTypes()

total_element_count = 0
edge_count = 0
face_count = 0

for elem_type in element_types:
    # 获取单元类型的维度（0=点,1=边,2=面,3=体）
    elem_dim, _, _, _, _, _ = gmsh.model.mesh.getElementProperties(elem_type)
    # 获取该类型的所有单元标签
    elem_tags, _ = gmsh.model.mesh.getElements(elem_type)
    
    if elem_tags:
        count = len(elem_tags[0])
        total_element_count += count
        
        if elem_dim == 1:
            edge_count += count
        elif elem_dim == 2:
            face_count += count

print(f"总单元数: {total_element_count}")
print(f"边的数量: {edge_count}")
print(f"面的数量: {face_count}")

4. 清理Gmsh环境

gmsh.finalize()

pygmsh虽然主打网格生成，但底层依赖Gmsh，也支持读取msh文件并提取基础信息，代码更简洁：

import pygmsh

with pygmsh.open("your_mesh_file.msh") as mesh:
    # 节点数量
    node_count = len(mesh.points)
    print(f"节点总数: {node_count}")
    
    # 总单元数：累加所有类型单元的数量
    total_element_count = sum(len(cells) for cells in mesh.cells.values())
    print(f"总单元数: {total_element_count}")
    
    # 边的数量：对应"line"类型单元
    edge_count = len(mesh.cells.get("line", []))
    print(f"边的数量: {edge_count}")
    
    # 面的数量：统计所有2D单元类型（如triangle、quad等）
    face_types = ["triangle", "quad", "polygon"]
    face_count = 0
    for ft in face_types:
        face_count += len(mesh.cells.get(ft, []))
    print(f"面的数量: {face_count}")

总结 #

Gmsh Python API功能更全面，能处理复杂的网格类型和物理组信息；pygmsh则胜在代码简洁，适合快速提取基础统计数据，可根据你的需求选择，完全不需要必须使用meshio。

内容的提问来源于stack exchange，提问作者abukaj