Gmsh教程

于 2024-07-18 05:00:00 发布
阅读量64 收藏
点赞数
分类专栏： Gmsh 文章标签： python Gmsh
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44897632/article/details/140492877
版权

13、在没有底层CAD模型的情况下重新擦除STL文件

import gmsh  # 导入Gmsh库，用于几何建模和网格划分  
import math  # 导入数学库，用于计算  
import os    # 导入操作系统库，用于处理文件路径  
import sys   # 导入系统库，用于处理命令行参数  
  
gmsh.initialize()  # 初始化Gmsh环境  
  
def createGeometryAndMesh():  
    # 清除之前的模型和数据  
    gmsh.clear()  
      
    # 获取当前脚本的目录，并合并上级目录中的't13_data.stl'文件路径  
    path = os.path.dirname(os.path.abspath(__file__))  
    gmsh.merge(os.path.join(path, os.pardir, 't13_data.stl'))  
  
    # 从Gmsh的ONELAB接口获取用户输入的角度参数  
    angle = gmsh.onelab.getNumber('Parameters/Angle for surface detection')[0]  
  
    # 获取是否创建可参数化曲面的用户输入  
    forceParametrizablePatches = gmsh.onelab.getNumber(  
        'Parameters/Create surfaces guaranteed to be parametrizable')[0]  
  
    # 设置包含边界和曲线角度的参数  
    includeBoundary = True  
    curveAngle = 180  
  
    # 对模型表面进行分类，准备创建几何体  
    gmsh.model.mesh.classifySurfaces(angle * math.pi / 180., includeBoundary,  
                                     forceParametrizablePatches,  
                                     curveAngle * math.pi / 180.)  
  
    # 创建几何体  
    gmsh.model.mesh.createGeometry()  
  
    # 获取所有曲面，并创建一个曲面环和体积  
    s = gmsh.model.getEntities(2)  
    l = gmsh.model.geo.addSurfaceLoop([e[1] for e in s])  
    gmsh.model.geo.addVolume([l])  
  
    # 同步几何模型  
    gmsh.model.geo.synchronize()  
  
    # 添加一个网格大小字段，根据用户选择应用不同的网格大小函数  
    f = gmsh.model.mesh.field.add("MathEval")  
    if gmsh.onelab.getNumber('Parameters/Apply funny mesh size field?')[0]:  
        gmsh.model.mesh.field.setString(f, "F", "2*Sin((x+y)/5) + 3")  
    else:  
        gmsh.model.mesh.field.setString(f, "F", "4")  
    gmsh.model.mesh.field.setAsBackgroundMesh(f)  
  
    # 生成三维网格  
    gmsh.model.mesh.generate(3)  
    # 将网格写入文件  
    gmsh.write('t13.msh')  
  
# 设置用户可调的参数  
gmsh.onelab.set("""[  
  # 角度参数，用于表面检测  
  {  
    "type":"number",  
    "name":"Parameters/Angle for surface detection",  
    "values":[40],  
    "min":20,  
    "max":120,  
    "step":1  
  },  
  # 是否创建可参数化曲面的选项  
  {  
    "type":"number",  
    "name":"Parameters/Create surfaces guaranteed to be parametrizable",  
    "values":[0],  
    "choices":[0, 1]  
  },  
  # 是否应用有趣的网格大小字段的选项  
  {  
    "type":"number",  
    "name":"Parameters/Apply funny mesh size field?",  
    "values":[0],  
    "choices":[0, 1]  
  }  
]""")  
  
# 调用函数创建几何模型和网格  
createGeometryAndMesh()  
  
def checkForEvent():  
    # 检查是否有用户操作事件，如果有则重新创建几何模型和网格，并刷新图形界面  
    action = gmsh.onelab.getString("ONELAB/Action")  
    if len(action) and action[0] == "check":  
        gmsh.onelab.setString("ONELAB/Action", [""])  
        createGeometryAndMesh()  
        gmsh.graphics.draw()  
    return True  
  
