应力集中与应力奇异：别再傻傻分不清
CAE摸鱼指南 · 第07期
做仿真的朋友，应该都遇到过这种情况：算完一看结果，某个地方应力特别高，红得吓人。这时候问题来了——这到底是真实的应力集中，还是数值上的应力奇异？
很多刚入行的工程师，经常把这俩混为一谈。结果要么是把奇异当集中，白担心一场；要么是把集中当奇异，漏掉了真正的风险点。
今天咱们就把这两个概念掰扯清楚，从理论到实践，争取看完就能用。
一、应力集中：真实存在的物理现象
先说应力集中。这个是真实存在的，物理上说得通。
简单讲，当构件的几何形状发生突变时，比如有孔、有槽、有台阶，应力流线会在这些地方发生弯曲、密集，导致局部应力显著高于周围区域。这就像河流遇到桥墩，水流会在桥墩附近变急一样。
图1 带圆孔平板拉伸应力云图：孔边应力明显集中
1.1 应力集中系数
衡量应力集中程度的指标，叫理论应力集中系数（Theoretical Stress Concentration Factor），用 Kt 表示。定义很简单，就是最大应力和名义应力的比值：
Kt = σmax / σnom
式中：σmax — 局部最大应力；σnom — 截面名义应力
这里的名义应力，一般是按净截面面积算的平均应力。Kt 永远大于1，越大说明集中越严重。
1.2 经典案例：无限大平板圆孔
最经典的例子，就是无限大平板中间开个圆孔，单向受拉。这个是有解析解的，由基尔施（Kirsch）在1898年推导出来。
结果是：孔边垂直于拉伸方向的环向应力最大，等于远场应力的3倍。也就是说：
σθ,max = 3σ∞
即：Kt = 3
这个3倍，是工程上非常经典的一个数值。很多时候估算孔边应力，直接拿名义应力乘3，差不了太多。
当然，如果不是圆孔而是椭圆孔，情况又不一样。椭圆孔的应力集中系数和长短轴比有关：
Kt = 1 + 2a/b
式中：a — 椭圆长半轴（垂直于载荷方向）；b — 椭圆短半轴
可以看到，孔越扁（a/b越大），应力集中越严重。这也解释了为什么裂纹尖端应力会那么大——因为裂纹可以看成是极扁的椭圆。
1.3 应力集中的特点
总结一下，应力集中有这么几个特点：
物理真实
：是实际存在的现象，实验可以测量到 峰值有限
：最大应力是有限值，有明确的上限 局部性
：只发生在几何突变附近，离远了就恢复均匀了 可收敛
：有限元计算中，随着网格加密，应力会收敛到一个稳定值
这最后一点很重要，是区分应力集中和应力奇异的关键。后面会详细说。
二、应力奇异：数学上的假象
说完应力集中，再说说应力奇异。这个就不一样了，它不是真实的物理现象，而是数学模型理想化的副产品。
什么意思呢？就是说，我们在建立有限元模型的时候，做了一些理想化的假设，这些假设在数学上会导致某些点的应力趋于无穷大。但在真实世界里，这种情况是不存在的。
图2 L形截面内尖角应力奇异：尖角处应力极高
2.1 应力奇异的来源
常见的应力奇异来源，主要有这么几类：
第一类：几何尖角
这是最常见的。比如两个面垂直相交，形成一个完美的90度内角。在弹性力学里，这种尖角的应力理论上是无穷大的。因为曲率半径为零，应力集中系数就变成了无穷大。
但实际加工中，不可能有真正的尖角，多少都会有一点圆角。只是我们在建模的时候，为了省事，直接画成了尖角。
第二类：点约束/刚性约束
在一个节点上施加固定约束，或者在一个点上施加集中力，也会产生应力奇异。因为载荷被施加在一个无穷小的面积上，应力自然就是无穷大。
同样，真实情况中，约束和载荷都是作用在一定面积上的，不会是真正的"点"。
第三类：材料界面尖角
两种不同材料的界面，如果在边界处形成尖角，也会产生奇异。这个相对少见一些。
第四类：裂纹尖端
这个比较特殊。裂纹尖端确实是应力奇异的，但这个奇异是有物理意义的，是断裂力学研究的对象。我们可以用应力强度因子K来描述，而不是直接看应力值。
2.2 应力奇异的特点
应力奇异的特点，刚好和应力集中对着来：
非物理
：是数学模型的产物，真实世界不存在 峰值无限
：理论上应力趋于无穷大 不收敛
：网格越细，应力越大，没有上限 高度局部化
：只发生在个别节点上，范围极小
⚠️ 划重点：应力奇异点的应力值本身没有物理意义，不能用来评估结构强度。拿奇异点的应力去和屈服强度比，那是完全错误的。
三、核心区别：一张表讲清楚
说了这么多，可能有人还是有点晕。没关系，咱们列个表对比一下，一目了然：
四、怎么判断？黄金标准是网格加密
道理都懂，但实际拿到一个计算结果，怎么判断是集中还是奇异呢？
最可靠的方法，也是行业公认的黄金标准，就是网格加密测试。
做法很简单：把你怀疑的区域的网格加密一倍，再算一遍，然后看应力怎么变。
图3 网格收敛对比：左图应力集中收敛，右图应力奇异发散
如果是应力集中：网格加密后，应力会先升一升，然后逐渐趋于稳定，最后收敛到某个值。这个值就是真实的应力集中峰值。
如果是应力奇异：网格越细，应力越大，而且看不到收敛的趋势。理论上，只要网格足够细，应力可以无限大。
💡 实用建议：
至少做两级网格加密，也就是三个不同密度的网格，看应力的变化趋势。只看两级有时候不够准确。
另外，还要注意看应力的分布范围。如果高应力只集中在一两个节点上，周围很快就降下来了，那大概率是奇异；如果高应力有一定的区域范围，那更可能是集中。
五、怎么处理？对症下药
判断清楚了，接下来就是怎么处理。应力集中和应力奇异的处理思路完全不一样，千万别搞混了。
5.1 应力集中的处理
应力集中是真实存在的，所以不能回避，得正面应对。常用的方法有：
1. 几何优化
这是最根本的办法。比如把尖角改成圆角，把直角台阶改成斜坡，把小孔改成大孔。总的原则就是让几何变化尽量平缓，不要太突兀。
圆角半径越大，应力集中系数越小。一般来说，r/d（圆角半径/直径）大于0.1以后，改善效果就比较明显了。
2. 局部加强
如果几何改不了，可以在应力集中的地方局部加厚、加筋、或者用强度更高的材料。
3. 合理评估
不是所有应力集中都要处理。要看具体工况：如果是静载，材料塑性好，应力集中的影响就小一些，因为塑性变形会让应力重新分布；如果是疲劳载荷，那应力集中的影响就很大，必须认真对待。
5.2 应力奇异的处理
应力奇异是假的，所以处理思路就是"去假存真"——把导致奇异的理想化假设去掉，让模型更接近真实情况。
图4 应力奇异处理：尖角变圆角，应力分布更合理
1. 几何修正
最常用的就是给尖角加圆角。只要加一个很小的圆角，奇异就消失了。至于圆角半径取多大，要看实际加工情况，不能瞎加。如果实际就是尖角，那你可以加一个工艺圆角，比如R0.5、R1，然后说明这是建模处理。
2. 约束等效
把点约束改成面约束，或者用耦合、约束方程之类的方式，把约束分散到一定区域。不要在单个节点上钉死。
3. 载荷分散
把集中力改成分布力，作用在一定的面积上。真实的载荷本来就不是作用在一个点上的。
4. 直接忽略
如果奇异点不在你关注的区域，而且离关键位置有一定距离，那可以直接忽略它。根据圣维南原理，局部的扰动不会影响远处的应力分布。
这时候你可以看奇异点旁边一两个单元的应力，那个更接近真实值。
📌 圣维南原理：
简单说就是：如果把物体一小部分边界上的力，换成静力等效的另一组力，那么只会在这一小部分附近产生明显的应力变化，对远处的应力影响很小。
这个原理是我们忽略局部奇异的理论依据。
六、总结与建议
🎯 核心要点总结
本质不同
：应力集中是物理真实，应力奇异是数学假象 判断方法
：网格加密测试，看应力是否收敛（黄金标准） 处理思路
：应力集中要正面应对，应力奇异要修正模型或忽略 关键公式
：Kt = σmax / σnom；圆孔Kt=3；椭圆孔Kt=1+2a/b
最后给大家几点工程上的建议：
第一，拿到结果先别急着下结论。看到高应力，先想想这地方是不是有尖角、点约束之类的，会不会是奇异。
第二，养成做网格收敛性验证的习惯。不是说每个项目都要做，但至少对关键位置、关键结果，要心里有数。
第三，建模的时候尽量贴近实际。不要图省事到处画尖角、加点约束，不然后面处理起来更麻烦。
第四，区分清楚"应力大"和"应力奇异"。应力大不一定是奇异，奇异一定是应力大，但反过来不成立。
应力集中和应力奇异，是有限元分析里最基础也最容易混淆的概念。搞清楚这俩，才算真正入了仿真的门。
希望今天这篇文章，能帮大家把这个知识点夯实。以后再遇到高应力，心里就有底了。
🐟 摸鱼工程师说两句
其实做仿真这行，最怕的不是算错，而是不知道自己算错了。应力集中和应力奇异，就是最容易踩坑的地方之一。
刚入行的时候，我也经常对着一个红彤彤的应力云图发呆，不知道这到底是真的有问题，还是模型本身的毛病。后来踩的坑多了，慢慢就有感觉了——看到尖角先打个问号，看到点约束先留个心眼，拿不准就加密网格试试。
仿真这东西，说到底就是"模型近似+数值近似"，没有绝对的对，只有相对的合理。搞清楚哪些是真实的物理现象，哪些是数学上的假象，才能对结果心里有数。
好了，今天就聊到这儿。下期咱们接着摸鱼~
🐟 摸鱼工程师的职场感悟
做仿真这行，最怕的就是"知其然不知其所以然"。很多人算完了，结果一看应力大，就慌了，要么直接报给领导说不行，要么就瞎调参数。
其实很多时候，应力大不一定是结构真的有问题，可能只是你模型建得有问题。尖角、点约束、集中力，这些随手一画的东西，最后都可能变成应力奇异，给你虚假的警报。
所以说，前处理是根，后处理是果。根不正，果肯定歪。多花点时间把模型建扎实了，比后面瞎调半天强得多。
毕竟，摸鱼也要摸得有技术含量，对吧？
— — — END — — —
