mesh - 在 PyGmsh / Gmsh 中添加内部边界或裂缝

Question

我正在尝试使用以下代码使用 PyGmsh 生成有限元网格：

import pygmsh

geom = pygmsh.opencascade.Geometry(
    characteristic_length_min=0.1,
    characteristic_length_max=0.1,
    )

rectangle = geom.add_rectangle([-1.0, -1.0, 0.0], 2.0, 2.0)
disk1 = geom.add_disk([-1.2, 0.0, 0.0], 0.5)
disk2 = geom.add_disk([+1.2, 0.0, 0.0], 0.5)
disk3 = geom.add_disk([0.0, -0.9, 0.0], 0.5)
disk4 = geom.add_disk([0.0, +0.9, 0.0], 0.5)

union = geom.boolean_union([rectangle, disk1, disk2])
diff = geom.boolean_difference([union], [disk3, disk4])

mesh = pygmsh.generate_mesh(geom, dim=2)

我可以生成以下网格：

但是，我想在网格上添加一个裂缝，例如：

这里的裂缝只是一个例子，它需要在网格划分之前定义。我尝试创建 2 个点 ( geom.add_point()) 和一条线 ( geom.add_line())，然后 geom.boolean_difference()在最终的几何图形和线/裂缝之间做一个，但这不起作用。

任何帮助将不胜感激。

编辑

这种类型的网格生成的目的是模拟身体中的物理裂缝。在网格划分过程中，可以通过网格的元素连通性对裂纹进行建模（即元素必须具有不同的节点才能创建裂纹面）。例如，在施加任何载荷之前，裂缝是闭合的：

施加载荷后，裂缝打开，因为单元连通性允许这样做：

标签： meshtriangulationfinite-element-analysisopencascade

Answer 1

您可以通过在该区域建模一个非常窄的矩形来实现这一点。你可以很容易地给出像 1e-10 这样的尺寸。我还模拟了裂纹尖端，通过模拟一个非常小的圆圈将节点折叠在一个点上。它工作得很好。

现在也有一个插件。它会自动分离指定裂缝线/表面处的节点。