二维有限元分析

我做过一个含圆孔无限大板的平面应力分析项目，客户想验证有限元结果和经典弹性理论中的应力集中系数是否一致。我这个项目虽然模型简单，但包含了平面应力、网格收敛和应力集中三个核心知识点。

左侧图是含孔板在单向拉伸 σ₀ 作用下的 von Mises 应力云图。板左右两端受均匀拉应力，中间有一个圆孔。颜色显示孔边应力明显高于远场，孔边左右两侧出现红色高应力区，应力集中系数 K_t 标注为约 3.0。这和弹性理论中无限大板圆孔应力集中系数 K_t=3 完全吻合，说明有限元建模和网格足够精细。

右侧图是沿板中心线的应力分布 σ/σ₀ 随位置 x 变化。远离孔的位置（x<2 和 x>6）应力比为 1，等于远场应力；在孔边（x≈3 和 x≈5）应力比突升到 3.0，然后快速衰减回 1。图中标注 K_t = σ_max/σ₀ = 3.0，这就是应力集中系数。对于椭圆孔，K_t 会更大，长轴垂直于载荷方向时 K_t = 1 + 2a/b，a、b 分别是半长轴和半短轴。

建模要点：用 ANSYS Mechanical 建 2D 平面应力模型，材料弹性模量 200 GPa，泊松比 0.3；一侧固定 x 方向位移，另一侧施加均匀拉应力。关键是在孔边加密网格，因为应力梯度大。我用的单元是 PLANE182 或 PLANE183，孔边单元尺寸控制在孔径的 1/20 以内。

踩坑经验：第一个是边界条件。如果左右两侧不加约束只加拉力，刚体位移会导致求解奇异，必须至少约束一个节点的 x 和 y 位移。第二个是平面应力假设：板厚方向应力为零，适合薄板；如果是厚板，应该用平面应变。第三个是网格收敛：孔边粗网格会低估应力集中，我用 3 种网格密度做了收敛测试，K_t 从 2.6 收敛到 3.0。

工程经验：应力集中是结构疲劳和脆性断裂的根源。设计时如果不能避免开孔，就要通过加倒角、补强圈或选择更合理的开孔方向来降低 K_t。这个案例也适合教学生：有限元不是黑箱，经典理论解是验证模型正确性的金标准。

更多 ANSYS 平面应力、平面应变和应力集中分析的实战案例，我汇总在 [二维有限元分析](https://www.keyanxueshu.com/category/fea/) 栏目里。[二维有限元分析](https://www.keyanxueshu.com/)

含孔板分析还有一个工程延伸：开孔补强。当板上必须开孔时，可以通过加补强圈、加厚孔边局部壁厚或采用椭圆孔长轴平行于载荷方向来降低 K_t。右侧应力集中曲线显示，孔边高应力区非常局限，在 x≈2-3 倍孔径以外就衰减回远场应力。这意味着补强只需要覆盖孔边 2-3 倍孔径范围即可，不需要整体加厚。有限元网格收敛上，我要求孔边至少 3-5 层单元、单元长宽比接近 1，K_t 才能收敛到 3.0。粗网格下 K_t 可能只有 2.2-2.5，会误导设计。把这些细节写进报告，客户自己也能判断仿真结果是否可信。