对于传统的有限元法和有限元法,为什么不能用网格来表示几何,而直接用表示几何的网格来进行计算呢?

计算科学 计算几何 网格生成
2021-12-22 17:00:30

等几何分析 [ 1 ] 具有使用 NURBS 或 Spline 集成几何模型和网格模型的优势。同时我想问朋友一个问题:对于传统的FEM和FVM,为什么不能用mesh来表示几何,直接用表示几何的mesh来做计算呢?

具体来说,根据离散微分几何,如果能找到连续几何与离散几何的等价定义,就可以用网格来表示连续几何。我们知道,对于CG软件(例如Blender),建模的主要方式是使用网格来表示连续的几何体,我们通过修改网格来修改几何体。但是对于科学计算网格生成软件(例如Gmsh),我们通常需要先建立几何模型,然后设置一些参数来进行网格划分。他们之间有很大的不同。我知道 CG 主要使用 3D 空间中的表面网格。科学计算主要使用实体网格。但我认为CG使用网格来表示几何的方式更简单,适合形状变化(优化)。所以我想知道为什么不能 我们在科学计算中使用网格来表示像 CG 这样的几何图形。几何图形直接由网格表示,然后您可以直接在其上进行计算。我认为对于 FEM 和 FVM 的先驱者来说,将几何和网格结合的想法之前一定有考虑过。但是现在几乎所有的仿真软件都没有结合几何和网格模型。所以我想知道不这样做的主要困难是什么。所以想问一下朋友们:对于传统的FEM和FVM,为什么不能用mesh来表示几何,直接用表示几何的mesh来做计算呢?之前必须考虑过结合几何和网格的想法。但是现在几乎所有的仿真软件都没有结合几何和网格模型。所以我想知道不这样做的主要困难是什么。所以想问一下朋友们:对于传统的FEM和FVM,为什么不能用mesh来表示几何,直接用表示几何的mesh来做计算呢?之前必须考虑过结合几何和网格的想法。但是现在几乎所有的仿真软件都没有结合几何和网格模型。所以我想知道不这样做的主要困难是什么。所以想问一下朋友们:对于传统的FEM和FVM,为什么不能用mesh来表示几何,直接用表示几何的mesh来做计算呢?

这个问题的灵感来自我上一个问题的答案:FVM和FEM中的网格生成方法是否可以完全基于计算机图形学中的网格生成理论知识?

参考书目

1 . 休斯,TJR;科特雷尔,JA;Bazilevs, Y.,等几何分析:CAD,有限元,NURBS,精确几何和网格细化,计算机。方法应用程序。机甲。英。194, No. 39-41, 4135-4195 (2005)。ZBL1151.74419

3个回答

潜在的问题是网格实际上不是几何体。你想模拟一座桥?它具有一定的几何形状,您可以使用网格进行近似,但网格不是精确的几何形状。大多数其他弯曲物体也是如此。

当然有一些方法可以将几何图形集成到有限元网格中。特别是,我想向您指出这个预印本主要障碍是历史性的:传统的工作流程是使用 CAD 软件描述精确的几何图形;该 CAD 几何图形被馈送到创建网格的网格生成器;然后将网格提供给有限元(“分析”)软件,但它不再可以访问 CAD 几何。链接到的预印本详细描述了将 CAD 几何图形一直路由到分析引擎需要做什么,以及为什么这样做很有用。

“所以我想知道为什么我们不能在科学计算中使用网格来表示像CG这样的几何体。几何体直接用网格表示,然后你可以直接在上面进行计算。” ...

这正是我们在 FEM 中使用等参数元素所做的;我们使用相同的网格离散几何和场变量。等参元素的概念不限于拉格朗日族元素或 NURBS。理论上,可以使用任何适当的多项式空间来离散几何和场变量,例如细分曲面、T-Splines、Box-Splines、Chebyshev 多项式等。

在 FVM 中,几何变量和场变量也使用相同的网格。但是等参数表示的概念不适用于 FVM,尤其是对于以细胞为中心的 FVM 方法。如果我错了,请纠正我。

乍一看,IGA 听起来很有希望。与传统的 FEM 相比,它当然有它的优势。IGA 已经在 LS-DYNA 中可用。但是,我认为在实际应用中,IGA 目前的劣势超过了它的优势。还有其他非技术性的问题,见第 5 点。

1.) IGA 并非不受空间离散化的影响。仍然必须离散几何体(如果这有意义的话,将其与“网格”一词区分开来)。就 IGA 而言,需要为 NURBS 生成一个新的节点空间。CAD 软件中尚不提供此类功能。它主要是手动执行或在学术研究小组中使用自定义脚本执行。

2.) IGA 植根于跨元素边界的更高连续性的想法。虽然在有限元法中较高的连续性对于某些特定问题是有利的,特别是对于具有平滑解的问题,但这一点对于场变量不连续的问题并不是那么有利。此外,更高的连续性仅限于没有尖角的单个贴片。

此外,元素间的不连续性是以细胞为中心的 FVM 所固有的。因此,我认为直接将 IGA 用于 FVM 是没有意义的。(我在 FVM 方面没有太多经验来进一步详细说明)。

3.) 与传统 FEM 相比,IGA 成本高。更高的连续性意味着(有效)刚度矩阵中有更多的非零条目。虽然我们确实可以使用具有高阶元素的粗网格获得准确的结果,但在实际应用中使用高于二次多项式的效果并不明显。根据我的经验,在考虑准确性和效率(运行时)时,二次多项式是最佳选择。

4.) 将 IGA 扩展到固体力学中的高级问题也很困难,例如塑性和不可压缩超弹性,这需要复杂的公式。

  • 基于投影的方法,例如 F-bar(或有人称之为 B-bar)公式,需要对块状质量矩阵进行反转。-> 昂贵。
  • Taylor-Hood 类型单元(压力场的多项式比位移(或速度)场的多项式低一阶)不是 inf-sup 稳定的。需要使用细分稳定来确保 inf-sup 稳定性。(类似于 Q1-iso-Q1 元素)。这有实际的局限性。

如果我们在使用 IGA 时不必求助于复杂的配方,那就太棒了。但是,不幸的是,我们仍然必须使用那些先进的配方。

5.) 基于使用传统 FEM 和 FVM 的软件工具,业界的设计和分析工作流程已经很成熟。除非我们向他们展示一些快速的回报或显着的节省,否则很难说服行业改变他们的工作流程以使用新的模拟范例。

可以肯定地说,IGA 还不够成熟,无法取代传统的 FEM。很难说什么时候会。但是可以结合 IGA 和传统 FEM 的概念和工具来改进现有方案,参见论文1论文 2

是的,您可以使用相同的网格来表示域的几何形状并求解 PDE,这是您可以做到的。例如,一个正方形可以完全由两个三角形描述,您可以使用此网格来求解您的 PDE。如果您的 PDE 的解是恒定的或线性的,则此网格就足够了,但如果您的解中有更高的梯度,则需要更多的元素。

通常,会考虑一系列网格,并且有一个足够精细的网格(适用于您的应用程序),它会在您所需的容差下给您一个错误。例如,求解

2u=80sin(4x)cos(8y),

为了(x,y)[0,π]2需要一个正方形的网格。然而,仅使用两个元素不会为您提供接近分析结果的解决方案u=sin(4x)cos(8y).