CAE底层网格算法的无穷美学【附专用网格处理工具】

7天前浏览804

网格划分技术作为有限元仿真中的核心一环，历来是工程师们头疼且费时较多的一个环节，目前主流商用软件都提供通用的网格划分功能。作为软件用户群体，我们很少去探究网格划分背后的算法问题，而且也无需去做过多的探索，不过今天小编将截取部分论文研究成果，为大家展示从底层计算机图形学方面提供的算法实践成果，这些成果涵盖了三角形、四边形、四面体、六面体，或是粒子群等离散形式，同时为大家推荐几个网格专用处理工具，感兴趣的可以作进一步的深入了解和研究，在此特别膜拜一下国内在计算机图形方面的大牛团队，浙大鲍虎军、山大陈宝权、MSRA团队以及国外亚琛工业大学Prof. Dr. Leif Kobbelt团队，相关图片均来自上述团队论文成果。

所以从层级上讲，网格划分是有限元开展的基础工作，而这些底层图形学算法又是网格划分的基础，也就是基础中的基础，咱们是用户层，人家是软件开发层，还是人家更厉害，不得不服。

1 部分论文成果图

下面是一些看起来极度舒适的离散图片，对于强迫症患者而言更是有一种心灵上的解脱，一起观赏一下（由于作者水平和专业所限，尚无法为大家一一解读论文，感兴趣的可以自行检索，图文附引用）。

[1] Frame Field Generation through Metric Customization. Tengfei Jiang, Xianzhong Fang, Jin Huang, Hujun Bao, Yiying Tong, Mathieu Desbrun ACM Transaction on Graphics(TOG),(conditional accept),2015.

[2] Perceptually Guided Rendering of Texture Point-based Models. Lijun Qu, Xiaoru Yuan, Minh X. Nguyen, Gary Meyer and Baoquan Chen In Proceeding of the 3rd IEEE/Eurographics Symposium on Point-Based Graphics (PBG06), Pages 95-102. Boston, MA, USA. July 29-30, 2006.

[3] All-Hex Meshing using Closed-Form Induced Polycube. Xianzhong Fang, Weiwei Xu, Hujun Bao, Jin Huang，ACM Transactions on Graphics (TOG), 35(4), 2016.

[4] Integrating Mesh and Meshfree Methods for Physics-Based Fracture and Debris Cloud Simulation. Nan Zhang, Xiangmin Zhou, Desong Sha, Xiaoru Yuan, Kumar K. Tamma, and Baoquan Chen. In Proceeding of the 3rd IEEE/Eurographics Symposium on Point-Based Graphics (PBG06), Pages 145-154. Boston, MA, USA. July 29-30, 2006.

[5] A Convolutional Decoder for Point Clouds using Adaptive Instance Normalization. Isaak Lim, Moritz Ibing, Leif Kobbelt.Eurographics Symposium on Geometry Processing 2019.

[6] All-Hex Meshing using Singularity-Restricted Field. Yufei Li,Yang Liu,Weiwei Xu,Wenping Wang,Baining Guo. ACM Transactions on Graphics (TOG) - Proceedings of ACM SIGGRAPH Asia 2012.

[7] Structured Volume Decomposition via Generalized Sweeping.Gao X, Martin T, Deng S, Cohen E, Deng Z, Chen G.IEEE Trans Vis Comput Graph. 2016 Jul;22(7):1899-911. doi: 10.1109/TVCG.2015.2473835. Epub 2015 Aug 27.

[8] Surface Smoothing and Quality Improvement of Quadrilateral/Hexahedral Meshes with Geometric Flow. Zhang Y, Bajaj C, Xu G. Computational Visualization Center, Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA. jessica@ices.utexas.edu.

2 软件工具推荐

最后，推荐几个轻型的专用网格处理软件（非Hypermesh、ANSA），如果自己研究底层程序编写或者自定义一些模型求解过程，可以作为网格生成工具使用。

2.1 MsehLab

MeshLab是一个开源的处理三角形网格的C++处理框架，提供了三角网格的数据结构和算法，诸如曲面重建、编辑、修复、光顺、编辑等算法。MeshLab并没有集成太多独特的算法，但是作为一个三维网格数据的显示工具和框架已被学术界广泛使用，并作为科研的必备程序库之一。

2.2 Libigl

Libigl是由瑞士Ethz大学的Olga Sorkine研究小组开发的C++网格处理库，使用说明文档比较齐全，很容易上手使用。它具有较好的构造稀疏离散微分算子和有限元稀疏方程组等功能。有很多研究人员对其有贡献，包括：Alec Jacobson, Daniele Panozzo, Christian Schüller, Olga Diamanti, Qingnan Zhou, Nico Pietroni, Stefan Bruggerr , Kenshi Takayama, Wenzel Jakob, Nikolas De Giorgis, Luigi Rocca, Leonardo Sacht, Olga Sorkine-Hornung。

