openfoam中为什么一个单元用27个点表示，代表什么图形（由27个节点组成的三维立方体单元，在有限元方法（FEM）中被称为“三次立方体单元”）

问题：

近期在做openfoam项目的时候，发现openfoam中固体的点为什么一个单元用27个点表示，想着代表什么图形呢？如果以顶点表示的话好像图形就复杂了，然后查询一下资料，结果如下

解答：

在OpenFOAM中，使用27个点来表示一个单元通常指的是一种高阶单元。这种单元类型在有限元方法（FEM）中被称为“三次拉格朗日立方体单元”或“三次立方体单元”。它是由27个节点组成的三维立方体单元，用于更精确地表示场变量，如位移、应力和应变。

上图中展示了一个27节点三次拉格朗日立方体单元的节点分布。这种高阶单元可以更准确地模拟物理现象，特别是在固体力学和结构分析中。

27个节点的三次立方体单元

这种单元有以下特点：

• 8个顶点节点：位于立方体的每个顶点。
• 12个边节点：位于立方体的每条边的中点。
• 6个面节点：位于立方体的每个面的中心。
• 1个中心节点：位于立方体的中心。

节点编号如下：

1. 顶点节点（8个）：

   • 底面顶点：0, 1, 2, 3
   • 顶面顶点：4, 5, 6, 7

2. 边节点（12个）：

   • 底面边节点：8, 9, 10, 11
   • 顶面边节点：12, 13, 14, 15
   • 垂直边节点：16, 17, 18, 19

3. 面节点（6个）：

   • 底面中心：20
   • 顶面中心：21
   • 四个侧面中心：22, 23, 24, 25

4. 中心节点（1个）：

   • 立方体中心：26

优点

使用27节点三次立方体单元的主要优点是：

1. 更高的精度：由于这种单元具有更多的节点，可以更精确地表示复杂的几何形状和物理场。
2. 更好的曲面拟合：三次形函数能够更好地拟合曲面，适用于非线性问题。
3. 减少网格依赖性：相比于低阶单元（如线性单元或二次单元），高阶单元可以在较少的网格情况下提供更好的精度。

应用

在固体力学模拟中，27节点单元适用于高精度应力分析和复杂几何结构的模拟。例如，在需要精确捕捉应力集中区域或非线性行为的场景中，这种高阶单元能够提供更好的模拟结果。