# 如果命令行参数中没有'-nopopup'，则初始化图形界面，并循环检查用户事件  
if "-nopopup" not in sys.argv:  
    gmsh.fltk.initialize()  
    while gmsh.fltk.isAvailable() and checkForEvent():  
        gmsh.fltk.wait()  
  
# 清理并结束Gmsh环境  
gmsh.finalize()

14、同源和共同源计算

import gmsh  # 导入gmsh库  
import sys  # 导入sys库，用于访问与Python解释器紧密相关的变量和函数  
  
gmsh.initialize(sys.argv)  # 初始化gmsh，并传递命令行参数  
  
gmsh.model.add("t14")  # 添加一个新的模型，并命名为"t14"  
  
# 定义网格大小和几何高度  
m = 0.5  # 网格大小  
h = 2  # 几何体在z方向的高度  
  
# 添加模型的点，并指定它们的坐标和网格大小  
gmsh.model.geo.addPoint(0, 0, 0, m, 1)  
# ...（此处省略了其他点的添加，类似上面的格式）  
  
# 添加模型的线，连接之前定义的点  
gmsh.model.geo.addLine(1, 9, 1)  
# ...（此处省略了其他线的添加，类似上面的格式）  
  
# 添加曲线环和平面表面  
gmsh.model.geo.addCurveLoop([...], 13)  
gmsh.model.geo.addCurveLoop([...], 14)  
gmsh.model.geo.addPlaneSurface([13, 14], 15)  
  
# 对平面进行拉伸操作，创建一个三维体  
e = gmsh.model.geo.extrude([(2, 15)], 0, 0, h)  
  
gmsh.model.geo.synchronize()  # 同步几何模型  
  
# 添加物理组，用于后续的有限元分析  
domain_tag = e[1][1]  # 获取域的标签  
domain_physical_tag = 1001  # 定义域的物理标签  
gmsh.model.addPhysicalGroup(3, [domain_tag], domain_physical_tag, "Whole domain")  
  
# ...（此处省略了其他物理组的添加，类似上面的格式）  
  
# 添加同调请求，用于后续的网格划分和同调分析  
gmsh.model.mesh.addHomologyRequest("Homology", [...], [...], [...])  
# ...（此处省略了其他同调请求的添加，类似上面的格式）  
  
# 生成网格  
gmsh.model.mesh.generate(3)  
  
# 将模型写入文件  
gmsh.write("t14.msh")  
  
# 如果命令行参数中没有'-nopopup'，则运行Gmsh的FLTK图形界面  
if '-nopopup' not in sys.argv:  
    gmsh.fltk.run()  
  
# 清理Gmsh资源  
gmsh.finalize()

15、嵌入式点、线和面

import gmsh  # 导入gmsh库  
import sys  # 导入sys库，用于访问与Python解释器紧密相关的变量和函数  
  
gmsh.initialize()  # 初始化gmsh  
  
lc = 1e-2  # 定义网格特征长度  
# 添加四个点，并指定它们的坐标和特征长度  
gmsh.model.geo.addPoint(0, 0, 0, lc, 1)  
gmsh.model.geo.addPoint(.1, 0, 0, lc, 2)  
gmsh.model.geo.addPoint(.1, .3, 0, lc, 3)  
gmsh.model.geo.addPoint(0, .3, 0, lc, 4)  
# 添加四条线，连接之前定义的点  
gmsh.model.geo.addLine(1, 2, 1)  
gmsh.model.geo.addLine(3, 2, 2)  
gmsh.model.geo.addLine(3, 4, 3)  
gmsh.model.geo.addLine(4, 1, 4)  
# 添加一个曲线环，用于定义表面  
gmsh.model.geo.addCurveLoop([4, 1, -2, 3], 1)  
# 添加一个平面表面  
gmsh.model.geo.addPlaneSurface([1], 1)  
  
lc = lc * 4  # 更新网格特征长度  
# 设置特定线的网格大小  
gmsh.model.geo.mesh.setSize([(0, 1), (0, 2), (0, 3), (0, 4)], lc)  
  
# 添加一个新的点  
gmsh.model.geo.addPoint(0.02, 0.02, 0., lc, 5)  
  
