当使用Hexpress网格生成器时,可以获得六面体单元格,可能带有悬挂节点。由于平滑步骤,悬挂节点可以移动:它们不一定在连接角的直线上,如下所示。

使用 FVM,可以用 5 个节点对面进行三角剖分,以近似其表面、面中心等,例如参见 Ferziger 和 Peric 所著的《流体动力学计算方法》,Springer,2001 年。如何在 FEM 中处理这种情况? 我熟悉使用线性形状函数将 hexa 8 映射到主元素,但这似乎还不够。
当使用Hexpress网格生成器时,可以获得六面体单元格,可能带有悬挂节点。由于平滑步骤,悬挂节点可以移动:它们不一定在连接角的直线上,如下所示。

使用 FVM,可以用 5 个节点对面进行三角剖分,以近似其表面、面中心等,例如参见 Ferziger 和 Peric 所著的《流体动力学计算方法》,Springer,2001 年。如何在 FEM 中处理这种情况? 我熟悉使用线性形状函数将 hexa 8 映射到主元素,但这似乎还不够。
最好的方法是考虑为额外节点添加等参数元素。例如,参见
http://www.edwilson.org/book-wilson/05-iso.pdf
第 5.5 节对 2D 双线性元素进行了讨论。3D 三线性的扩展相当简单。
请注意,这需要代码在局部元素空间方面具有一定的灵活性——这个元素有 9 个元素而不是 8 个。关于这样做的数学和这样做的软件设计的最佳参考是这本书:
http://www.amazon.com/Finite-Element-Method-Mechanical-Engineering/dp/0486411818