力学有限元分析中的网格划分策略与应力集中处理

网格决定结果的80%

在力学有限元分析中，网格划分不是前置步骤，而是决定结果质量的核心环节。本项目在多个工程项目中发现：同样的几何模型和边界条件，粗网格和细网格的应力结果可以相差300%以上——特别是在应力集中区域，差异更为悬殊。

一个实际案例：某法兰连接结构，初始用自动网格划分（平均尺寸5mm），最大应力为185 MPa。在螺栓孔周围加密到1mm后，最大应力跳到420 MPa——超过了材料屈服强度。这个差异意味着结构从”安全”变成了”危险”，完全是网格密度改变的结果。

关于有限元分析的完整方法体系，站内有系列文章。本篇聚焦网格层面的工程实践。

四面体 vs 六面体：不是随便选的

混合网格策略

实际工程中，本项目采用混合网格策略：

• 远场区域（应力梯度小）：粗四面体，尺寸10-20mm
• 过渡区域：中等四面体，尺寸5-10mm
• 应力关注区域：细六面体或四面体，尺寸1-3mm
• 应力集中点（孔、缺口、圆角）：局部加密到0.5-1mm

ANSYS中的实现方式：

1. 全局用四面体自动划分
2. 在关注面/体上施加Body Sizing
3. 在应力集中点施加Refinement
4. 在圆角处施加Face Meshing（转六面体）

应力集中区域的网格加密

孔洞周围

平板中心圆孔的应力集中系数K_t=3（理论值），有限元要准确捕捉这个系数，孔边的网格密度有明确要求：

D/d=孔径/孔边网格尺寸。本项目要求至少D/d=20，重要结构要求D/d=40。

圆角/倒角

圆角是另一个应力集中高发区。网格要求：

• 圆角沿弧长方向至少8-12个单元
• 圆角法向方向至少3-5层边界层网格
• 第一层网格厚度根据应力梯度估算

ANSYS中的操作：在圆角面上施加Edge Sizing，指定Number of Divisions=12，Behavior=Hard（强制执行）。

裂纹尖端

裂纹尖端的应力奇异性需要特殊处理：

网格无关性验证

标准流程

本项目执行网格无关性验证的标准步骤：

收敛曲线解读

应力奇异性的处理

当网格无限加密时应力持续增大（不收敛），通常意味着应力奇异性——几何上的尖角（零半径圆角）、点载荷、固定边界的角点等都会产生数学上的应力无穷大。

处理策略：

高阶单元与低阶单元的选择

线性vs二次的精度差异：以悬臂梁弯曲为例，Hex8在粗网格下会出现剪切锁定（shear locking），位移偏小50%以上。改用Hex20或增强应变公式(ANSYS中的SOLID185+keyopt(2)=1)可以消除这个问题。

本项目的默认选择：Tet10作为通用单元，Hex8/Hex20用于规则区域和高精度需求。

收敛性问题的网格因素

力学有限元分析中常见的非线性收敛问题，很多与网格有关：

实际工程案例

法兰螺栓连接结构

网格无关性验证结果：

中网格与细网格偏差4%，满足<5%的工程要求，选定中网格(40,000单元)作为最终计算网格。关于科研学术网中更多有限元分析案例，站内有专题系列。

力学有限元分析中网格划分不是一个可以自动完成的步骤——它需要分析者理解应力分布的物理图景，在关键区域投入网格资源，在非关键区域节省计算量。本项目始终坚持”先粗后细、先整体后局部、先验证后交付”的原则，确保每一份分析报告的数值结果都经得起网格无关性检验。