如何在 3D 4 节点元素上集成多项式表达式?

计算科学 有限元 正交
2021-11-27 02:49:51

我想在 3D 的 4 节点元素上集成多项式表达式。几本关于 FEA 的书籍涵盖了在任意平面 4-noned 元素上执行积分的情况。在这种情况下,通常的程序是找到 Jacobi 矩阵并使用它的行列式将积分基础更改为归一化的,其中我有更简单的积分限制 [-1;1] 并且高斯-勒让德正交技术很容易使用。

换句话说Sf(x,y) dxdy 被简化为1111f~(e,n) |det(J)|dedn

但在 2D 情况下,我将平面任意元素更改为平面但形状良好的正方形 2×2。

3D 4 节点元素通常不是平坦的,但我想它仍然可以使用与笛卡尔坐标系相关的 2D 坐标系进行映射。我不知道如何用 {e,n} 来表达 {x,y,z} 以及在这种情况下 Jacobi 矩阵的大小是多少(它应该是正方形)。

2D 和 3D 域

1个回答

的二维流形上集成一个函数流形分析方面的书籍(如 Munkres 的易读书籍或 Lee 的流形书籍)有助于讨论定义这种积分类型的理论。R3

假设上定义的实值函数,它是您的 4 节点 3-D 元素。fM

您要计算:

MfdS.

假设映射的函数。然后φ[1,1]2M

MfdS=[1,1]2f(φ(x,y))(det(DφT(x,y)Dφ(x,y)))1/2dxdy

(我用这套笔记来刷新我的记忆。)上面,是 \varphi 的雅可比矩阵是它的转置。DφφDφT

上写出积分,你就可以使用数值方法来评估它。[1,1]2

一些评论:

  • 我很确定您的 4 节点 3-D 元素是一个流形。如果是,则函数存在(根据定义),是分段连续的(对于拓扑流形),并且是可逆的。您可以自行决定是否找到具有这些属性的函数。φ
  • 上述论证假设是一个光滑流形,这意味着存在一个连续可微在您的情况下,您描述的元素可能不是连续可微的。如果这是真的,你可能仍然可以将你的流形划分为两个光滑的流形,然后上面的论点仍然成立。同样,你必须找到满足可逆性和连续可微性的性质的Mφφ