Gmsh有限元网格剖分(Python)---点的创建、连点成线，围线成面

已于 2025-01-12 22:44:25 修改
阅读量1.6k 收藏 13
点赞数 10
分类专栏： 电磁计算 网格剖分 文章标签： python 平面 开发语言
于 2024-12-18 22:10:13 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44584198/article/details/144569131
版权

Gmsh有限元网格剖分(Python)—点的创建、连点成线，围线成面

最近在学习有限元的网格剖分算法，主要还是要参考老外的开源Gmsh库进行，写一些博客记录下学习过程，方便以后回忆嘞。

Gmsh的官方英文文档可以参考：gmsh.pdf

但咋就说，这里面东西太多了，不太适合初学者去看，因此分享下自己学习的过程吧。

目录

0、Gmsh的Python接口基础

使用Gmsh的Python接口一般有下面几个步骤：

Step1：导入包

import gmsh
import sys

Step2：初始化gmsh引擎、给整个结构命名(此处命名为my_t1)

gmsh.initialize()
gmsh.model.add("my_t1")

Step3：设置创建结构使用的引擎
------gmsh.model.geo是GMSH自带的几何定义接口，适用于简单的几何定义和修改
------gmsh.model.occ基于Open CASCADE，模型修改更加高级

# gmsh.model.geo是GMSH自带的几何定义接口，适用于简单的几何定义和修改；
# gmsh.model.occ基于Open CASCADE，模型修改更加高级
# kernel = gmsh.model.geo
kernel = gmsh.model.occ

Step4：创建结构，如点、线、面等，此处给出例子

mesh_size = 1e-2
point1_tag = kernel.addPoint(0, 0, 0, meshSize=mesh_size)
point2_tag = kernel.addPoint(0, 1, 0, meshSize=mesh_size)
point3_tag = kernel.addPoint(1, 0, 0, meshSize=mesh_size)
point4_tag = kernel.addPoint(1, 1, 0, meshSize=mesh_size)

Step5：模型同步，并将结构赋予实际的物理分组

# GMSH模型同步
kernel.synchronize()
# 创建点集合，命名为My point
gmsh.model.addPhysicalGroup(0, [point1_tag, point2_tag, point3_tag, point4_tag], name="My point")

Step6：进行剖分、保存、调用gui显示、退出

# 2维网格剖分
gmsh.model.mesh.generate(2)
# 保存一下
gmsh.write("my_t1.msh")
# 默认启动gui界面
if '-nopopup' not in sys.argv:
    gmsh.fltk.run()
# 结束gmsh
gmsh.finalize()

1、Gmsh中点的创建

gmsh.model.geo.addPoint函数可以进行点的创建，函数的解释如下：

  1. 前三个参数是点的坐标（x, y, z） 输入mesh_size是目标网格大小，它指定了在点附近期望的网格尺寸
  2. 可以手动设置创建的点的标签，如gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size, tag=1)
  3. 不手动设置标签，创建时会自动返回标签，如point1_tag = gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size)

实际的使用代码：

point1_tag = kernel.addPoint(0, 0, 0, meshSize=mesh_size)
point2_tag = kernel.addPoint(0, 1, 0, meshSize=mesh_size)
point3_tag = kernel.addPoint(1, 0, 0, meshSize=mesh_size)
point4_tag = kernel.addPoint(1, 1, 0, meshSize=mesh_size)

# GMSH模型同步
kernel.synchronize()
# 创建点集合，命名为My point
gmsh.model.addPhysicalGroup(0, [point1_tag, point2_tag, point3_tag, point4_tag], name="My point")

运行代码可以调用gmsh的UI界面进行显示，可以查看到四个点被成功创建：
在这里插入图片描述
双击空白位置，打开坐标轴：
在这里插入图片描述
结构更加清晰明了了：
在这里插入图片描述
在菜单栏的Tool中可以查看结构中各个元素的可视化选项：
在这里插入图片描述
My point元素组刚好对应四个点的集合：
在这里插入图片描述
全部代码：