gmsh.model.geo.synchronize()  # 同步几何模型  
  
# 将点5嵌入到表面网格中  
gmsh.model.mesh.embed(0, [5], 2, 1)  
  
# 添加更多的点和线  
gmsh.model.geo.addPoint(0.02, 0.12, 0., lc, 6)  
gmsh.model.geo.addPoint(0.04, 0.18, 0., lc, 7)  
gmsh.model.geo.addLine(6, 7, 5)  
  
gmsh.model.geo.synchronize()  # 同步几何模型  
# 将线5嵌入到网格中  
gmsh.model.mesh.embed(1, [5], 2, 1)  
  
# 对之前的平面进行拉伸操作，创建一个三维体  
gmsh.model.geo.extrude([(2, 1)], 0, 0, 0.1)  
  
# 添加一个新的点，并获取其ID  
p = gmsh.model.geo.addPoint(0.07, 0.15, 0.025, lc)  
  
gmsh.model.geo.synchronize()  # 同步几何模型  
# 将点p嵌入到网格中  
gmsh.model.mesh.embed(0, [p], 3, 1)  
  
# 添加更多的点和线  
gmsh.model.geo.addPoint(0.025, 0.15, 0.025, lc, p + 1)  
l = gmsh.model.geo.addLine(7, p + 1)  
  
gmsh.model.geo.synchronize()  # 同步几何模型  
# 将线l嵌入到网格中  
gmsh.model.mesh.embed(1, [l], 3, 1)  
  
# 添加更多的点，用于定义一个新的平面  
gmsh.model.geo.addPoint(0.02, 0.12, 0.05, lc, p + 2)  
gmsh.model.geo.addPoint(0.04, 0.12, 0.05, lc, p + 3)  
gmsh.model.geo.addPoint(0.04, 0.18, 0.05, lc, p + 4)  
gmsh.model.geo.addPoint(0.02, 0.18, 0.05, lc, p + 5)  
  
# 添加四条线，连接之前定义的点  
gmsh.model.geo.addLine(p + 2, p + 3, l + 1)  
gmsh.model.geo.addLine(p + 3, p + 4, l + 2)  
gmsh.model.geo.addLine(p + 4, p + 5, l + 3)  
gmsh.model.geo.addLine(p + 5, p + 2, l + 4)  
  
# 添加一个曲线环  
ll = gmsh.model.geo.addCurveLoop([l + 1, l + 2, l + 3, l + 4])  
# 添加一个平面表面  
s = gmsh.model.geo.addPlaneSurface([ll])  
  
gmsh.model.geo.synchronize()  # 同步几何模型  
# 将表面s嵌入到网格中  
gmsh.model.mesh.embed(2, [s], 3, 1)  
  
# 生成网格  
gmsh.model.mesh.generate(3)  
  
# 将模型写入文件  
gmsh.write("t15.msh")  
  
# 如果命令行参数中没有'-nopopup'，则运行Gmsh的FLTK图形界面  
if '-nopopup' not in sys.argv:  
    gmsh.fltk.run()  
  
# 清理Gmsh资源  
gmsh.finalize()

16、构造性立体几何，OpenCASCADE几何内核

import gmsh  # 导入Gmsh库  
import math  # 导入数学库，虽然在这段代码中未使用  
import sys   # 导入系统库，用于处理命令行参数  
  
gmsh.initialize()  # 初始化Gmsh环境  
  
gmsh.model.add("t16")  # 创建一个名为"t16"的新模型  
  
gmsh.logger.start()  # 启动日志记录  
  
# 添加两个盒子，一个边长为1，一个边长为0.5  
gmsh.model.occ.addBox(0, 0, 0, 1, 1, 1, 1)  
gmsh.model.occ.addBox(0, 0, 0, 0.5, 0.5, 0.5, 2)  
  
# 使用布尔操作"cut"从大盒子中减去小盒子  
gmsh.model.occ.cut([(3, 1)], [(3, 2)], 3)  
  
