1. 引言
重载铁路运输是铁路货运的发展方向。通过多年来对列车重载技术的开发及发展,有效地提高了重载货运运输能力,取得了较大的经济效益。列车重载技术在助力铁路建设的同时也不断带来新的挑战。新挑战、新问题制约着铁路货运朝着“重载与提速”的方向发展。作为当前重载列车主要代表的2万t列车面临着列车编组长、开行密度大、线路复杂以及循环制动缓解后加速缓解作用不明显等问题,这些问题使得列车在循环制动时产生较大的纵向冲动,加剧钩缓装置的磨耗损伤,甚至导致列车断钩、脱轨等重大安全事故的发生,因此列车循环制动时制动减压量的选取尤为重要,研究2万t列车制动减压量对列车纵向动力学性能影响也就成为了必然选择。
探究2万t列车制动减压量可以采用列车纵向冲动试验和仿真分析两种方式。列车纵向冲动试验耗资大、试验周期长、列车纵向冲动规律研究效果较差。与之相反,仿真分析的优点是周期短、耗资少、分析效果好、还可以与试验形成互补并能模拟某些特殊工况加以比较研究。
Janarthanan建立了重载列车纵向动力学瞬态分析仿真模型,并使用该模型对重载列车在不同轨面上的加速度和制动性能进行了仿真研究 [1];澳大利亚的Cole [2] 教授带领其团队研发了列车运行模拟系统,该系统被用于确保列车的运行安全。与之相比国内的纵向冲动仿真工作也在不断发展。杨亮亮等通过Matlab/Simulink对2万t重载组合列车在紧急制动工况、不同编组情况、不同车辆参数及不同运行条件对于列车纵向冲动的影响进行了仿真比较分析 [3],但是其工况的分析单一不具备代表性;魏伟等通过对制动工况下列车纵向冲动作用机理的研究,发现列车纵向冲动是由冲击和挤压作用二者共同造成的 [4];张东芹分析了万吨列车缓解工况下列车管减压量、缓解初速度对列车纵向冲动的影响 [5];魏伟、赵连刚对1 + 2 + 1编组两万吨列车制动特性进行了仿真分析,并得到了2万t列车车钩力的分布特性 [6]。目前有关2万t列车制动减压量对列车制动特性、缓解特性以及列车车钩力等列车纵向动力学性能的研究较少。因此本文采用列车空气制动系统与纵向动力学联合仿真系统(TABLDSS),对不同制动减压量下列车制动及缓解的同步性及长大坡道列车一次制动缓解产生的最大拉钩力进行比较分析,进而探究制动减压量对列车纵向动力学性能的影响,相关研究结果可为2万t列车在长大坡道制动时减压量的选取提供参考,从而减小列车在长大坡道循环制动时产生较大的纵向冲动,避免列车断钩、脱轨等重大安全事故的发生。
2. 仿真模型的建立
2.1. 列车纵向动力学模型
列车由一系列机车与车辆组成,车辆(或机车)间由车钩和缓冲器连接。建立模型时将每个车辆(或机车)视为一个集中质量,车辆间由弹簧阻尼单元连接。其中任意一个车辆(或机车)的i受力如图1所示。其中
、
、
为第i辆车瞬时位置、速度和重力;
、
、
、
、
、
、
分别为第i辆车的瞬时惯性力、车钩力、牵引力或动力制动力、运行阻力、空气制动力、曲线阻力、坡道阻力。
![](//html.hanspub.org/file/3-2760509x19_hanspub.png)
Figure 1. The force diagram of a single vehicle
图1. 单个车辆受力图
对一节车辆i,可以列出运动方程:
(1)
对于首、尾车辆
(2)
式中:n为列车中机车车辆数总和。
通过求解方程(1) (2)组成的方程组可计算出各个车辆的加速度,再利用加速度与速度、位移间的关系,便可得出各个车辆的速度、位移,进而获得各车辆在牵引与制动过程中所有运动参数。
2.2. 空气制动系统模型
计算管路内气体流动是建立空气制动系统模型同时保证模型准确性的关键。列车空气制动系统内各管路均采用细长管,因此气体在列车管内流动时管内气体径向流动效应与纵向流动效应相比可忽略不计。假设管路内壁是刚性的,管路内为一维、非等熵不定常的气体流动,在不计空气重力的情况下,根据气体流动的连续性、动量守恒定律和能量守恒定律建立描述气体状态的方程组 [7] 如下:
(3)
式中:
、u、p、a、k、D、F、f、q、x、t分别为气体密度、流速、压力、音速、比热比、管路直径、管截面积、管内壁摩擦系数、传热率、距离和时间。上述方程组可应用特征线法求解,具体求解过程可参阅文献 [8]。
3. 列车空气制动与纵向动力学联合仿真系统
本文运用列车空气制动与纵向动力学联合仿真系统(TABLDSS),该系统实现了列车制动系统与纵向动力学的同步仿真。通过不断地分析与调试测试,证明了该系统具备较高程度的准确性 [9]。该系统能够模拟司机根据列车所行驶的线路状况、列车编组情况、运行速度等多种因素控制列车安全平稳运行的过程,并且具有良好的制动特性、车钩加速度等参数的同步技术、数值显示及图形显示功能,为分析列车纵向冲动机理,优化制动、缓冲器系统,模拟列车驾驶提供了基本理论和实用工具 [10]。该仿真系统于2017年参加了VSD组织的国际纵向动力学评测,在计算精度和计算速度上获得优异成绩 [11]。
4. 制动减压量对列车制动同步性的影响
本节以2万t列车为计算对象,列车编组为1HXD1+108辆C80车+1HXD1+108辆C80车+可控列尾,牵引总重21,600吨,列车初速度均为70 km/h,计算工况分别为大闸减压50 kPa、60 kPa及70 kPa。线路形式为平直线路。
列车制动同步性是影响列车纵向动力学性能的关键因素。由于种种因素的限制,列车各车辆的制动是不同步的,同时这种不同步性也是列车纵向冲动产生的主要根源之一,各车辆制动的不同步主要表现在两方面:1) 制动波传播受其有限的制动波速的限制,造成各车辆制动时间不一致,靠近机车的车辆相较于远离机车的车辆更早制动。2) 制动缸上升速率不同造成的制动缸升压曲线形状的差异,这种差异是由于不同位置列车管减压速度不同造成的,靠近机车的车辆列车管减压快,制动缸压强上升快,远离机车的车辆正好相反 [12]。两种因素叠加造成了制动不同步的产生。为探究制动减压量与列车制动同步性间的关系,现使2万t列车以相同初速度制动待首尾列车管压力稳定后,比较全列车的制动缸活塞开始伸出时间,10个车输出一个制动缸压强,其比较结果如图2所示。
如图2所示,从制动缸活塞开始伸出分布规律看,三种工况下列车制动波传播规律均是主控机车由列车前部向后部传播,从控机车向列车前后两个方向传播,其中减压60 kPa与70 kPa下部分车辆的制动缸活塞开始伸出时间略早于减压50 kPa时对应时间,但二者总体差异不明显,说明三种工况列车制动波速差别不大。三种工况下最后发生制动的是列车中后部车辆,并非列车尾部车辆。原因是可控列尾装置的排气作用加快了尾部车辆列车管减压速度,对应尾部车辆制动缸升压速度加快,导致尾部车辆相较于中后部车辆提前制动。从首车与最后制动车辆之间的时间差(最大时间差)来看,三种工况下最大时间差相同,差值均为6.9 s。以上情况表明三种工况下列车制动的同步性差异不大,说明常用制动下制动减压量对于列车制动同步性的影响较小。
![](//html.hanspub.org/file/3-2760509x26_hanspub.png)
Figure 2. The curve of the train’s braking characteristic under three operation conditions
图2. 三种工况下列车制动特性曲线
5. 制动减压量对列车缓解同步性的影响
制动减压量主要通过120-1阀加速缓解作用的强弱进而影响列车缓解的同步性。随着列车加长,列车管总容积增大,列车缓解时列车管升压速度较慢,列车的缓解性能也越来越差。为改善缓解性能,仿照供制动用的副风缸模式和局部减压模式,给每个车辆配备了一个加速缓解风缸,实现加速缓解。2万t重载列车普遍应用的120-1阀便采用了加速缓解阀和加速缓解风缸的结构形式,其作用是当120-1阀主活塞动作下移到缓解位,让准备排入大气的制动缸压力空气作为压力信号先引到加速缓解膜板处、使加速缓解阀产生动作,让本车加速缓解风缸的压力空气通过加速缓解阀中被顶开的夹心阀充入列车管。这样,列车管除了有来自机车供风系统的压力空气充入以外,还有来自本车加速缓解风缸的压力空气充入,即加速缓解作用 [13]。加速缓解作用使列车管增压速度不会沿列车长度严重衰减,使列车前后缓解作用的时间差得以减小。制动减压量的增加,势必导致加速缓解作用的变化,进而影响列车缓解的同步性。为探明列车缓解同步性受制动减压量的影响程度,现使2万t列车以相同初速度制动待首尾列车管压力稳定后以相同速度缓解,比较全列车的开始缓解时间,每个车输出一个制动缸压强,其比较结果如图3所示:
![](//html.hanspub.org/file/3-2760509x27_hanspub.png)
Figure 3. The curve of the train’s releasing characteristic under three operation conditions
图3. 三种工况下列车缓解特性曲线
由图3可知,从列车开始缓解时间分布规律看,减压50 kPa、60 kPa、70 kPa下列车缓解时主控机车由列车前部向后充风,从控机车分别向列车前后部两个方向充风,导致前部120辆左右的车辆与后部剩余车辆之间存在较大的缓解时差,其中前120辆车在8 s内全部开始缓解,而后部剩余车辆分别在25 s、21 s、18 s内产生缓解作用,后部120辆车开始缓解平均时长分别为15 s、13 s、11 s。从三种工况缓解波传播曲线比较看,减压50 kPa、60 kPa下列车缓解波速的传播是不均匀的,具有随机性。同时减压50 kPa缓解曲线的斜率最大,即缓解波速最慢。减压70 kPa缓解曲线的斜率最小,即缓解波速最快。从全列车缓解动作顺序看,三种工况下最后发生缓解作用的均是是列车中后部车辆并非列车尾部车辆。由图3可知减压50 kPa下列车加速缓解作用不明显,导致存在中后部部分车辆开始缓解时间相对周围车辆严重滞后的现象(缓解异常现象)。从首车和最后缓解车辆之间的时间差(最大时间差)看,减压60 kPa、70 kPa下缓解传递时间相比于减压50 kPa分别减少16.5%和28.3%。且减压60 kPa下缓解异常现象有所改善,减压70 kPa下缓解异常现象完全消失。说明相较于减压50 kPa,减压60 kPa、70 kPa下列车加速缓解作用逐渐增强,列车缓解波速逐渐加快,列车缓解的同步性逐步提高。故相对减压50 kPa而言减压60 kPa、70 kPa下列车缓解的同步性有所提高,说明常用制动下减压量的增大对列车缓解同步性的增强效果较明显。
6. 制动减压量对长大坡道列车最大拉钩力的影响
2万t列车在长大坡道循环制动缓解后会由于种种因素不可避免地产生较大的拉钩力。列车一次制动缓解后产生的最大拉钩力与列车制动减压量的变化息息相关。本研究以某重载线路长大下坡道为例,在50 kPa、60 kPa、70 kPa三种减压工况下对2万t列车以相同速度制动和缓解的过程进行仿真,比较不同减压工况下列车最大拉钩力,进而探究不同制动减压量对于2万t列车在长大坡道缓解后产生的最大拉钩力的影响。
图4所示线路是某重载线路K16+589至K23+500的一段长大下坡道。图中绘有三种工况下列车的速度仿真曲线、坡道变化及操纵指令。由图可知在仅改变制动减压量的前提下三种工况下列车缓解时坡道的坡度均相同且整个过程中列车再生力保持不变。三种工况下列车均于K20+648处以64 km/h的速度开始制动并以40 km/h的速度缓解,整个过程各车辆最大拉钩力分布情况如图5所示。
![](//html.hanspub.org/file/3-2760509x28_hanspub.png)
Figure 4. The curve of the locomotive operating instructions, ramp changes and speed change under three operation conditions
图4. 三种工况下列车操纵指令、坡道变化及速度变化曲线
![](//html.hanspub.org/file/3-2760509x29_hanspub.png)
Figure 5. The distribution of the maximum coupler force under three operation conditions
图5. 三种工况下列车最大拉钩力分布情况
图5为三种减压工况下列车的最大车钩力沿车长分布曲线。由图所示三种工况下整列车最大拉钩力均发生于109车。减压50 kPa、60 kPa、70 kPa三者列车最大拉钩力分别为1536 kN、1855 kN、1963 kN。相比减压50 kPa,减压60 kPa、70 kPa最大拉钩力分别增大了20.8%和27.8%。说明长大坡道列车最大拉钩力会随着制动减压量的增大而显著增大,即制动减压量对于长大坡道列车最大拉钩力影响较大。
为探究2万t列车最大拉钩力随制动减压量增大的根本原因,本节整理了列车最大拉钩力的主要影响因素并展开研究,其过程如下:缓解时列车的最大拉钩力的主要影响因素有坡道种类、机车电制动力、列车管减压量、列车缓解初速度、车钩间隙等,本研究采取控制变量的方式仅改变列车管减压量而控制其他因素不变,因此影响最大拉钩力的主要因素是列车管减压量。列车管减压量影响着各车辆制动缸最终平衡压强值,后者又与列车的空气制动力密切相关。相邻车辆间的制动力差恰恰是造成车辆间产生车钩力的根本原因 [14]。2万t列车缓解时沿列车长度方向制动力缓解的差异是造成列车纵向冲动的根源,而这种差异又与制动缸降压特性以及列车缓解波速息息相关。因此分别对三种工况下列车缓解波速以及制动缸降压特性二者对列车最大车钩力的影响进行探究,进而分析最大拉钩力增大的原因。
1) 列车缓解波速
由上节可知相较于减压50 kPa,减压60 kPa、70 kPa下列车缓解波速逐渐加快,列车缓解的同步性逐步提高,应该有助于最大拉钩力的降低,但减压60 kPa与70 kPa下整车最大拉钩力反而增大。说明列车缓解波速并非是影响列车最大拉钩力主要因素。
2) 制动缸降压特性
同一编组列车,其制动缸降压特性会随着机车减压量的改变而产生差异。制动缸降压特性越陡,越会加剧列车的纵向冲动。制动缸降压特性变缓,将会有助于减少列车的纵向冲动 [15]。故对三种减压工况下制动缸降压特性的比较如下表1~3所示。
![](Images/Table_Tmp.jpg)
Table 1. Decompression characteristics of the Brake cylinder under 50 kPa braking reduction
表1. 减压50 kPa制动缸降压特性
![](Images/Table_Tmp.jpg)
Table 2. Decompression characteristics of the Brake cylinder under 60 kPa braking reduction
表2. 减压60 kPa制动缸降压特性
![](Images/Table_Tmp.jpg)
Table 3. Decompression characteristics of the Brake cylinder under 70 kPa braking reduction
表3. 减压70 kPa制动缸降压特性
由上表计算可知减压50 kPa、60 kPa、70 kPa三者制动缸平均降压速率分别为8.7、10、11.1。说明相较于减压50 kPa,减压60 kPa、70 kPa制动缸压力下降速度加快,制动缸的降压特性变陡,因此加剧了列车纵向冲动,导致最大拉钩力增大。同时也证明了制动缸降压特性对列车纵向冲动的影响比列车缓解波速的影响显著。
7. 结论
1) 制动工况下减压50 kPa、60 kPa、70 kPa三种工况列车制动的同步性差异不明显,说明常用制动下制动减压量对列车制动的同步性影响较小。
2) 缓解工况下减压60 kPa、70 kPa相比减压50 kPa,120-1阀加速缓解作用增强,列车缓解波速加快,缓解的同步性也有所提高,说明常用制动下列车缓解同步性随制动减压量的增大而提升且影响程度较大。
3) 减压50 kPa到减压70 kPa,列车缓解波速加快,缓解同步性提高,降低了纵向冲动,而制动缸降压特性的变陡,增大了纵向冲动。而最终相对减压50 kPa,减压70 kPa纵向冲动增大,证明制动缸降压特性对列车纵向冲动的影响比列车缓解波速的影响显著。
4) 制动减压量由50 kPa增至70 kPa,制动缸降压特性逐渐变陡,加剧了列车纵向冲动,导致缓解后列车最大拉钩力逐渐增大。因此在保证速度的前提下应尽量采用较小的制动减压量以降低长大下坡道循环制动缓解后列车最大拉钩力。