应力集中是工程结构强度失效的主要诱因之一。当构件存在孔洞、缺口、截面突变等几何不连续时，应力流线在局部弯折聚集，使局部应力远高于名义应力。经典弹性理论中，无限大平板中心圆孔在单向拉伸下的应力集中系数 K_t = 3.0（Kirsch 解），是应力集中分析的基准案例。孔边应力分布、应力集中系数的准确预测，对含孔结构（法兰、机身蒙皮、压力容器开孔）的强度设计与疲劳寿命评估至关重要。
有限元分析（FEA）可数值求解含孔平板的应力场，量化应力集中系数及其网格收敛性。本项目以中心开孔矩形板受单向拉伸为对象，用有限元方法分析孔边应力分布，研究应力集中系数随网格细化的收敛行为，并与 Kirsch 理论解对比。
对无限大平板中半径为 a 的圆孔，受远场单向拉伸应力 σ₀，孔边切向应力（环向应力）为 σ_θθ(θ) = σ₀(1 + 2cos2θ)。在 θ = 90°（垂直加载方向）处达最大值 3σ₀，对应应力集中系数 K_t = 3.0；在 θ = 0° 处 σ_θθ = −σ₀（压缩）。孔边应力沿周向呈”拉-压-拉”交替分布，应力集中在孔边最剧烈，随距离迅速衰减。
应力集中系数定义为 K_t = σ_max/σ_nominal，其中 σ_max 为局部最大应力，σ_nominal 为远离孔口的名义应力。对有限宽板，K_t 高于 3.0（宽度约束效应），本案例宽径比 W/D = 5 时 K_t 约 2.85。应力集中系数是疲劳设计的核心参数——疲劳裂纹几乎总是在应力集中处萌生。
有限元解的精度依赖网格密度，尤其在应力集中区域需要细密网格。网格收敛性验证通过逐步加密网格观察关键量（K_t）的变化：随着单元数增加，K_t 逐渐收敛到理论值，相邻两次网格的 K_t 变化小于 1% 时可认为网格无关解已获得。本案例中 FEA 收敛值 2.85 与理论值 3.0 的差异主要来自有限宽度效应。
图 18 的左面板给出应力集中系数 K_t 随网格细化的收敛曲线（橙色圆点，对数横轴）。单元数从 500 增至 24000，K_t 从 2.55 单调上升并趋于平台，最终收敛至 2.85（与理论 Kirsch 解 3.0 的红线仅差约 5%）。收敛曲线的形态清晰展示了有限元方法的收敛特性：粗网格低估应力集中（应力梯度未被解析），细网格逐步逼近真实解。2.85 与 3.0 的残余差异主要源于本案例为有限宽板而非无限大平板（宽度约束使孔边应力略有重分布），属预期偏差而非计算误差。该收敛性验证是有限元分析质量保证的标准环节，证明后续基于该网格的应力结果可靠。
图 18 的右面板给出孔周环向应力的等值线云图（应力色标从深蓝 −100 MPa 到深红 285 MPa）。应力场呈典型的”上下红、左右蓝”的蝴蝶分布：在孔边水平位置（θ = 90°，垂直加载方向）出现最大拉应力约 285 MPa（红色区域），在孔边竖直位置（θ = 0°）出现约 −100 MPa 的压应力（蓝色区域）。应力自孔边向外快速衰减，远离孔口区域恢复至均匀的 100 MPa（绿色）。云图直观复现了 Kirsch 解的”3σ₀ 峰值 + 拉压交替”特征，验证了有限元模型的正确性。285 MPa 的局部应力远高于名义应力，说明若按名义应力设计而忽略应力集中，将严重低估局部应力水平。
两面板结合，左面板验证了网格收敛性（数值可靠性），右面板给出了应力分布的全貌（物理正确性），共同完成含孔平板静力学分析，为含孔结构的强度与疲劳设计提供了关键应力集中数据。
本项目对中心开孔平板进行了有限元静力学分析，应力集中系数收敛于 2.85（理论 3.0，有限宽度修正），孔边最大应力 285 MPa（名义 100 MPa）。案例图以 K_t 收敛曲线与孔周应力云图双面板呈现。
后续可拓展：(1) 对不同孔径比 W/D 系列计算 K_t，建立设计曲线；(2) 结合 S-N 曲线开展含孔构件疲劳寿命预测；(3) 研究孔边补强（加厚环、加强筋）对应力集中的缓解效果。
