制动盘热散逸与致阻力研究

2021-05-13 20:42陈磊
科技风 2021年13期
关键词:空气阻力流场温度场

陈磊

摘 要:制动盘为制动摩擦副中承担热负荷的重要部件。在制动工况时期望制动盘通风能力越强越好,在非制动工况时制动盘通风现象会引起空气阻力,消耗牵引功率。研究了制动盘热散逸和致阻力的计算方法,基于叶片型散热筋结构制动盘进行了仿真计算,得到了制动盘热散逸和致阻力的规律,为进一步研究制动盘热散逸和致阻力均衡匹配研究提供研究路径。

关键词:制动盘;热散逸;空气阻力;流场;温度场

中图分类号:U270.35

Abstract:Disc is the key component bearing the thermal load which is expected for a better ventilating ability during braking.Whereas,the ventilation causes air resistance and consumes traction power in other time.The calculation method of thermal dissipation and air resistance for the disc are researched in this study.A simulation was conducted based on the vanetype brake disc and the conclusions about the thermal dissipation and air resistance were obtained,which provide an approach for the investigation of the match between the thermal dissipation and air resistance.

Keywords:brake disc;thermal dissipation;air resistance;flow field;temperature field

任何載运工具都离不开制动系统,对轨道车辆而言,制动系统最重要的使命是要确保安全。随着技术进步,轨道车辆采用的制动方式越来越丰富,从传统纯机械驱动的踏面制动、盘形制动到越来越依赖电能(机)的电阻制动、再生制动、磁轨制动、涡流制动等[1]。这些涉“电”制动方式往往由于其本身的特点,使得它们不能或难以成为轨道车辆的安全制动方式。时至今日,踏面制动或盘形制动这种摩擦制动仍是被普遍接受的轨道车辆安全制动方式,制动时通过摩擦副摩擦生热将列车动能转化为热能,实现列车能量的转移,产生列车制动减速作用[2]。摩擦副在摩擦生热过程中会产生热累积,热累积到一定程度会产生热疲劳,致使制动盘失效[3][6]。目前大多采用带有通风能力的加强筋制动盘,制动盘处于一定流速的流场环境中,空气流动会带走一部分制动热量[7][8]。因此在制动时期望制动盘通风能力越强越好,关于制动盘通风能力散热筋的设计与分析已经有学者进行了研究[9][13]。列车非制动工况时制动盘通风现象仍然存在,必定会引起空气阻力,消耗牵引功率[14][15]。这就造成了制动工况和非制动工况对制动盘通风能力的要求不同。基于此,本文研究了制动盘热散逸和致阻力的计算方法,为进一步研究制动盘热散逸和致阻力均衡匹配研究提供研究路径。

1 制动盘热逸散和致阻力的研究方法

高速流场环境对制动提出了制动温度低,散热效果好的要求。然而在提高制动盘热逸散效果的同时,也会相应地增加了非制动工况制动盘产生的阻力。表征制动盘的热累积效应可采用制动温度,表征制动盘的热逸散效应可采用对流散热系数或者散热功率。制动盘产生的阻力采用泵风功率来表征。本论文旨在通过计算和仿真来研究高速流场作用下制动盘的热累积和热逸散效应,图1是研究的技术路线。基于给定的高速制动盘模型,分别建立制动盘的流场计算模型和有限元热容计算模型,通过计算得出泵风功耗、散热功率和制动温度。得出高速流场作用下制动盘的热累积和致阻力效应规律。

2 制动盘泵风效应及致阻力仿真计算

以叶片型散热筋结构制动盘(φ640×110mm)为对象,研究其制动工况下泵风功率和散热功率的变化规律。其速度随时间变化规律为:

2.1 制动盘模型和流场模型

制动盘参数如下表所示。制动盘三维模型如图2所示。制动盘散热筋结构如图3所示。制动盘流场模型如图4所示。

2.2 计算方法与结果

2.2.1 对流散热系数

