货车车轮踏面损伤有限元分析

货车车轮踏面损伤温度场与应力场的有限元分析货车车轮踏面损伤温度场与应力场的有限元分析摘要：本文通过对 21t 轴重、120km/h 的货车车轮在一次紧急制动过程中的温度场和应力场分布进行有限元模拟，探讨分析了温度场和应力场分布与货车车轮货车车轮踏面损伤有限元分析- 2 -踏面损伤的关系，为车轮的热疲劳损伤机理研究提供了技术储备和参考。关键词：货车车轮，温度场，应力场，有限元模拟，热疲劳损伤。正文：提速和重载是提高铁路运输能力的有效措施，已成为铁路货车发展的趋势。我国货车目前制动方式仍然是踏面制动，列车车轮在强摩擦、高热负荷以及大轮轨作用力等恶劣条件下工作。列车在制动过程中，动能逐渐转变为制动装置产生的热能，对于采用踏面制动的高速重载铁路货车，这样的制动过程非常严苛，由此产生的热疲劳损伤已成为车轮失效的主要形式之一。车轮经过多次制动后，会在车轮与铁轨的接触踏面上产生均匀分布的横向裂纹，周围会伴随剥离、掉块等现象。因此，在国家倡导货运列车提速的前提下，现有的踏面制动正面临的严峻的挑战，也对车轮的抗热损伤能力和疲劳寿命提出了更高的要求。由于热损伤和疲劳损伤都与车轮在紧急制动过程中的温度场和应力场分布有密切的关系，本文以 21t 轴重、120km/h 的货运列车车轮为研究对象，拟结合具体货车车轮的结构，利用建模软件对其建模，通过有限元模拟其紧急制动过程中的温度场和应力场分布，并针对实际踏面损伤情况对其模拟准确性给予评估，为进一步研究车轮的热疲劳损伤提供技术参考。1、车轮紧急制动温度场模拟货车车轮的轮径为 840mm，轮辋内侧内径为 710mm，轮毂孔直径为170mm，轮辋外径为 273mm，理论重量 351 kg。车轮材料为 CL60，材料各项物热参数如下：弹性模量 E =2.05105 MPa，泊松比 =0.3，密度 =7800 kg/m3，热膨胀系数 =10.3 10-6 -1，比热容 c=470 J/(kgK)，热传导率 k =51W/(mK)，对流换热系数 h=40W/(m2K)。由于车轮是周向对称的结构，在考虑热流输入车轮踏面和车轮的热耗散时，可以认为在车轮的周向是无变化的，即温度场是轴对称的，因此，选取车轮的1/18 进行分析，即取周向 20的模型。车轮的三维模型如图 1。图 1 车轮的三维模型货车车轮踏面损伤有限元分析- 3 -1）热流密度的确定：为简化问题，可以认为在高速行进过程中，踏面的温度在周向是均匀分布的，且热流输入也是均匀的。初速度为 120km/h 时，各项制动参数如表 1 所示。表 1 制动参数制动初速度（km/h）轴重（t）减速度（m/s2）制动距离（ m）闸瓦压力（ kN） 摩擦系数热量分配系数120 21 0.556 1000 21 0.278 0.91根据热流密度计算公式（1-1）fdStptq)()(其中： 闸瓦在踏面上扫过的面积； 输入到车轮的热量分配系fS数； 热生成功率。)(tpd计算得热流密度 1)/(83.1459)()( 2mkWtStptqfd考虑到车轮周围的空气流动状态较为简单，本模拟中，取对流换热系数=40W/(m 2K)，并忽略车轮的热辐射影响。本次计算对车轮进行一次紧急制动，并在制动结束后空冷 10 分钟。2）模型的网格划分、约束、对流换热面的确定将在 SolidWorks 中建好的模型导入 ANSYS 软件，进行网格划分。在热分析中选取 Solid90 单元，应力分析选择 Solid186 单元，网格大小 8mm，踏面附货车车轮踏面损伤有限元分析- 4 -近的网格为 6mm，网格划分的结果如图 2 所示。对于热分析，选取车轮的踏面为热流输入面，选取与空气接触的面，设置为对流换热面。如图 3。1XYZFile: F:Whel01whel.x_t MAY 12 201016:51:1AREASTYPE NUM1XYZFile: F:Whel01whel.x_t MAY 12 201017:02:12ELEMENTS图 2 车轮模型的网格划分图 3 热流输入面（左）和对流换热面（右）货车车轮踏面损伤有限元分析- 5 -在制动盘的剖切面上施加周向对称耦合，如图 4。对模型施加载荷，输入载荷命令流，设置输出每个子步选项，并开始计算。3）计算结果及分析图 5、图 6 分别显示了制动结束时刻和制动结束空冷 10 分钟后的温度场分布。制动结束时刻的最高温度为 224，冷却 10 分钟后，最高温度为 97.5。图 4 在剖切面施加周向对称耦合图 5 制动结束时刻(60s)的温度场分布货车车轮踏面损伤有限元分析- 6 -分别取踏面上如图 7 所示的位置的节点，读取其制动过程中的最高温度，并绘制其制动过程中的节点温度变化曲线，如图 8 所示。各点在制动过程中的最高温度如表 2 所示：表 2 踏面上不同位置节点在制动过程中的最高温度节点编号 1 2 3 4 5温度（） 201 308 315 262 87图 6 制动完毕冷却 10 分钟后的温度场分布图 7 在车轮踏面上选取生成曲线的节点51243货车车轮踏面损伤有限元分析- 7 -3分别取车轮径向上如图 9 所示的位置的节点，读取其制动过程中的最高温度，并绘制其制动过程中的节点温度变化曲线，如图 10 所示。12453图 8 踏面不同位置节点在制动过程中的温度变化曲线12345图 9 在车轮径向上选取生成曲线的节点货车车轮踏面损伤有限元分析- 8 -各点在制动过程中的最高温度如表 3 所示：表 3 径向上不同位置节点在制动过程中的最高温度节点编号 1 2 3 4 5温度（） 315 167 113 110 91综合以上的模拟结果，经过分析后可以得到如下结论：（1） 在整个制动过程中，温度的最高点都集中在闸瓦与车轮的接触摩擦面部位；（2） 随着制动过程的进行，货车的运行速度也在下降，闸瓦与车轮踏面的摩擦剧烈程度也发生着变化，导致制动过程的温度不是一直上升，而是达到峰值后缓缓降低；（3） 踏面上，只有摩擦面与闸瓦接触的部位温度较高，其余部分的温度较低；（4） 车轮的温度是由踏面向轮轴位置逐渐降低的，越靠近轮轴位置，轮的温度越低，出现最高温度的时间也越晚。12345图 10 径向不同位置节点在制动过程中的温度变化曲线货车车轮踏面损伤有限元分析- 9 -2、车轮紧急制动应力场模拟1）输入载荷的确定对于货车车轮在紧急制动过程中应力场的模拟，可以利用间接偶合法，将紧急制动全程车轮的温度场的模拟结果作为应力场的载荷，施加在模型的所有节点上，利用命令流，调用温度场计算结果文件（*.rth 文件）中的节点温度值到结构分析中，并在结构分析时设置与热分析相对应的每一个载荷步，以确保耦合的正确性。2）结构约束的施加由于车轮与车轴是压装的，属于过盈装配，因此需要在轮毂的内圈施加沿轴向和径向的零位移约束，即全约束。由于模型选取的是周角对称的 1/18 模型，需要在剖切面上施加对称耦合，以满足边界连续的条件。结构约束示意图如图11。3）计算结果及分析图 12、图 13 分别显示了制动结束时刻和制动结束空冷 10 分钟后的应力场分布。制动结束时刻的最大应力为 370MPa，冷却 10 分钟后，最大应力为126MPa。图 14 显示了制动过程中车轮中有最大应力值时刻（制动开始后的第 25 秒）的应力分布云图，最大应力值为 460MPa。图 11 模型结构约束的施加方式货车车轮踏面损伤有限元分析- 10 -分别取如图 15 所示的位置的节点，读取其制动过程中的最大热应力，并绘制其制动过程中的热应力变化曲线，如图 16 所示。图 12 制动结束时刻的应力场分布图 13 制动空冷 10 分钟后的应力场分布图 14 制动过程中出现最大应力时刻的应力场分布货车车轮踏面损伤有限元分析- 11 -综合以上的模拟结果，经过分析后可以得到如下结论：（1） 在整个制动过程中，随着踏面温度的迅速升高，成为车轮的最高温部位，有强烈的热膨胀，因此此过程应力值最大的点都在闸瓦与车轮的接触摩擦面部位；（2） 随着制动过程的进行，车轮踏面的热量逐渐传入车轮内部，使得踏5 4 3 2 1图 15 在车轮径向上选取生成应力曲线的节点1 2345图 16 径向上不同节点在制动过程中的 Von Mises 应力曲线货车车轮踏面损伤有限元分析- 12 -面温度与轮辋和轮毂温度的差值缓缓减小，从而制动结束时，踏面的温度迅速下降，而使内部温度高于踏面，最大应力部位产生在车轮踏面之下；（3） 踏面的最大热应力随着与内部的温度差的减小，也并不是一直递增的，而是达到峰值后缓缓下降；（4） 踏面附近较大的应力为周向的拉应力。3、踏面疲劳损伤分析通过对上述车轮紧急制动过程中的温度场和应力场的分析，可以看出车轮在制动过程中会由于温度变化产生较大的热应力，使得踏面附近极易发生屈服，而车轮原本存在的加工缺陷和热处理留下的参与应力与之相互作用，会促进热裂纹的萌生。在短时间的连续制动情况下，踏面会来不及降温而反复接受摩擦，使得温度不断上升，踏面附近由原本的无应力或残余的周向拉应力变为受到周向的压应力，而制动结束，温度降低后，踏面附近则会恢复到残余的周向拉应力状态，再次制动又会使周向受拉，如此往复循环，会产生低周疲劳，是车轮在服役时最严苛的服役状态，极易产生疲劳裂纹。因此可以推断，在货车的制动过程中，由于其温