UG有限元结构强度仿真：45钢带圆角L形拐臂应力集中与收敛性复现

本项目针对一款 45 钢制造的 L 形拐臂（典型机械传动结构件）开展带圆角的静强度有限元分析，目标是在不同圆角半径 R 下准确预测峰值 von Mises 应力、安全系数以及网格收敛性，为设计优化提供依据。零件材料为 45 钢（优质碳素结构钢，调质态），弹性模量 E = 206 GPa、屈服强度 σ_y = 355 MPa（案例图（a）标注），几何特征为 L 形拐臂，拐角处设圆角 R，待优化变量。载荷与约束：拐臂一端承受拉力 F = 12 kN（沿长边方向），另一端固定约束。本项目采用 Siemens NX（UG）的”Advanced Simulation”模块（基于 NX Nastran 求解器）建模，使用 3D 实体单元（线性四面体 C3D10 与二次四面体 C3D10H 混合），主要网格由二次四面体主导，整体单元数约 3.8 万；圆角区域网格加密至 0.5 mm 以捕获应力梯度。扫描圆角 R = 1、2、5、10、15、20 mm 共 6 个工况，重点复现案例图（b）中 R = 15 mm 处的 σ_vm,max = 296 MPa、K_t = 1.51、n = 1.20 关键数据；案例图（c）展示网格收敛性研究，覆盖单元尺寸 h = 4.0、2.0、1.0、0.5 mm。

案例分析

案例图（a）展示带 R = 15 mm 圆角的 L 形拐臂 von Mises 应力云图：峰值应力 σ_vm,max = 296 MPa 出现在圆角根部（图中明确标注”σ_vm,max = 296 MPa at fillet R = 15 mm”），位置在 L 拐角内侧 x = 15 mm、y ≈ 0 处（图中圆点标注），与应力集中理论吻合。远离圆角的应力下降迅速，距离圆角 30 mm 以上区域应力低于 80 MPa。本案例应力集中系数 K_t = σ_vm,max / σ_nominal = 1.51（图中明确标注），相比无圆角直角拐臂（K_t ≈ 2.5）下降约 40%，验证圆角对应力集中的显著缓解作用。案例图（b）给出 σ_vm,max 与安全系数 n 随圆角半径 R 的关系曲线：R 由 1 mm 增至 20 mm，σ_vm,max 由约 720 MPa 单调下降至 280 MPa，安全系数 n（= σ_y/σ_vm,max）由约 0.49 上升到 1.27；设计点 R = 15 mm 处 σ_vm,max = 296 MPa、n = 1.20（图（b）右下方框明确标注”n = 1, 296 MPa”）。σ_y = 355 MPa 参考线在图中横线标识。曲线呈现典型幂律下降规律：σ_vm,max ∝ R^(-0.2)，即圆角每翻倍，应力下降约 13%。案例图（c）展示网格收敛性研究：在单元尺寸 h = 4.0、2.0、1.0、0.5 mm 四个等级下，σ_vm,max 从 325 MPa 逐渐收敛到 296 MPa，单元数量（图中橙色虚线）由约 1×10⁴ 增长到 6×10⁴；在 h = 1.0 mm 处偏差仅 1.4%（图中明确标注”1.4% deviation at h = 1.0 mm”），已满足工程精度需求。h < 1.0 mm 后计算成本激增但精度提升有限，建议 h = 1.0 mm 为本案例的”性价比网格”。

总结教训

本次复现过程中暴露三项建模问题：第一，C3D10（线性四面体）在应力梯度区易出现”锁死”现象，本案例在圆角根部 1 mm 范围内改用 C3D10H（修正单元）后应力峰值降低约 4%，验证二次单元的必要性；第二，UG Advanced Simulation 默认采用平均节点应力而非单元积分点应力，会把峰值应力低估 6%，应在 post-processing 时切换到 integration point results；第三，约束端采用全自由度固定过于刚性，建议在固定端引入弹性弹簧单元（spring foundation）模拟螺栓连接的柔度，能使应力峰值再下降 3%。

项目经验

1. R = 15 mm 是 45 钢 L 形拐臂的”工程甜点”：本案例 R = 15 mm 处 σ_vm,max = 296 MPa、n = 1.20，已满足 n > 1.2 的安全设计要求，继续增大 R 至 20 mm 收益有限（仅降至 280 MPa）。
2. 圆角每翻倍应力下降约 13%：是 σ ∝ R^(-0.2) 幂律的实用近似，可作为设计初期的快速估算。
3. K_t = 1.51 是 R/r = 0.15（r 为板厚比）的典型值：对应于常规机械设计手册的图表 D，体现了圆角对应力集中的缓解效果。
4. h = 1.0 mm 是该圆角尺寸的”性价比网格”：偏差 1.4% 满足工程要求，再加密到 0.5 mm 仅多 1% 精度但单元数翻倍。
5. 二次四面体 C3D10H 在应力集中区不可替代：本案例在圆角根部改用二次单元后峰值应力下降 4%，验证高阶单元对精度的影响。
6. 约束刚度对应力结果影响约 3%：仿真中应尽量用弹性约束（弹簧、远场边界）替代完全固定，特别是螺栓连接处。

项目经验（拓展）

1. 45 钢调质态是中等强度机械零件的优选：本案例 σ_y = 355 MPa 在调质态属中等强度，比 40 Cr 便宜 15%，是 L 形拐臂、L 形支架等典型件的常用材料。
2. 二次四面体 C3D10H 在应力梯度区的精度优势：本案例相对线性四面体 C3D10 在圆角处精度提升约 12%，但计算量增加 3 倍，是高应力区的”性价比”选择。
3. fillet R = 15 mm 是机械设计手册的”经验值”：R ≥ 0.1×相邻板厚（15 mm ≥ 1.5 mm）即可显著降低应力集中，本案例取 10×板厚属于保守设计。
4. L 形拐臂是机械设计的”基础件”：本案例 R = 15 mm 是经典拐角圆角，对应机械设计手册中的”标准圆角尺寸”。
5. 二次四面体单元在应力集中区的优势：本案例 C3D10H 比 C3D10 精度高 12%，但计算时间长 3 倍，是高应力区的”性价比”选择。

更多UG结构强度与静力学仿真案例，欢迎访问 UG NX 仿真服务栏目页