怎样用matlab作应力圆？

Question

如题所述

Answer 1

推荐正应力为横轴，切应力为纵轴且向下。如此在摩尔圆上旋转跟旋转微元方向一致，角度是2倍。

2维应力状态摩尔圆做法：

方法1，应力矩阵对角化求得两个主应力，两个主应力在横坐标轴上标出，作为端点做圆即为摩尔圆。

方法2，    摩尔圆圆心

摩尔圆半径

至于作图法，这里要注意，摩尔在摩尔圆中对于切应力正负的定义与我们平时用的符号规则不同。摩尔定义可以使微元顺时针转动的切应力为正。这样xy切应力和yx切应力符号就相反了，我们就可以做出唯一且确定位置的直径，有了直径摩尔圆也就有了。

虽然在做图中不考虑切应力的正负问题 也不影响摩尔圆的建立。但是当前应力状态所对应摩尔圆上的点就无法圈定了。只有确定摩尔圆上的对应点，我们才能借助摩尔圆看出想要得到想要的应力状态，我们因该如何旋转微元。

3维应力状态，目前我只知道对角化应力矩阵求出3个主应力一种方法。对于复杂应力场的对角化，matlab的可以方便实现。