可视化 | 三角形常应变单元应力后处理

有限元计算结果主要包括位移和应力两个方面。在位移方面一般无需进行整理工作。但是由于应力在单元间不连续，一般要做一定的修正。比如，三角形常应变单元也是常应力单元。

这个CST单元建立的有限元模型得到的X方向的正应力（sigma_xx）分布如图所示：

为了由计算结果得到弹性体内某一结点的接近实际的应力值，通常可采用绕结点平均法或两单元平均法。所谓绕结点平均法，就是把环绕某一结点的各单元常应力加以平均，用以表示该结点的应力。为了使由这样的平均得到的应力能够较好地表示结点处的实际应力，环绕该结点的各个单元的面积不应相差太大。单元之间的尺寸相差较多时，采用单元面积作为权系数的加权平均更为合理。绕结点平均法计算出来的结果，在内结点处较好，而在边界结点处可能很差。

以下是绕结点平均法的代码：

计算结果如下：

得到结点应力后，用MATLAB画sigma_xx云图