等参单元的雅可比矩阵行列式与单元面积的关系

今天给大家分享的是等参单元中的雅可比矩阵行列式与单元面积的关系。

OK，先来普及一下何为雅可比矩阵，维基 百科给出的解释是这样的：

在向量分析中，雅可比矩阵（也称作Jacobi矩阵，英语：Jacobian matrix）是函数的一阶偏导数以一定方式排列成的矩阵。—维基 百科

身为非数学专业的我们是不是感觉上述概念有些许抽象，在有限元的世界里，雅可比矩阵多出现与等参变换的过程中。

坐标系转换

大家可看一眼下方的图示，简要的表述了等参单元的自然坐标系与物理坐标系变换的关系。

 

左侧为等参单元在自然坐标系下的表示，为一个规则的正方形单元，边长为2；右侧为物理坐标系下的任意形状四边形单元。两者坐标系通过特殊的变换，即可由其中一个坐标系代表另一个坐标系。

雅可比矩阵马上登场！

其中，雅可比矩阵可表示为：

形函数对两个坐标系坐标的偏导，可由雅可比矩阵进行转换。至此，雅可比矩阵的概念已成功引出，再回到今天的主题：等参单元中的雅可比矩阵行列式与单元面积的关系

面积微元转换

既然两者坐标系可以相互转化，那坐标系下的单元是不是也可以相互转化？

对于单元的面积：自然坐标系下的等参单元的微元面积可表示为    ：

在这里补充一个数学知识：两个向量叉乘数值等于与两向量共线的平行四边形的面积。

物理坐标系下的    、    方向可表示为    、    ，两方向的向量可表示为：

单元面积进而可以表示为：

雅可比行列式就此登场！

比如，常应变单元中的单元面积可在三节点等参单元区域内进行二重积分求解：

对于等参四边形单元：

相同的道理，对于等参六面体单元：

通过上述变换，公式左侧为物理坐标系下的单元面积，1/2、4、8均为相应等参单元的面积、体积，两者通过雅可比行列式相连接。

声明：文中的公式部分参考了张雄老师的《有限元法基础》。