制动过程中,对流换热系数随时间的变化曲线如图5所示,t=43s时刻制动盘表面对流换热系数如图6所示。

2.2.2 散热功率

将制动盘表面对流换热系数变化曲线代入ABAQUS进行制动过程中制动盘温度场计算,得到紧急制动下制动盘温度随时间变化曲线如图8所示。制动盘整体散热功率如图9所示。t=43s时刻温度场如图10所示。

2.2.3 泵风功率

在紧急制动的过程中,制动盘结构受到的风阻力与速度成二次关系,泵风功率与速度呈三次关系,制动盘泵风功率随速度变化的关系如图11所示。

3 结语

本文研究制动盘热散逸和致阻力的计算方法,以叶片型散热筋结构制动盘进行了仿真计算,总结如下:(1)随着制动过程的发展,对流换热系数值逐渐减小,其CFD计算值与经验公式值误差小于10%;(2)随着速度逐渐减小,制动盘表面温度逐渐增加;(3)随着速度逐渐减小,泵风功率逐渐减小;散热功率呈现出先上升后下降的趋势。

参考文献:

[1]饶忠.列车制动[M].第2版.北京:中国铁道出版社,2010:912.

[2]夏寅荪.速度达350km/h及以上的高速列车非黏着制动装置的研究[J].城市轨道交通研究,2016,19(10):36.

[3]王文静,谢基龙,李强,等.铁路列车制动盘常用材料的热疲劳性能研究[J].机械工程材料,2005,29(02):4058.

[4]杨月,谢基龙.高速客车SiCp/A356铝基制动盘材料的热疲劳裂纹形成与扩展试验研究[J].铁道学报,2007(05):4347.

[5]Li ZQ,Han JM,Yang ZY,et al.The effect of braking energy on the fatigue crack propagation in railway brake discs[J].Engineering Failure Analysis,2014,44:272284.

[6]Wu SC,Zhang SQ,Xu ZW.Thermal crack growthbased fatigue life prediction due to braking for a highspeed railway brake disc[J].International Journal of Fatigue,2016,87(Jun.):359369.

[7]金星,張永恒,王良璧,等.高速列车制动盘表面对流传热特性[J].科学通报,2015,60(23):22452252.

[8]任慧芳.高速列车制动盘传热特性及热负荷的数值研究[D].兰州交通大学,2019.

[9]Tirovic' M.Energy thrift and improved performance achieved through novel railway brake discs[J].Applied Energy,2009,86(3):317324.

[10]GalindoLpez CH,Tirovic' M.Maximising heat dissipation from ventilated wheelhubmounted railway brake discs[J].Proceedings of the Institution of Mechanical Engineers,Part F:Journal of Rail and Rapid Transit,2012,227(3):269285.

[11]Rajagopal TKR,Ramachandran R,James M,et al.Numerical investigation of fluid flow and heat transfer characteristics on the aerodynamics of ventilated disc brake rotor using CFD[J].Thermal Science,2014,18(2):667675.

[12]潘利科,韩建民,李志强,等.列车制动盘通风散热的数值仿真[J].北京交通大学学报,2015,39(01):118124.

[13]黄晨.高速列车扇形柱直肋制动盘内部通道传热特性的实验研究[D].兰州交通大学,2018.

[14]左建勇,罗卓军.高速列车制动盘泵风效应分析[J].交通运输工程学报,2014,14(02):3440.

[15]Luo ZJ,Zuo JY.Conjugate heat transfer study on a ventilated disc of highspeed trains during braking[J].Journal of Mechanical Science and Technology,2014,28(5):18871897.

猜你喜欢
空气阻力流场温度场
基于并联双离心送风机内外流场的数值模拟研究
基于机器学习的双椭圆柱绕流场预测
能源桩群温度场分布特征数值仿真研究
“金雕”展翅
真实流场中换热管流体诱导振动特性研究
不怕摔的蚂蚁
降落伞
基于内外流场仿真分析对整车性能研究和优化
复合式空调系统的数值模拟研究
“牛顿第一定律”练习