CAE仿真中怎么处理应力集中和应力奇异？.docxVIP

1.理解概念：应力集中vs.应力奇异

应力集中：

???定义：由于几何形状的突变（如孔、槽、台阶、尖角）或载荷的突变（如集中力、约束点），导致局部区域的应力显著高于周围区域的现象。

物理本质：是真实的物理现象，存在于实际结构中。应力峰值是有限的，但可能远高于名义应力。

原因：几何不连续处会改变力的传递路径，导致局部应力流线密集。

例子：带孔平板的孔边应力、轴肩过渡处的应力、焊缝根部应力。

应力奇异：

定义：在理论模型或理想化模型中，由于几何或边界条件的数学奇异性（如理想尖角、理想点约束/点载荷），导致有限元解在奇点处的应力值理论上趋于无穷大的现象。

物理本质：不是真实的物理现象。实际结构中不存在绝对尖角（总有微小圆角）或真正的点载荷/点约束（载荷/约束总有一定的分布面积）。

???原因：数学模型本身的不完善或过度简化。在尖角或点接触处，弹性理论的解在数学上发散。

???例子：两块板粘接的理想尖角处应力、在单个节点施加集中力或完全固支约束处的应力。

关键区别：?应力集中是真实的、有限的高应力区，需要关注并可能影响设计；应力奇异是数值计算在模型缺陷点产生的虚假无穷大应力，其值本身无物理意义，需要识别并妥善处理。

2.处理应力集中（真实物理现象）

处理目标是准确预测这个有限的峰值应力值，并据此评估结构的安全性（是否超过材料强度）。

精细化网格：

???核心方法：在应力集中区域（孔边、圆角根部、缺口等）使用足够精细的网格。

???原理：有限元法通过离散化近似连续解。应力梯度大的区域需要更小的单元来捕捉应力变化。

操作：使用网格细化技术（局部加密、映射网格、边界层网格等），确保应力集中区域有3-5层单元（尤其是二阶单元）。对比不同网格密度下的结果，进行网格收敛性研究，确保峰值应力结果趋于稳定。

???注意：过度加密奇异点附近的网格对解决应力奇异无效且浪费计算资源。

使用高阶单元：

推荐：在应力集中区域优先使用二阶单元。

???原理：二阶单元（如SOLID186,SHELL181withKEYOPT(3)=2）具有更高的形函数阶次，能更准确地模拟应力梯度变化和弯曲变形，通常比一阶单元给出更精确的应力结果（尤其在粗网格下）。

正确的几何建模：

关键：准确反映实际的几何细节，特别是过渡圆角。

操作：不要忽略制造过程中必然存在的微小圆角（即使图纸未标注），在模型中用合理的半径建模。即使是小圆角，也能显著改变应力峰值的位置和大小，使其更接近真实情况。

子模型技术：

适用场景：大型复杂模型，全局模型难以在关键细节处使用足够细的网格。

方法：

先用较粗的网格运行全局模型分析。

在感兴趣的高应力区域（应力集中区）切割出一个更小的子模型。

将全局模型在该区域切割边界上的位移解（作为边界条件）插值到子模型的边界上。

在子模型中对局部细节（如小圆角、焊缝细节）进行非常精细的网格划分并重新分析。

优点：在不过度增加整体计算成本的前提下，获得关键区域高精度的应力结果。

应用疲劳或断裂力学分析：

如果应力集中可能导致疲劳失效或脆性断裂，不能仅看峰值应力是否低于屈服强度。

???使用基于局部应力应变法的疲劳分析工具。

???对于含裂纹或类裂纹缺陷的结构，应用断裂力学（计算应力强度因子K,J积分等）来评估其安全性。

3.处理应力奇异（数值计算假象）

处理目标是识别奇异点，避免误读该点应力值，并从合理距离处提取有意义的应力用于评估。

识别奇异点：

位置：理想几何尖角（如未倒圆的90°内角）、点载荷施加点、点约束（完全固支的单个节点）处。

特征：随着网格加密，奇异点处的应力值持续增大而不收敛（发散）。这是判断奇异性的最可靠方法（网格收敛性研究）。

避免基于奇异点应力做判断：

最重要原则：绝对不要将奇异点处的计算应力值（无论看起来多高或多低）作为结构强度评估的依据！这个值是无物理意义的。

修改模型以消除奇异性：

添加微小圆角：这是处理几何尖角奇异性的最有效、最符合工程实际的方法。即使添加一个很小的工艺圆角（如R0.1mm），也能使奇异性消失，应力变为有限值（虽然可能仍然很高，但这是真实的应力集中）。然后按处理应力集中的方法分析。

分布载荷/约束：避免施加点载荷或点约束。

载荷：将集中力施加在一块小面积上（如washer），或通过耦合/约束方程分配到附近多个节点上。

约束：避免将完全固支约束施加在单个点上。应施加在具有一定面积的