import gmsh
import sys

gmsh.initialize()
# 新建模型，命名为t1，t1指的是整个结构
gmsh.model.add("my_t1")

# gmsh.model.geo是GMSH自带的几何定义接口，适用于简单的几何定义和修改；
# gmsh.model.occ基于Open CASCADE，模型修改更加高级
# kernel = gmsh.model.geo
kernel = gmsh.model.occ

mesh_size = 1e-2
# 前三个参数是点的坐标（x, y, z）
# 输入mesh_size是目标网格大小，它指定了在点附近期望的网格尺寸
# 可以手动设置创建的点的标签，如gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size, tag=1)
# 不手动设置标签，创建时会自动返回标签，如point1_tag = gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size)
point1_tag = kernel.addPoint(0, 0, 0, meshSize=mesh_size)
point2_tag = kernel.addPoint(0, 1, 0, meshSize=mesh_size)
point3_tag = kernel.addPoint(1, 0, 0, meshSize=mesh_size)
point4_tag = kernel.addPoint(1, 1, 0, meshSize=mesh_size)

# GMSH模型同步
kernel.synchronize()
# 创建点集合，命名为My point，第一个输入是0代表是0维的点
gmsh.model.addPhysicalGroup(0, [point1_tag, point2_tag, point3_tag, point4_tag], name="My point")

# 2维网格剖分
gmsh.model.mesh.generate(2)
# 保存一下
gmsh.write("my_t1.msh")
# 默认启动gui界面
if '-nopopup' not in sys.argv:
    gmsh.fltk.run()
# 结束gmsh
gmsh.finalize()

2、Gmsh中线的创建-连点成线

使用下面的案例代码，line1是点1和点2相连得到，line2是点2和点4相连得到，其余类似：

line1_tag = kernel.addLine(point1_tag, point2_tag)
line2_tag = kernel.addLine(point2_tag, point4_tag)
line3_tag = kernel.addLine(point4_tag, point3_tag)
line4_tag = kernel.addLine(point3_tag, point1_tag)

使用下面的代码创建物理集合，第一个输入参数是1代表输入的是1维的线：

gmsh.model.addPhysicalGroup(1, [line1_tag, line2_tag, line3_tag, line4_tag], name="My line")

创建后运行代码，可以看到gmsh识别到全部的点和线：
在这里插入图片描述
全部代码：

import gmsh
import sys

gmsh.initialize()
# 新建模型，命名为t1，t1指的是整个结构
gmsh.model.add("my_t1")

# gmsh.model.geo是GMSH自带的几何定义接口，适用于简单的几何定义和修改；
# gmsh.model.occ基于Open CASCADE，模型修改更加高级
# kernel = gmsh.model.geo
kernel = gmsh.model.occ

mesh_size = 1e-2
# 前三个参数是点的坐标（x, y, z）
# 输入mesh_size是目标网格大小，它指定了在点附近期望的网格尺寸
# 可以手动设置创建的点的标签，如gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size, tag=1)
# 不手动设置标签，创建时会自动返回标签，如point1_tag = gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size)
point1_tag = kernel.addPoint(0, 0, 0, meshSize=mesh_size)
point2_tag = kernel.addPoint(0, 1, 0, meshSize=mesh_size)
point3_tag = kernel.addPoint(1, 0, 0, meshSize=mesh_size)
point4_tag = kernel.addPoint(1, 1, 0, meshSize=mesh_size)

line1_tag = kernel.addLine(point1_tag, point2_tag)
line2_tag = kernel.addLine(point2_tag, point4_tag)
line3_tag = kernel.addLine(point4_tag, point3_tag)
line4_tag = kernel.addLine(point3_tag, point1_tag)


# GMSH模型同步
kernel.synchronize()
# 创建点集合，命名为My point
gmsh.model.addPhysicalGroup(0, [point1_tag, point2_tag, point3_tag, point4_tag], name="My point")

gmsh.model.addPhysicalGroup(1, [line1_tag, line2_tag, line3_tag, line4_tag], name="My line")

