奇异值计算mathematica_应力奇异问题的处理及子模型应用

在进行应力计算时，我们通常会在应力集中的位置加密网格来看该处的应力是否有变化，若应力变化较大，则认为计算结果没有收敛。我们需要继续加密网格再次求解，直到应力不再随网格数量的增加而变化。这样的应力结果才是mesh independent 的。

但是有的时候，在逐渐加密网格的时候，我们会发现应力会一直增加。这时候，我们需要注意应力集中的位置是否出现了应力奇异。

应力奇异是受力体由于几何关系，在求解应力函数的时候出现应力无穷大的情况。

在实际情况下，任何物体都是有一定的强度的，因此不会出现应力无穷大的现象。

因此，应力奇异是不符合实际的，我们不能够使用一个奇异的应力值去评定零件的强度。我们在计算的时候要能够判断应力集中位置是否出现了应力奇异现象，如果出现了要如何去解决。

一般在有限元计算中，应力奇异通常会出现在几何不连续的位置，如没有圆角过渡的拐角处，受力的边缘等。

下面通过一个简单的阶梯轴案例来说明应力奇异的现象与它的处理方式

1. 背景

阶梯轴两端受向下的300KN的力，轴上两端黄色线上添加远端位移约束，计算得到的阶梯轴的连接处有应力集中