ANSYS建模：法兰接头应力分析中几何简化、接触定义与网格收敛的工程决策

一个压力容器法兰接头的应力分析项目，CAD模型从供应商处直接导入ANSYS Workbench——螺纹、倒角、测压孔全部保留，自由度超1500万。第一次求解跑了两天后，von Mises峰值应力出现在螺纹根部的奇异点上——但物理上法兰根部环向应力才应该是限制载荷。这个模型需要彻底的”减法”。

几何简化的判断逻辑遵循尺寸比数准则：特征尺寸小于整体模型最小方向尺寸的1/20时，它对全局应力分布的影响可忽略。按这个标准，M20螺栓螺纹牙根半径约0.3 mm——相对于法兰厚度40 mm远小于阈值，删除螺纹后将螺纹段替换为光滑圆柱，法兰根部应力计算偏差不超过2%。倒角半径1 mm的外缘圆角同样删除。但垫片槽底部的R1.5圆角必须保留——该处是垫片压缩过程中应力梯度最大的区域，删除后局部峰值应力会被人为抹平约15-20%。

两个特征构成了例外：法兰根部的过渡圆角（R5，特征尺寸占法兰厚度的12.5%）和螺栓孔的配合面。前者的应力集中系数在1.5至2.5之间，是结构完整性分析的核心目标，删除等于废掉了整个分析。后者的孔壁接触压力分布直接决定螺栓的预紧力传递路径——如果把配合面上的螺纹孔简化成通孔，螺栓预紧力在法兰之间的分配比例会偏离实际20%以上。

几何简化完成后，自由度从1500万降到了约320万——缩减了近80%，后续网格加密的算力余量完全释放了出来。

接触定义是装配体建模中应力偏差量级最大的来源。ANSYS Mechanical默认对导入装配体的全部接触面施加Bonded约束，但法兰端面和垫片之间绝非”粘接”——螺栓预紧力下垫片发生弹塑性压缩，法兰面之间允许微量分离和切向滑移。正确的选择是Frictional接触（μ=0.2，对钢-石墨垫片界面），法向行为设为”允许分离”（Augmented Lagrange formulation）。该项目将Bonded改为Frictional后，法兰根部的环向应力从88 MPa降至62 MPa——偏差约30%，远超几何简化引入的2-5%误差。

网格收敛性验证按”逐级加密、结果追踪”的标准流程执行。全局网格尺寸依次设为8 mm、4 mm、2 mm和1 mm，观察法兰根部von Mises应力和环向应力。从4 mm到2 mm，峰值应力变化约8%；从2 mm到1 mm，变化收窄至3.2%——进入网格无关解的收敛区间（变化率小于5%）。对于法兰根部应力过渡区，额外施加局部细化控制，确保过渡圆角弧线上至少6个二次单元覆盖。

单元类型的选择在当前算力条件下已不是需要折中的决策。十节点四面体（SOLID187 Tet10）在每棱边中点插入额外节点，位移场插值从线性提升到二次，应力精度提升一个量级，而自由度仅增加约50-80%。对比四节点四面体（SOLID185 Tet4），Tet4对弯曲主导应力计算偏硬15-25%——在法兰根部弯曲应力占主导的位置，这个偏差直接传递到疲劳安全系数判定中。该项目的默认选择在所有载荷步中统一采用Tet10。

最终模型求解后在法兰根部给出的环向应力峰值与ASME VIII-2规范的标准解析解偏差在4%以内。项目收尾时，这套”几何先减后加——网格逐级收敛——接触严格定义”的三步模型准备流程被写入了内部的结构分析建模规范。