我正在尝试使用雅可比行列式作为度量来确定我的网格元素的质量。我的算法采用三角形网格中的节点顶点并将它们移动以形成扭曲的三角形。让我们假设我的顶点 P(x,y,z) 移动到 P'(x',y',z') 如何计算网格质量的最终变化?一种方法是计算网格中每个三角形的雅可比,但我不明白如何计算仅给定节点位移的偏导数。
所有关于雅可比计算的现有文献似乎都使用一个函数来做到这一点。谁能简单的解释一下?
我正在尝试使用雅可比行列式作为度量来确定我的网格元素的质量。我的算法采用三角形网格中的节点顶点并将它们移动以形成扭曲的三角形。让我们假设我的顶点 P(x,y,z) 移动到 P'(x',y',z') 如何计算网格质量的最终变化?一种方法是计算网格中每个三角形的雅可比,但我不明白如何计算仅给定节点位移的偏导数。
所有关于雅可比计算的现有文献似乎都使用一个函数来做到这一点。谁能简单的解释一下?
您可以为网格质量定义多个度量。例如,gmsh使用衡量元素形状的标准。它基于元素半径比,即内接圆与外接圆之间的比率(参见参考文献 1)。我认为您也可以使用雅可比来描述元素的质量,尽管您可能需要查看参考 2。
关于雅可比,首先要提到的是它是高阶元素中元素坐标的函数。然后,对于这些元素,您需要计算每个元素的优值,例如,行列式的平均值或规范的平均值。对于线性元素,您将进行如下图所示的变换
对于这个简单的情况,转换由
其中
和。这是考虑转换的最简单方法,您可以发现它解决了线性方程组,它仅取决于节点的位置。
对于您的问题,您想计算雅可比行列式,然后将位移应用于每个节点,然后再计算一次。
参考