2.3 Trimmesh

Trimesh是由美国Princeton大学的Thomas Funkhouser研究小组开发的C++网格处理库。非常容易上手，使用不难，适合初学者。但是其中的相关算法实现比较少。

2.4 OpenMesh/OpenFlipper

OpenMesh是由德国RWTH Aachen大学的Leif Kobbelt研究小组开发的C++网格处理库。OpenFlipper是基于OpenMesh基础上架构的网格处理框架，使用非常广泛。

2.5 TetGen

(A Quality Tetrahedral Mesh Generator and a 3D Delaunay Triangulator)：是最有名的空间四面体网格生成库，由华人学者Hang Si博士所开发。该算法库获得了2012年SGP会议的最佳开源软件奖。

2.6 CGAL (Computational Geometry Algorithms Library)

CGAL是一套开源的C++算法库，提供了计算几何相关的数据结构和算法，诸如三角剖分（2D约束三角剖分及二维和三维Delaunay三角剖分），Voronoi图（二维和三维的点，2D加权Voronoi图，分割Voronoi图等），多边形，多面体（布尔运算），网格生成（二维Delaunay网格生成和三维表面和体积网格生成等），几何处理（表面网格简化，细分和参数化等），凸壳算法，搜索结构（近邻搜索，kd树等），插值，形状分析，拟合等。

2.7 VEGA FEM

Vega is a computationally efficient and stable C/C++ physics library for three-dimensional deformable object simulation. It is designed to model large deformations, including geometric and material nonlinearities, and can also efficiently simulate linear systems. Vega contains about 145,000 lines of code, and is open-source and free. It is released under the 3-clause BSD license, which means that it can be used freely both in academic research and in commercial applications.

以上，是关于有限元程序开发以及使用过程中可能涉及的底层离散算法。这些算法充满艺术美学，让人眼前一亮，有门槛又非常脑洞有趣，希望对大家产生一些帮助，拓展一些技术边界，如果觉得文章内容不错，欢迎“在看”和分享~

特别致谢:

a. 中国科学技术大学刘利刚教授(Prof. Dr. Ligang Liu) || b. CSDN会员hjwang1

来源：纵横CAE

网格划分，画个“球”啊！

球体零件在工程中应用广泛，如深沟球轴承、滚珠丝杠等，CAE仿真时经常会进行球类结构的网格划分。该类结构看似简单，却蕴含着高质量六面体网格单元划分技巧。滚珠丝杠中的“滚珠”深沟球轴承的“滚动体”本文联合运用Solidworks和Workbench，详细介绍球类结构高质量网格划分方法。难度不大，技巧满满。看文本文，你就会了，尤其是对于初学人员会有很大的帮助。球类结构高质量六面体网格实体单元在SolidWorks环境下，在球体中心区域切除一个边长为球径一半的正方体，然后再拉伸一个正方体不合并，最后利用基准平面、分割等功能进行切分，如下图所示。球体结构智能切分完成球体智能切分后，点击SolidWorks菜单栏中的“工具“，展开Ansys，选择Ansys Workbench，将模型传输至Ansys Workbench的Germetry中，如下图所示。球体结构模型传输进入Workbench界面后，设置Units为mm，右击Geomrtry，选择Edit Geometry in DesignModeler....，进入DM界面，右击Attacxh1，选择Generate生成几何模型。导入并生成球体结构关闭DM界面，拖拉Component Systems中的Mesh至Geometry上，创建球体结构网格划分流程，如下图所示。点击File，选择Save，命名并保存分析项目。球体结构网格划分流程双击Mesh，进入Meshing界面，采用自动网格划分，右击Mesh，选择Generate Mesh生成网格，如下图所示。当然可以根据实际计算需求，设置局部网格控制精细化网格模型。球体网格模型右击Model，插入两个Named Selection，图形区分别选择1/2球体模型和1/8球体模型，点击下方列表中的Apply，便于详细查看球体结构内部网格模型，如下图所示。1/8球体网格模型1/2球体网格模型点击Mesh，展开下方列表中的Quality，Mesh Metric中进行选择网格质量评价指标，如Element Quality。具体参见文章：一文搞懂Ansys Workbench网格质量评价。球体结构网格质量评价虽然球体拓扑结构规则、模型对称简单，但是生成高质量的六面体网格实体单元，还是需要用到一些网格划分技巧。对于初学者来说，是一个很好的练习题目。喜欢作者，请点赞和在看来源：纵横CAE