ANSYS强度仿真：SA-516-70圆柱壳接管开孔与内压应力分析

项目过程

本项目为某压力容器设计院评估圆柱壳接管区域的强度裕量。容器主体为圆柱壳，内径D=1000 mm、壁厚t=12 mm，材料SA-516-70（压力容器用碳钢板），屈服强度S_y=260 MPa。壳体上开有Ø170 mm径向接管，设计内压最高2.6 MPa（测试压力），工作压力更低。客户需要确认开孔处是否满足强度与疲劳要求，为开孔补强设计提供依据——此前一台类似容器在开孔接管焊缝处出现泄漏，怀疑应力集中未被充分评估。团队采用ANSYS建立轴对称整体模型+局部三维实体子模型，计算环向应力、纵向应力及von Mises等效应力。

材料设为线弹性，E=200 GPa、ν=0.3（高温下可降，本项目按常温考核）。内压施加在壳体内表面及接管内表面，边缘按对称边界处理。整体模型用壳单元快速得到薄膜应力，再用实体单元子模型（submodeling）在接管与壳体相贯线附近加密网格（局部尺寸0.5 mm），提取局部峰值应力。为捕捉几何不连续应力集中，接管根部需足够细的网格与恰当的单元类型（二阶单元）。求解后输出：壳体中面环向/纵向薄膜应力随内压变化曲线，以及接管开孔区域von Mises应力云图（含峰值位置与数值），用于按JB 4732或ASME VIII-2进行应力分类（一次薄膜、一次弯曲、二次应力）与校核。团队还对比了”无补强”与”加补强圈”两种方案下的峰值应力。

案例分析

案例图分左右两幅。左图是壳体薄膜应力随内压变化：环向应力σ_θ（红）斜率最高，2.6 MPa时约104 MPa；纵向应力σ_z（蓝）约52 MPa（约为环向一半）；von Mises等效应力（绿虚）介于两者之间约90 MPa。右图是接管开孔区域von Mises应力云图，壳体薄膜约67 MPa，接管开孔边缘（相贯线）峰值应力达160 MPa。

左图符合薄壁圆筒理论：σ_θ=pr/t、σ_z=pr/2t。代入p=2.6 MPa、r=500 mm、t=12 mm，σ_θ=2.6×10⁶×0.5/0.012≈108 MPa（图中104 MPa含封头边缘效应修正），σ_z≈52 MPa，von Mises≈90 MPa，均远低于屈服260 MPa。右图160 MPa峰值出现在接管与壳体相贯线——典型的几何不连续应力集中，集中系数约160/67≈2.4。按应力分类，该峰值属一次局部薄膜+弯曲应力，许用值可取1.5S_y=390 MPa（ASME一次应力限制），160 MPa远未超限，但疲劳评定需按交变应力幅另行考核。进一步量化：峰值应力/屈服≈160/260≈0.62，强度裕量充足；但开孔削弱了壳体，若不开补强圈，应力集中系数会更高（可达3.5以上）。客户据此在接管处加补强圈（按等面积法），将相贯线峰值压到约130 MPa，进一步提高安全系数，并改善了焊缝处的高应力状态，直接回应了此前泄漏隐患。

总结教训

第一，团队最初用全壳单元建模，接管根部无法捕捉峰值应力（只得到67 MPa薄膜值）；改用实体子模型后才得到160 MPa峰值，说明薄壳理论不适用局部不连续区，子模型是标准做法。第二，未考虑焊接残余应力时峰值被低估20%以上，后续疲劳评定需叠加焊接拉应力（可达0.3S_y），否则疲劳寿命预测过于乐观。第三，SA-516-70在低温下韧性下降，除屈服校核外还需冲击韧性（CVN）验证，尤其北方或低温介质工况。第四，网格在相贯线加密不足会让峰值偏低30%，团队用0.5 mm局部尺寸并做网格无关性验证（加密至0.3 mm变化<3%）。第五，接管与壳体厚度比、坡口形式影响应力流线，详细设计需建真实坡口几何而非理想化直角，否则峰值位置与数值偏差。

项目经验

1. 薄壁圆筒环向应力≈pr/t：2.6 MPa时σ_θ≈108 MPa，与图中104 MPa吻合，是壳体强度速算基础，纵向应力约为其半。
2. 接管开孔必须用实体子模型：壳单元严重低估局部峰值（67 vs 160 MPa），应力集中不可忽略，子模型法是规范路径。
3. 峰值/屈服≈0.62仍有裕量：但疲劳需按交变幅单独评定，开孔区建议加补强圈降应力，回应焊缝泄漏隐患。
4. 网格无关性必做：相贯线局部尺寸需加密至峰值稳定（变化<3%），否则强度结论不可信，易漏高风险点。
5. 焊接残余应力须计入疲劳：静强度够不代表疲劳够，焊缝处叠加残余拉应力后交变幅显著升高，需单独评。

更多ANSYS压力容器与强度仿真经验，可参考科研学术网ANSYS仿真栏目。