# 2维网格剖分
gmsh.model.mesh.generate(2)
# 保存一下
gmsh.write("my_t1.msh")
# 默认启动gui界面
if '-nopopup' not in sys.argv:
    gmsh.fltk.run()
# 结束gmsh
gmsh.finalize()

3、Gmsh中面的创建-连线成面

首先要把线的集合变成封闭的曲线（CurveLoop函数）：

line1_tag = kernel.addLine(point1_tag, point2_tag)
line2_tag = kernel.addLine(point2_tag, point4_tag)
line3_tag = kernel.addLine(point4_tag, point3_tag)
line4_tag = kernel.addLine(point3_tag, point1_tag)

curve_loop_tag = kernel.addCurveLoop([line1_tag, line2_tag, line3_tag, line4_tag])

随后从封闭的曲线创建平面：

surface_tag = kernel.addPlaneSurface([curve_loop_tag])

创建的封闭线、平面全部全部添加到实际物理Group：

gmsh.model.addPhysicalGroup(1, [curve_loop_tag], name="My CurveLoop")
gmsh.model.addPhysicalGroup(2, [surface_tag], name="My Surface")

运行代码，案例如下：
在这里插入图片描述
全部代码如下：

import gmsh
import sys

gmsh.initialize()
# 新建模型，命名为t1，t1指的是整个结构
gmsh.model.add("my_t1")

# gmsh.model.geo是GMSH自带的几何定义接口，适用于简单的几何定义和修改；
# gmsh.model.occ基于Open CASCADE，模型修改更加高级
# kernel = gmsh.model.geo
kernel = gmsh.model.occ

mesh_size = 1e-1
# 前三个参数是点的坐标（x, y, z）
# 输入mesh_size是目标网格大小，它指定了在点附近期望的网格尺寸
# 可以手动设置创建的点的标签，如gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size, tag=1)
# 不手动设置标签，创建时会自动返回标签，如point1_tag = gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size)
point1_tag = kernel.addPoint(0, 0, 0, meshSize=mesh_size)
point2_tag = kernel.addPoint(0, 1, 0, meshSize=mesh_size)
point3_tag = kernel.addPoint(1, 0, 0, meshSize=mesh_size)
point4_tag = kernel.addPoint(1, 1, 0, meshSize=mesh_size)

line1_tag = kernel.addLine(point1_tag, point2_tag)
line2_tag = kernel.addLine(point2_tag, point4_tag)
line3_tag = kernel.addLine(point4_tag, point3_tag)
line4_tag = kernel.addLine(point3_tag, point1_tag)

curve_loop_tag = kernel.addCurveLoop([line1_tag, line2_tag, line3_tag, line4_tag])

surface_tag = kernel.addPlaneSurface([curve_loop_tag])


# GMSH模型同步
kernel.synchronize()
# 创建点集合，命名为My point
gmsh.model.addPhysicalGroup(0, [point1_tag, point2_tag, point3_tag, point4_tag], name="My point")

gmsh.model.addPhysicalGroup(1, [line1_tag, line2_tag, line3_tag, line4_tag], name="My line")

gmsh.model.addPhysicalGroup(1, [curve_loop_tag], name="My CurveLoop")

gmsh.model.addPhysicalGroup(2, [surface_tag], name="My Surface")

# 2维网格剖分
gmsh.model.mesh.generate(2)
# 保存一下
gmsh.write("my_t1.msh")
# 默认启动gui界面
if '-nopopup' not in sys.argv:
    gmsh.fltk.run()
# 结束gmsh
gmsh.finalize()

4、简单调整平面剖分精度

gmsh.model.geo.addPoint函数可以进行点的创建，函数的解释如下：

  1. 前三个参数是点的坐标（x, y, z） 输入mesh_size是目标网格大小，它指定了在点附近期望的网格尺寸
  2. 可以手动设置创建的点的标签，如gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size, tag=1)
  3. 不手动设置标签，创建时会自动返回标签，如point1_tag = gmsh.model.geo.addPoint(x=0, y=0, z=0, meshSize=mesh_size)