# 在剩余的大盒子顶部添加5个小球  
x = 0  
y = 0.75  
z = 0  
r = 0.09  
holes = []  
for t in range(1, 6):  
    x += 0.166  
    z += 0.166  
    gmsh.model.occ.addSphere(x, y, z, r, 3 + t)  
    holes.append((3, 3 + t))  
  
# 使用fragment操作从剩余的大盒子中减去所有小球  
ov, ovv = gmsh.model.occ.fragment([(3, 3)], holes)  
  
# 打印fragment操作产生的体积和前后关系  
print("fragment produced volumes:")  
for e in ov:  
    print(e)  
  
print("before/after fragment relations:")  
for e in zip([(3, 3)] + holes, ovv):  
    print("parent " + str(e[0]) + " -> child " + str(e[1]))  
  
gmsh.model.occ.synchronize()  # 同步几何模型  
  
# 为每个小球和剩余的体积添加物理组  
for i in range(1, 6):  
    gmsh.model.addPhysicalGroup(3, [3 + i], i)  
  
gmsh.model.addPhysicalGroup(3, [ov[-1][1]], 10)  
  
# 设置网格大小  
lcar1 = .1  
lcar2 = .0005  
lcar3 = .055  
  
gmsh.model.mesh.setSize(gmsh.model.getEntities(0), lcar1)  # 设置整个模型的网格大小  
  
gmsh.model.mesh.setSize(gmsh.model.getBoundary(holes, False, False, True),  
                        lcar3)  # 设置小球边界的网格大小  
  
eps = 1e-3  
ov = gmsh.model.getEntitiesInBoundingBox(0.5 - eps, 0.5 - eps, 0.5 - eps,  
                                         0.5 + eps, 0.5 + eps, 0.5 + eps, 0)  
gmsh.model.mesh.setSize(ov, lcar2)  # 设置特定区域的网格大小  
  
gmsh.model.mesh.generate(3)  # 生成三维网格  
  
gmsh.write("t16.msh")  # 将网格写入文件  
  
log = gmsh.logger.get()  # 获取日志记录  
print("Logger has recorded " + str(len(log)) + " lines")  
gmsh.logger.stop()  # 停止日志记录  
  
if '-nopopup' not in sys.argv:  
    gmsh.fltk.run()  # 如果命令行参数中没有'-nopopup'，则运行图形界面  
  
gmsh.finalize()  # 清理并结束Gmsh环境

