结合案例图分析，我把这个项目的完整过程、参数设置、结果解读和踩坑经验整理如下。
客户做的是带孔平板在单向拉伸下的应力集中分析。他要的是孔边最大应力、应力集中系数 K_t，以及应力随距离衰减的规律，用于校核开孔板强度。
我用的软件是 ANSYS Mechanical。模型是二维平面应力板，尺寸 200×100 mm，中心圆孔直径 20 mm（孔径与板宽比 0.2），远场单向拉伸应力 σ_∞=100 MPa，材料为线弹性钢（E=210 GPa，ν=0.3）。
参数设置：单元 PLANE183（8 节点四边形），孔边网格加密到 0.2 mm，远场粗化到 5 mm；由于对称性，取 1/4 模型，对称面上施加对称边界条件。
第一张案例图是板的 Mises 应力云图。颜色从蓝到红代表应力从低到高。可以看到孔边出现明显的应力集中，最大应力区位于孔边与拉伸方向垂直的位置（即 θ=90° 附近），颜色为黄色/红色。图上标注了应力集中系数 K_t = σ_max / σ_∞ = 3.0，说明孔边最大应力是远场应力的 3 倍。
第二张案例图是沿孔边径向的 Mises 应力衰减曲线。横轴是归一化径向距离 r/a（a 为孔半径），纵轴是 Mises 应力（MPa）。在 r/a=1.0（孔边）处应力为 300 MPa；随着 r/a 增大，应力迅速下降；到 r/a≈3 时已接近远场应力 100 MPa；r/a>4 后基本恒定在 100 MPa。这个衰减规律与弹性理论中的 Kirsch 解一致。
第一，孔边网格必须极密。我最初用 1 mm 网格，K_t 只有 2.7；加密到 0.2 mm 后才得到 3.0，和理论解吻合。应力集中区的梯度极大，网格不加密会严重低估峰值。
第二，1/4 模型的对称边界不能错。我有一次把对称面设成了固定约束，相当于给板加了额外约束，导致应力分布不对称。正确做法是只约束法向位移。
第三，K_t 只和几何有关，和远场应力大小无关。这个结论对客户很重要——只要孔径/板宽比不变，无论载荷多大，K_t 都是 3.0；破坏是否发生取决于材料屈服强度和峰值应力的比较。
下次再做应力集中分析，我会固定这个流程：几何简化 → 1/4 或 1/2 对称模型 → 孔边局部加密 → 线弹性解 → 和理论 K_t 对标。如果峰值应力超过材料屈服强度，我会建议客户做弹塑性分析或大变形分析，因为线弹性结果会给出不真实的无限大应力。
计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验 需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。