注意到，调整函数的meshSize输入参数可以设置网格精度，上面最后代码按照0.1进行剖分，如果我们修改上面代码设置：

mesh_size = 5e-1

结果为：
在这里插入图片描述

用gmsh创建一个立方体、物理分组、网格划分
m0_68889767的博客
06-24 2954
点击Physical groups下的surface(这里我选择的是对面进行分组），会跳出以下界面，输入你要设置的物理组名称，我设置的是surface，在模型中选中你需要的面，选择完毕后按下e，物理组即可设置完成，如果所需的物理组大于等于2，后面的设置方法也是如此。为避免创建面的时候出现忘记那个面已经创建过，所以建议创建的线足够构成一个一个面时即可创建一个面，然后再接着创建线，可以构成一个面时再创建一个面，依次类推，直至创建完所有的面。点击Mesh下的2D（因为这里我用到的是2D网格，按自己所需进行设置）
gmeh搭配python建模，嵌入式模型，学习心得
fum1125的博客
04-06 907
当使用语法my_list[index]访问一个列表中的元素时，会得到列表中索引位置的元素。在本例中，point 是一个包含单个元组的列表（4.905, 5.0075, 0）。因此，point[0] 返回元组（4.905, 5.0075, 0），point[0][0] 返回元组的第一个元素，也就是 4.905。建立两个实体，volume a 是盒子，从盒子中挖出矩形体 volume b，在盒子中嵌入颗粒 volume a。points[i][0]指的是列表point中第一个点的x坐标。
使用Gmsh给一个长方体划分规整的四边形网格
m0_61467485的博客
03-31 531
使用Gmsh给一个长方体划分规整的四边形网格。
pygmsh: Python下的Gmsh应用库
最新发布
gitblog_06717的博客
04-30 726
pygmsh: Python下的Gmsh应用库 【下载地址】pygmshPython下的Gmsh应用库 pygmsh是一个专为Python开发者设计的开源库，它巧妙地将Gmsh的强大功能与Python的灵活性相结合。通过提供对Gmsh Python接口的高效抽象，pygmsh使得创建复杂几何图形变得异常简便。无论是进行科...
Gmesh的网格剖分模块
Chinamming的专栏
11-22 9997
第6章 网格剖分模块 Gmsh的网格剖分模块集中了几种1维，2维，3维的网格剖分算法，他们都能产生适合于有限元的网格： 2维的非结构化算法产生三角单元或者三角单元与四方形（当使用拼接曲面时）。3维的非结构化算法产生四面体单元。 3维的结构化算法（超限或拉伸）默认产生三角单元，但是通过recombine命令，能得到四方形。3维结构化算法产生四面体，六面体，三菱柱、椎体，取决于曲面网格的类型。
python之gmsh划分网格
T20151470的博客
04-15 2396
对于Python用户，可以使用Pygmsh，这是Gmsh的Python接口，它改进了直接使用Gmsh的一些缺点，例如在Gmsh脚本中每个几何实体必须手动分配唯一编号，而Pygmsh允许使用变量名来命名这些几何实体，避免了对编号唯一性的检查。在使用Gmsh时，用户可以通过设置基本的几何实体和布尔操作来建立几何形状，然后进行网格划分。Gmsh遵循GNU通用公共许可证（GPL）的条款进行分发，这意味着任何人都可以自由使用和分发Gmsh，但不能将Gmsh集成到任何计划分发的封闭源代码软件中，除非获得不同的许可。
Gmsh剖一维网格教程附代码
guihunkun的博客
09-20 2844
文章目录前言一、点，线段，单元，多边形等的定义1. 点的定义2. 多边形的定义3. 线段的定义4. 拓扑关系的定义5. 单元的定义二、区域构造三、调用Gmsh的API剖网格并得到网格信息四、程序运行结果1. 程序输出2. 带有点标号的结果图3. 带有单元标号的结果图五、总结 前言 Gmsh的安装与入门教程可以去看我的前一篇文章有限元剖网格之Gmsh安装与使用入门。本篇文章主要介绍如何借助Gmsh的C++ API剖一维网格(即线段)。 一、点，线段，单元，多边形等的定义 因为我们需要自己构造区域内的几何信
Gmsh有限元网格剖分(Python)---点、直线、平面的移动和旋转
weixin_44584198的博客
12-22 1226
kernel = gmsh.model.geo时对形状的编辑耦合比较严重kernel = gmsh.model.occ时各个图形隔离度高，推荐使用occ引擎。
Gmsh教程
weixin_44897632的博客
07-18 509
【代码】Gmsh教程。
meshpy:适用于Python的2D3D简单网格生成器接口（Triangle，TetGen，gmsh）
03-08
这种类型的网格主要用于有限元仿真代码中，但也有许多其他应用，从计算机图形学到机器人技术。 为了产生二维和三维网格，MeshPy提供的Python接口到三个广泛关注的网格生成器，由J. Shewchuk， 由杭司前二者都包含...
python有限元网格剖分
02-22
### Python 中有限元网格剖分的方法与库 #### 使用 `meshpy` 进行三角形和四边形单元的网格生成 对于二维结构，可以采用 `MeshPy` 来创建高质量的三角形或四边形单元网格。该工具包提供了简单易用的应用编程接口...
pygmsh:适用于Python的Gmsh
04-28
适用于Python的Gmsh。 pygmsh将的功能与Python的多功能性结合在一起。 它提供了Gmsh自己的Python界面的有用抽象，因此您可以更轻松地创建复杂的几何图形。 要使用，请从安装Gmsh本身和pygmsh： [sudo] apt install python3-gmsh pip install pygmsh 本文档和目录包含许多小示例。 有关完整的文档，请参见。 平面形状 多边形 圆圈 （B-）花键 代码： import pygmsh with pygmsh . geo . Geometry () as geom : geom . add_polygon ( [ [ 0.0 , 0.0 ], [ 1.0 , - 0.2 ], [ 1.1 , 1.2 ],
1750260_gmsh-2.5.0-source_2.5_源码_gmsh_
09-30
gmsh不仅支持创建复杂的几何模型，还具备强大的网格剖分能力，能够生成高质量的有限元网格。它集成了用户界面、几何建模、网格化和后处理功能，适用于各种数值模拟的预处理工作。gmsh的源码开放，使得开发者可以深入...
tetgen：C ++ TetGen库的Python接口，用于生成任何3D多面体域的四面体网格
02-01
它将复杂的几何形状转化为离散的、结构化的单元，这些单元可以被进一步用于有限元分析、流体动力学模拟或其他数值方法。四面体网格因其良好的数学性质（如局部和全局的刚性）而被广泛采用，尤其适用于三维问题。 ...
有限元网格开源软件GMSH基础
Aaron
11-19 9138
脚本式建模例子 脚本建模的语法网上较多，不做介绍，附一个例子，帮助快速入门。 //////////dimensions of three cylinders r1=21.09; l1=27; r2=15; l2=4.68; r3=27; l3=20; lc=10; //////////
有限元剖网格之Gmsh安装与使用入门
热门推荐
guihunkun的博客
09-12 2万+
有限元剖网格之Gmsh安装与使用入门 有限元法（FEM，Finite Element Method）是一种求解偏微分方程问题的数值方法。随着电子计算机的发展，有限元方法迅速发展成一种现代计算方法，在固体力学、流体力学、热传导、电磁学、声学、生物力学等方面有着广泛的应用。今天就先介绍有限元方法之前处理，利用Gmsh剖网格。 Gmsh介绍 Gmsh是一个自动的三维有限元网格生成带有内置在CAD和后期处理器。当然其也可以生成一维和二维的有限元网格。 Gmsh安装 1. Gmsh图形界面安装 可以去官网下载安装包G
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

怡步晓心l

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值