网卡了
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
gmsh软件安装包
12-10
最后，"点击加入【ucbug交流反馈群】.url"、"ucbug软件下载站-www.ucbug.cc.url"和"ucbug游戏网-www.ucbug.com.url"是提供给用户交流和获取更多资源的链接，用户可以在这些平台上寻找更多关于gmsh的教程、案例分享...
gmsh.pdf最新官方版
10-08
软件还提供了详细的参考手册，包括如何阅读手册的指南、系统运行说明、各模块详细介绍、教程和选项列表等，确保用户能够充分理解和使用Gmsh的各项功能。 总的来说，Gmsh是一个强大而全面的工具，它集成了几何建模、...
Gmsh概述
weixin_44897632的博客
07-15 1788
Gmsh是一个内置CAD引擎和后处理器的三维有限元网格生成器。其设计目标是提供一个快速、轻量级且用户友好的网格划分工具，具有参数化输入和灵活的可视化功能。Gmsh围绕四个模块构建（几何、网格、求解器和后处理），这些模块可以通过图形用户界面（GUI；见Gmsh图形用户界面）、命令行界面（见Gmsh命令行界面）、使用Gmsh自己的脚本语言编写的文本文件（.geo文件；见Gmsh脚本语言）或通过C++、C、Python、Julia和Fortran应用程序编程接口（API；见Gmsh应用程序编程接口）进行控制。
Gmsh剖一维网格教程附代码
guihunkun的博客
09-20 2476
文章目录前言一、点，线段，单元，多边形等的定义1. 点的定义2. 多边形的定义3. 线段的定义4. 拓扑关系的定义5. 单元的定义二、区域构造三、调用Gmsh的API剖网格并得到网格信息四、程序运行结果1. 程序输出2. 带有点标号的结果图3. 带有单元标号的结果图五、总结 前言 Gmsh的安装与入门教程可以去看我的前一篇文章有限元剖网格之Gmsh安装与使用入门。本篇文章主要介绍如何借助Gmsh的C++ API剖一维网格(即线段)。 一、点，线段，单元，多边形等的定义 因为我们需要自己构造区域内的几何信
Gmsh剖三维网格教程附代码
guihunkun的博客
01-23 3245
文章目录前言一、程序代码二、程序结果1. 带有点标号的结果图2. 没有点标号的结果图三、总结 前言 前面写过三篇关于Gmsh的教程： 有限元剖网格之Gmsh安装与使用入门 Gmsh剖一维网格教程附代码 Gmsh剖二维网格教程附代码 以上文章对Gmsh都有一定的介绍，在此不再重复直接接剖三维网格的C++代码。 注：因有一个OpenCASCADE库暂时没安装好，所以目前只能用插点，连线，线连面，面拼体的方式构造三维几何体。 一、程序代码 #include <iostream> #include
gmsh边界条件
novanova2009的博客
03-02 563
How do I define boundary conditions or material properties in Gmsh? By design, Gmsh does not try to incorporate every possible definition of boundary conditions or material properties—this is a job best left to the solver. Instead, Gmsh provides a simple m
gmsh-4.6.0-Windows64_GMesh_gmsh_使用和开发_
10-03
gmsh的教程则更注重交互式建模和脚本编写，用户可以通过实例学习如何创建复杂的几何模型并生成网格。 在开发层面，两者都有开源的特性，这意味着开发者可以查看源代码，学习其内部工作原理，甚至为软件贡献新的功能...
Gmsh-4.0.6-source.zip_C++_FEM_Fortran_gmsh 有限元_gmsh-4.5
07-14
手册通常包含了软件的安装指南、基本操作教程、API文档以及一些示例，帮助用户快速上手并深入掌握Gmsh的各项功能。 在有限元分析中，前处理通常涉及几何建模和网格划分，后处理则涉及结果的可视化和分析。Gmsh不仅...
pygmsh:适用于Python的Gmsh
04-28
适用于Python的Gmsh。 pygmsh将的功能与Python的多功能性结合在一起。 它提供了Gmsh自己的Python界面的有用抽象，因此您可以更轻松地创建复杂的几何图形。 要使用，请从安装Gmsh本身和pygmsh： [sudo] apt install ...
cv2读取和保存图片
m0_53291740的博客
07-12 269
cv2.imwrite('b\\img1.png',img)#b存在就可以保存进去。如果b存在，c不存在，程序不会报错，也不会保存图片。如果目录b不存在就会存到当前文件夹。
《JavaSE》---6.＜基础语法（Java三大程序控制结构）＞
m0_73456341的博客
07-14 787
本篇博客主要讲解Java基础语法中的三大结构，一种顺序结构、两大分支结构i（if-else、swich-case）、四大循环结构（while、do while、fot、foreach）
华为od机试真题 — 代表团坐车(Python)
学习,你不是一个人在战斗 ～ 961302305（QQ,微信）
07-17 527
【华为od机试真题 】代表团坐车(Python) 某组织举行会议，来了多个代表团同时到达，接待处只有一辆汽车可以同时接待多个代表团，为了提高车辆利用率，请帮接待员计算可以坐满车的接待方案输出方案数量。
IPython 日志秘籍：%logstate 命令全解析
2401_85812026的博客
07-17 490
logstate命令是 IPython 日志系统的重要组成部分，它帮助我们掌控日志记录的状态。通过本文的介绍，你应该对如何在 IPython 中使用%logstate查看日志状态有了更深入的理解。记住，合理利用日志系统，可以让你的数据探索过程更加透明和可追溯。本文详细介绍了 IPython 的%logstate命令，包括其基本概念、基本用法、输出解析以及实战示例。希望本文能够帮助读者更好地利用 IPython 的日志系统，优化数据探索和科学计算的工作流。记住，掌握%logstate。
Python每日一题：回文数
xiongxiaoxuan的博客
07-16 235
回文数是指正序（从左向右）和倒序（从右向左）读都是一样的整数。例如，121 是回文，而 123 不是。给你一个整数 x ，如果 x 是一个回文整数，返回 true；否则，返回 false。解释：从左向右读, 为 -121。从右向左读, 为 121-。因此它不是一个回文数。先将整数x转化为字符串形式，判断x和它的反转形态是否完全相等。解释：从右向左读, 为 01。需注意：True和False首字母大写！提示：-231
N7翻译实战
最新发布
tjl521314_21的博客
07-17 225
本周完成项目实战用于训练一个简单的序列到序列（seq2seq）模型以实现英语到法语的翻译。数据预处理、模型构建、训练以及可视化损失的过程。本周学习了构建和训练一个简单的seq2seq模型用于英语到法语的翻译，对前面的知识做了一个综合的运用。
Gradio从入门到精通（1）---快速入门
疯狂飙车的蜗牛的博客
07-15 1360
Gradio 是一个开源 Python 包，允许您为机器学习模型、API 或任何任意 Python 函数快速构建演示或 Web 应用程序。然后，您可以使用 Gradio 的内置共享功能在几秒钟内共享指向演示或 Web 应用程序的链接。无需 JavaScript、CSS 或 Web 托管经验！Gradio最大的特点就是简单，便捷，短短几行代码就可以快速构建大模型web应用程序；当前很多人工智能模型的演示都是使用Gradio打造的；您会注意到，为了制作第一个演示，您创建了该类的实例。
sklearn基础教程：掌握机器学习入门的钥匙
AIGC绘图领域探索者，热衷于将创意与技术融合，致力于分享前沿AI绘图技巧与心得，让我们一起开启视觉艺术的新篇章。
07-13 970
是基于Python的一个开源机器学习库，它建立在NumPy、SciPy和matplotlib之上，提供了大量的算法和工具，用于数据挖掘和数据分析。无论是简单的线性回归，还是复杂的神经网络，sklearn都能提供直观易用的接口。作为Python中最为流行的机器学习库之一，以其简洁的API、丰富的算法和强大的功能赢得了广大用户的青睐。通过本文的介绍，相信你已经对sklearn有了初步的了解，并能够开始进行一些基础的机器学习任务。
文件在Linux下为binary格式在Windows下为utf-8格式
qq_38220334的博客
07-17 267
想着是不是之前遇到的问题，有可能是文件中含有特殊字符，然后用vim命令来编辑这个文件，发现，还真是含有一个特殊字符，然后将它删除掉，文件就变成utf-8格式的文件了。在Windows下和Mac下都是没有问题的，在华为云服务器上也是没问题的（操作系统不是centos），但是在自己搭建的Linux服务器（centos）上是有问题的。用shell或者Python脚本去将这个文件转为utf-8格式的，但是都失败了。
python gmsh
11-20
根据提供的引用内容，可以看出gmsh是一个用于生成三维有限元网格的软件，而Python可以通过gmsh提供的API来调用gmsh的功能。如果需要长期使用gmsh，建议学习使用Python编写，因为gmsh提供的API功能很全，使用Python更方便集成到其他项目中，能够使自动化程度更高。 安装gmsh可以使用pip进行安装，命令为：`pip install --upgrade gmsh`。安装完成后，可以在Python中使用gmsh的API来生成三维有限元网格。 具体使用方法可以参考gmsh官方文档，文档中提供了详细的API说明和示例代码。以下是一个简单的示例代码，用于生成一个球体的三维有限元网格： ```python import gmsh # 初始化gmsh gmsh.initialize() # 创建一个新模型 gmsh.model.add("sphere") # 定义一个球体 lc = 0.1 gmsh.model.geo.addSphere(0, 0, 0, 1, meshSize=lc) # 生成网格 gmsh.model.mesh.generate(3) # 将网格写入文件 gmsh.write("sphere.msh") # 退出gmsh gmsh.finalize() ```

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值