火车零部件有限元分析

轮轨接触是铁路系统里最关键的力学问题之一，接触应力集中在很小区域，却决定了车轮和钢轨的寿命。我这个项目做的是 火车零部件有限元分析 里的轮轨接触问题，用 Hertz 接触理论加有限元，算出接触应力分布，并解释最大剪应力为什么出现在表面下几毫米——这正是铁路疲劳裂纹的萌发位置。

理论起点是 Hertz 接触。两个曲面在法向压力下的弹性接触，接触区是个椭圆（轮轨接近二维线接触，退化为窄矩形带），接触应力呈半椭球分布，中心最高、向边缘衰减。我先按 Hertz 公式估算接触斑半宽和最大接触压力 p₀，作为有限元模型的校验基准。

有限元建模只取车轮和钢轨的局部（没必要建整车），车轮半径 460 mm，轴重按 17 t 换算成垂向载荷。接触区域网格加密到 0.2 mm，远离接触区逐渐变粗（兼顾精度和算时）。材料用线弹性钢，E=210 GPa、ν=0.3。边界条件：钢轨底部全约束、车轮施加垂向力。

结果分两部分看。左图 von Mises 应力云图显示，接触斑中心应力最高，约 350 MPa，向两侧快速衰减。右图是接触应力随深度衰减曲线：最大剪应力 τ_max 并不在表面，而是出现在表面下约 5 mm 处——这和 Hertz 理论的经典结论一致，也是铁路工程里”为什么裂纹从皮下起”的根本原因。客户据此在养护中把检测重点放在轨头皮下 5 mm 处，和后续探伤发现的裂纹位置吻合。

判读逻辑：接触应力是”表面压、皮下剪”，表面压应力不直接导致破坏，但皮下最大剪应力超限会萌生剪切裂纹，随后扩展到表面形成剥落（shelling）。这是个很典型的”应力位置 vs 破坏位置不一致”的案例，单看表面应力会误判。

踩坑记录：第一，接触收敛。轮轨接触是高度非线性的，初始接触对设不准会穿透或发散。我先按 Hertz 估算接触斑位置精确设初始间隙，再用增强拉格朗日法（augmented Lagrange）或纯罚函数求解；太松会穿透、太紧不收敛。第二，网格敏感性。接触应力对接触区网格尺寸极敏感，0.2 mm 是反复试出来的，太大（1 mm）会低估峰值 20% 以上。第三，摩擦影响。轮轨间有蠕滑和摩擦，我这个项目为简洁先做了无摩擦接触，但真实牵引工况下摩擦力会改变剪应力分布，后续补了库仑摩擦（μ=0.3）做对比。第四，材料非线性。高轴重下接触区可能局部屈服，我这个项目用线弹性，若要精确须加弹塑性。第五，子模型技术。全局部网格太粗看不清接触，我用子模型（submodeling）把接触区单独加密重算，效率高。

从我的工程经验看，火车零部件有限元分析 里轮轨接触最吃网格和接触算法，Hertz 解析解是必备的”标尺”。我给客户的建议是：用 Hertz 算出的 p₀ 和接触斑尺寸先校准 FEA，再做摩擦/塑性延伸。更多接触仿真，见 [火车零部件有限元分析](https://www.keyanxueshu.com/category/fea/)；需要深入接触分析，[火车零部件有限元分析](https://www.keyanxueshu.com/) 上有入口。