岩土数值计算的网格剖分一直是一个大难题。我们团队前面弄了一个MetaMesh程序,当时理想比较宏大,但是限于时间和精力,难以维护这么庞大的C++程序。所以期间我也没有闲着,沿C++代码的思路写了一个Python版本的网格生成工具,暂且叫做MeshWorks吧。
这几年感谢很多甲方爸爸对我们的支持,算是稍微手头有了一点点经费。但是基本上甲方给我的项目经费都养团队了,我自己现在开的是五菱,真的很穷。
我们都知道平面应变问题分析是一个比较危险的极端工况。如果最危险剖面的平面应变结果算过了,三维基本上问题也不大,所以很多时候二维计算还是有需要的。二维计算一般剖面要算的比较多,有时候地质体(断层、节理、透镜体)会导致就算是二维也比较繁琐,有时候一个模型要几个小时甚至两三天,为了解决这个痛点,我们开发了一个任意二维/假三维网格的自由剖分的功能。
为了最大限度的方便大家,我们做了一个一键生成的策略。只需要把一系列的线条准备好就可以了,然后点击一下按钮可以了。一个示例如下:
(1) 输入曲线清理容差
(2)框选所有的线条
(2)调整样条曲线的简化精度,根据研究结构的大小,对于岩土工程,一般最小选1,最大选2就差不多了,个人感觉不一定准。
(3)网格种子密度调整,为了让网格兼顾精度和控制数量,我们给外边缘设置最大尺寸,内部结构设置最小尺寸,这样网格看起来也会舒服一点。
(4)等待程序自动完成网格剖分
这是被多个生产项目验证过的程序,我们也随手剖分了一些案例给大家看看效果。我们的程序也支持通过设置容差过滤掉很小的边,防止出现那种畸形的小边。程序支持导出为ABAQUS的inp、FLAC的f2grid(2D格式)和f3grid(3D格式)。欢迎大家应用我们的工具解决具体的工程或者学术问题!
我们还有同款的三维任意模型的一键生成功能,堪称做项目神器,敬请期待!