用来量化三角形网格不规则性的常用指标


9

假设您在平面上有一个三角形网格。例如,这是为了最终解决力学中的某些问题而绘制的。

由于顶点之间和质心之间的距离始终相同,所以等边三角形的网格是最好的。这使得插值和梯度计算变得简单而准确。但是,由于种种限制和情况,并非总是能够对所有等边三角形的网格进行加工。

因此,问题涉及任意形状的三角形元素的网格。

关于单个网格元素。通常使用哪些度量来量化一个通用三角形与某些潜在的理想等边形状的相异性?

关于整个网格。哪些指标用于量化整体上任意三角形网格的不规则性?这些指标应指示网格的混乱程度。

感谢您的配合。

注意 来自有限元社区的所有贡献都得到了极大的赞赏。对于此问题,请注意,仅对几何形状(任意三角形与等边三角形)中的差异进行量化是很有意义的。对插值和条件误差的后续影响不在范围内。鉴于这些可以洞察力和相关性,它们使数学处理变得复杂。


4
你检查过这个问题了吗?从那篇文章开始:“什么是好的有限元?”。
nicoguaro

3
我认为,外接圆和外接圆的面积/半径之比可能有效。雅可比角,最小角和最大角的特征值之比也是如此。
nicoguaro

4
Shewchuck最著名的论文之一深入探讨了这个话题: 什么是好的线性有限元?
保罗

@nicoguaro谢谢。我对FEM并不是特别感兴趣,但是对量化元素形状的差异非常感兴趣。例如,您能否详细说明半径比?这与大小无关吗?换句话说,如果您可以在答案中列出您的选项供其他人继续使用,将不胜感激。
XavierStuvw

您也可以查看任何网格元素中的最小角度。这个想法是要尽可能的大
-KyleW

Answers:


4

就像@Nicoguaro和@Paul在问题贴的评论中说的那样,有很多方法可以做这种事情,而且我不确定是否有单一的“最佳”方法。


从伯克利的乔纳森·理查德·谢克Jonathan Richard Shewchuck)的评论研究中,答案是:

在此处输入图片说明

请参考原始文档(2002年12月31日版)以获取符号,术语,特殊功能以及更多内容(例如,四面体)。第6章是关于质量措施。链接到的文档是扩展版本,在JRS的网页中也有一个简短的文档。


就个人而言,我是“卷长”指标的粉丝。它是(各向同性)单纯形质量的很好的鲁棒标量指标,并且计算便宜。二维:

一个=433一个Ë[Rs2

哪里 一个是三角形和的有符号区域 是均方根边缘长度。理想元素达到,随着失真增加,该值朝零减小。方向相反的反向元素。Ë[Rs一个=1个一个<0

为了评估非结构化三角剖分的质量,通常需要查看此类元素质量指标的直方图。这类事情的实现有很多,但是我的一个简单MATLAB代码库就在这里

除体长分数外,默认情况下还计算元素角和顶点度的直方图。


您为什么喜欢这个指标?它能很好地预测您对网格所做的模拟的准确性吗?
BrunoLevy

@BrunoLevy:嗯,作为单纯形的简单“默认”选择:它可以可靠地推广到高维,计算便宜,在数值上条件良好,提供了一个“纠结”指标。方向,并且根据问题是一个简单的“仅几何图形”指示器。它是模拟质量的良好指标吗?好吧,这取决于您在做什么!如果您对各向同性网格物体感兴趣,我会说是的。高度依赖于方向的各向异性配置,然后不是,不是直接配置,尽管在这种情况下,仍可以在适当的坐标转换后使用它。
Darren Engwirda

1
也是平滑的,因此可以使用隐式任意拉格朗日-欧拉公式来运行它。您可以轻松地将其推广到各向异性网格。
likask

@likask:是的,很好-它可以是网格平滑和优化的良好成本函数。
Darren Engwirda '17

我添加了一份有关Shewcuck著作的摘录,该摘录扩展了Darren的回答范围。这也总结了几个评论。感谢这篇文章的所有贡献者。
XavierStuvw'7

4

我不认为这个问题总体上存在答案,因为这完全取决于网格的预期用途。例如,如果您要进行计算流体动力学,则可能需要在边界层附近有一个非常各向异性的网格。现在,如果您正在进行计算电磁学,则最佳网格可能会完全不同。

文献中有许多关于“网眼质量”标准的不同定义。他们中的大多数会喜欢带有等边三角形的网格。还可以提到最大化最小角度的想法(这是通过Delaunay三角剖分实现的,用于固定点集)。乔纳森·谢丘克(Jonathan Shewchuk)在评论中提到的分析证明了这一点的正确性,该分析将该角度与用P1元素离散化的拉普拉斯方程的刚度矩阵的条件数相关联,但同样,根据预期的用途,某人的良好网格可以是某人否则网眼不好。

我认为“纯粹量化几何差异(任意三角形与等边三角形)”是没有道理的:在测量三角形是否等边并确定哪种“偏差与等边性”是最好的三角形之前,有必要弄清楚我们是否想要“等边三角形”,并非总是如此!所有这些都来自您提到的“插值和条件”。是的,正如您所说的“它使数学处理变得复杂”,但如果没有它,就不可能在给定应用程序的客观标准与根本没有意义的标准之间有所区别。

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.