高斯坐标自然值计算_有限元简单科普之——高斯积分点的应用

上一篇我们已经介绍了，基于最小势能原理的单元平衡方程是有限元计算的根本。其形式就是单元刚度矩阵(element stiffness matrix)乘以位移向量等于右边的荷载向量(right hand side vector)：

其中，

为单元刚度矩阵。

为右侧的单元荷载向量。

这里的

是体积的微分，在二维问题里，写开来其实就是

（所以

当然是

啦）。 进一步放到母单元中可以写成

，那么

，右边的荷载向量也可以用Jacobian矩阵做mapping，然后积分。以上内容在上一篇已经讲得十分完备了，这里复读了一遍。

那么问题来了，这个刚度矩阵和荷载向量的积分要怎么算？譬如荷载的积分，我们知道，一个四边形单元（quadrilateral element）的每个点在两个方向上都有应力分布，而且考虑到实际情况，这个应力往往是一条曲线！譬如我们对着母单元这样横着切一刀，我作为好事者想看看这个截面上的竖直方向的正应力如何？

有同学可能要问为什么我要横着切一刀？为什么不是斜地切一刀……因为斜着切一刀的图像其实就是