FEA 中的一个基本问题-凹角应力奇异的本质

什么是凹角？

让我们用一个例子来说明这一点。这里我们有四个梁单元以一定角度排列，并标记了单元和节点。由于单元 3 沿 x 轴，其刚度矩阵涉及节点 3 和 4 处形状函数相对于 x (dN/dx) 的导数。然而，在节点 3 处，形状函数相对于 x 的导数为 0，因为元素 2 沿 y 轴。这导致刚度矩阵仅包含 0。

由于单元 2 沿 y 轴，其刚度矩阵涉及节点 2 和 3 处的形状函数相对于 y (dN/dy) 的导数。然而，在节点 3 处，形状函数相对于 y 的导数为 0，因为元素 3 沿 x 轴。这也会导致刚度矩阵仅包含 0。无论使用哪个元素来制定刚度矩阵，节点 3 处的刚度矩阵也将全部为 0，从而导致应力奇点。 

对于凹角该怎么办？

“在远离载荷应用区域的物体上的点处产生的应力和应变将与具有相同静态等效结果并施加到同一区域内的物体的任何施加载荷产生的应力和应变相同”。

如果横截面允许，对凹角处的截面力进行手动计算通常也会给出更准确的结果，特别是在运行线性静态分析时。

如果凹角位于所关注的高应力区域，另一种方法是添加小半径，从而将应力奇点转变为应力集中。这将使网格正确收敛。此外，在几乎所有实际应用中，由于制造限制，实际上并不存在完美尖锐的凹角，因此添加较小的半径是有效的。

在此示例中，横截面直径为 10 毫米、长度为 10 毫米的轴在 X 方向上加载 500 N 的力。该分析是在线性静态假设下进行的。

在下面的轮廓中，轴和较厚的部件（直径为 30 毫米）之间没有半径。就像本博客前面的张力示例一样，拐角处的应力值会发散，因为没有半径，因此存在凹角。

上表将 FEM 的应力值与距凹角 1mm 处手工计算的应力进行了比较。正如您所看到的，在距拐角 1mm 处，FEM 应力与手工计算的应力相当吻合，而 FEM 应力在拐角处发散。手动计算应力的公式为

作为直接从 FEM 探测应力的替代方法，您可以拉动截面力，这将为您提供更现实的答案，该答案将随着网格变得更精细而收敛。

上表将根据截面力计算出的应力值（使用与之前相同的公式）与手动计算的凹角处和距凹角 1mm 处的应力进行了比较。与上表中的有限元力相比，两个位置的截面力均趋向于手动计算值。然而，有一个非常重要的警告 - 尽管凹角处的应力在现实生活中不会是无限的，但这些角仍然会是应力集中，因此应力会升高。该升高的应力相当于手动计算的应力乘以应力集中系数 (Kt)。如果存在已知或假设的半径，则可以使用Peterson 应力集中因子等资源来量化应力集中因子并与截面力结合使用来查找真实应力，这在下面的示例中完成。

在下面的轮廓中，同一轴添加了 0.5 mm 的半径。

上图显示，当添加小半径时，应力值会收敛。在这种情况下，由于半径引起的应力集中，这些应力将与截面力不匹配；然而，这种应力集中可以使用Peterson 的应力集中因子进行量化。下表中，之前拐角处的手计算应力乘以 Kt = 2.39，这是使用Peterson计算得出的。