pygmsh合并带孔洞拉伸体失败，网格生成报错求因

核心原因 #

• 两个拉伸体的几何拓扑不匹配：带孔洞的体拉伸后，其与另一个体接触的面（包括孔洞内环）和对方的对应面是拓扑独立的元素（即使位置重合），gmsh无法自动关联两侧的拉伸节点，导致网格生成时找不到对应节点。
• 未通过布尔运算显式定义合并约束：简单创建两个体后直接合并，没有告诉gmsh哪些几何元素是重合的，系统无法处理拓扑冲突。

解决办法 #

1. 复用拓扑一致的几何元素创建拉伸体 #

不要分别创建两个位置相同但独立的面，而是复用同一个基础面来生成两个拉伸体，确保接触区域的拓扑完全一致：

import pygmsh

with pygmsh.geo.Geometry() as geom:
    # 创建带孔的基础面
    rect = geom.add_rectangle(0, 2, 0, 2, z=0.0, mesh_size=0.15)
    circle = geom.add_circle([1.0, 1.0, 0.0], 0.4, mesh_size=0.1)
    face_with_hole = geom.boolean_difference([rect], [circle])

    # 向上拉伸带孔面得到第一个体
    vol1 = geom.extrude(face_with_hole, [0, 0, 1.0])
    # 复用同一个矩形面向下拉伸得到第二个体
    vol2 = geom.extrude(rect, [0, 0, -1.0])

    # 显式执行布尔融合
    geom.boolean_union([vol1, vol2])

    mesh = geom.generate_mesh()
    mesh.write("merged_volume.msh")

2. 为布尔运算指定公差 #

如果几何元素存在微小位置偏差（比如浮点精度导致的不重合），在布尔运算中添加tolerance参数，帮助gmsh识别重合元素：

geom.boolean_union([vol1, vol2], tolerance=1e-6)

3. 先合并2D几何再整体拉伸（替代方案） #

如果两个体的接触面对应的2D面可以合并，先完成2D面的布尔运算，再拉伸成3D体，从根源避免拓扑冲突：

with pygmsh.geo.Geometry() as geom:
    # 创建两个2D面（一个带孔，一个完整）
    rect1 = geom.add_rectangle(0, 2, 0, 2, z=0.0, mesh_size=0.15)
    circle = geom.add_circle([1.0, 1.0, 0.0], 0.4, mesh_size=0.1)
    face_with_hole = geom.boolean_difference([rect1], [circle])
    rect2 = geom.add_rectangle(0, 2, 0, 2, z=0.0, mesh_size=0.15)

    # 合并2D面
    merged_face = geom.boolean_union([face_with_hole, rect2])
    # 拉伸合并后的面得到3D体
    vol = geom.extrude(merged_face, [0, 0, 2.0])

    mesh = geom.generate_mesh()
    mesh.write("merged_from_2d.msh")

4. 统一网格尺寸参数 #

确保所有几何元素的mesh_size一致，避免局部网格密度差异导致拉伸时节点无法对应。

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