CAE仿真中怎么处理应力集中和应力奇异？
基于《何为应力集中、何为应力奇异——材料力学中的关键现象》帖子的反馈需求，下面讲述下在CAE仿真中怎么处理应力集中和应力奇异？
1. 理解概念：应力集中 vs. 应力奇异
应力集中：
定义： 由于几何形状的突变（如孔、槽、台阶、尖角）或载荷的突变（如集中力、约束点），导致局部区域的应力显著高于周围区域的现象。
物理本质： 是真实的物理现象，存在于实际结构中。应力峰值是有限的，但可能远高于名义应力。
原因： 几何不连续处会改变力的传递路径，导致局部应力流线密集。
例子： 带孔平板的孔边应力、轴肩过渡处的应力、焊缝根部应力。
应力奇异：
定义： 在理论模型或理想化模型中，由于几何或边界条件的数学奇异性（如理想尖角、理想点约束/点载荷），导致有限元解在奇点处的应力值理论上趋于无穷大的现象。
物理本质： 不是真实的物理现象。实际结构中不存在绝对尖角（总有微小圆角）或真正的点载荷/点约束（载荷/约束总有一定的分布面积）。
原因： 数学模型本身的不完善或过度简化。在尖角或点接触处，弹性理论的解在数学上发散。
例子： 两块板粘接的理想尖角处应力、在单个节点施加集中力或完全固支约束处的应力。
关键区别： 应力集中是真实的、有限的高应力区，需要关注并可能影响设计；应力奇异是数值计算在模型缺陷点产生的虚假无穷大应力，其值本身无物理意义，需要识别并妥善处理。
2. 处理应力集中（真实物理现象）
处理目标是准确预测这个有限的峰值应力值，并据此评估结构的安全性（是否超过材料强度）。
精细化网格：
核心方法： 在应力集中区域（孔边、圆角根部、缺口等）使用足够精细的网格。
原理： 有限元法通过离散化近似连续解。应力梯度大的区域需要更小的单元来捕捉应力变化。
操作： 使用网格细化技术（局部加密、映射网格、边界层网格等），确保应力集中区域有3-5层单元（尤其是二阶单元）。对比不同网格密度下的结果，进行网格收敛性研究，确保峰值应力结果趋于稳定。
注意： 过度加密奇异点附近的网格对解决应力奇异无效且浪费计算资源。
使用高阶单元：
推荐： 在应力集中区域优先使用二阶单元。
原理： 二阶单元（如SOLID186, SHELL181 with KEYOPT(3)=2）具有更高的形函数阶次，能更准确地模拟应力梯度变化和弯曲变形，通常比一阶单元给出更精确的应力结果（尤其在粗网格下）。
正确的几何建模：
关键： 准确反映实际的几何细节，特别是过渡圆角。
操作： 不要忽略制造过程中必然存在的微小圆角（即使图纸未标注），在模型中用合理的半径建模。即使是小圆角，也能显著改变应力峰值的位置和大小，使其更接近真实情况。
子模型技术：
适用场景： 大型复杂模型，全局模型难以在关键细节处使用足够细的网格。
方法：
先用较粗的网格运行全局模型分析。
在感兴趣的高应力区域（应力集中区）切割出一个更小的子模型。
将全局模型在该区域切割边界上的位移解（作为边界条件）插值到子模型的边界上。
在子模型中对局部细节（如小圆角、焊缝细节）进行非常精细的网格划分并重新分析。
优点： 在不过度增加整体计算成本的前提下，获得关键区域高精度的应力结果。
应用疲劳或断裂力学分析：
如果应力集中可能导致疲劳失效或脆性断裂，不能仅看峰值应力是否低于屈服强度。
使用基于局部应力应变法的疲劳分析工具。
对于含裂纹或类裂纹缺陷的结构，应用断裂力学（计算应力强度因子K, J积分等）来评估其安全性。
3. 处理应力奇异（数值计算假象）
处理目标是识别奇异点，避免误读该点应力值，并从合理距离处提取有意义的应力用于评估。
识别奇异点：
位置： 理想几何尖角（如未倒圆的90°内角）、点载荷施加点、点约束（完全固支的单个节点）处。
特征： 随着网格加密，奇异点处的应力值持续增大而不收敛（发散）。这是判断奇异性的最可靠方法（网格收敛性研究）。
避免基于奇异点应力做判断：
最重要原则： 绝对不要将奇异点处的计算应力值（无论看起来多高或多低）作为结构强度评估的依据！这个值是无物理意义的。
修改模型以消除奇异性：
添加微小圆角： 这是处理几何尖角奇异性的最有效、最符合工程实际的方法。即使添加一个很小的工艺圆角（如R0.1mm），也能使奇异性消失，应力变为有限值（虽然可能仍然很高，但这是真实的应力集中）。然后按处理应力集中的方法分析。
分布载荷/约束： 避免施加点载荷或点约束。
载荷： 将集中力施加在一块小面积上（如washer），或通过耦合/约束方程分配到附近多个节点上。
约束： 避免将完全固支约束施加在单个点上。应施加在具有一定面积的区域（如面、边），或使用弹性支撑模拟实际连接刚度。如果必须在点约束，应远离高应力关注区域。
提取合理距离处的应力：
当无法完全消除奇异性时（有时是必要的简化），或者关心奇异点附近但不是奇异点本身的应力状态：
方法： 在奇异点附近（离开奇点一定距离）提取应力值。这个距离通常需要大于1-2个单元尺寸。
依据： 圣维南原理。该原理指出，载荷或约束的具体分布方式只影响其作用点附近局部区域的应力分布，在稍远的地方，应力分布主要由合力和合力矩决定。因此，离开奇异点一定距离后，应力解是可靠且有物理意义的。
操作： 绘制从奇异点沿某条路径的应力分布图，观察应力在离开奇点后如何趋于稳定。选择应力值趋于稳定的位置作为评估点。或者，提取该区域一定范围内的平均应力（如节点平均、面平均）。
关注应力趋势而非绝对值：
在奇异点附近，虽然绝对值不可信，但应力分布的趋势（如高应力区域的范围、应力梯度）仍然包含有用信息。
4.总结与关键原则
区分是关键： 首先判断你面对的是真实的应力集中还是数值的应力奇异。网格收敛性研究是判别的金标准。
处理应力集中（真实）： 精细化网格（收敛性验证）、高阶单元、准确几何（圆角！）、子模型技术、结合疲劳/断裂分析。
处理应力奇异（假象）：
识别它： 位置（尖角、点载/点约束）、行为（网格加密不收敛）。
忽略它： 绝对不用奇异点应力值做强度判断！
消除它（首选）： 添加微小圆角、使用分布载荷/约束。
规避它（次选）： 提取合理距离外的应力（利用圣维南原理）、关注应力趋势。
工程判断： CAE仿真是工具，结果解读离不开工程师对物理现象的理解和工程经验。理解模型简化的假设和局限性至关重要。
通过遵循这些原则和方法，可以有效地处理CAE仿真中的应力集中和应力奇异问题，获得更可靠、更有工程指导意义的分析结果。
