高速列车在风雪严寒环境运行时,道床上的积雪易被转向架与地面之间的强剪切流动卷入转向架区域,雪粒经历跟随、粘附和沉积等复杂过程,在转向架表面堆积
[1-2].在转向架发热部件(如制动夹钳、电机等)的作用下,部分积雪融化成水滴,并滴落在旋转的轮对和制动盘上,扩散至转向架各个区域,形成甩水现象
[3].四处飞溅的液滴致使转向架表面湿度增加,雪粒更容易粘附到转向架表面,加剧转向架积雪结冰.转向架结冰不仅增加列车轴重,严重时甚至引起制动系统失效,威胁列车运行安全
[4-5].
针对转向架积雪结冰问题,国内外学者开展了大量研究.Shishido等
[6]研究发现转向架区域的积雪主要来自于列车底部或侧面,且认为轮对旋转对转向架区域的积雪分布具有重要影响.Xie等
[7]基于离散相模型,预测了高速列车转向架区域雪粒的运动轨迹,揭示了雪粒在转向架上的分布情况.Wang等
[8]分析了雪粒进入转向架区域的运动过程和机理,探明了轮对旋转效应对风雪运动的影响.Bae等
[9]通过建立雪粒起跳的跃移模型,模拟了多车编组列车转向架区域的积雪,总结出积雪在第三节车厢后保持稳定的规律.Wang等
[10]采用滑移网格技术,再现了列车通过积雪道床的场景,并分析了列车风卷雪机理.冯永华等
[11]研究了高寒列车转向架区域空气流场特性,总结了形成积雪的可能原因,为转向架防积雪的优化工作提供依据和基础.何德华等
[12]和王东屏等
[13]通过设计流线型扰流板的防冰雪方案,成功减少了转向架区域的雪粒黏附数量,提升了转向架防结冰能力.
以上文献研究了雪粒在转向架区域的运动轨迹和分布,但是没有涉及结冰相变过程.近年来,关于高速列车轮对甩水致转向架结冰的研究才被报道.Liu等
[3]采用欧拉-拉格朗日耦合数值模拟方法研究了轮对甩水对高速列车转向架区域结冰的影响,发现水滴聚集在制动夹钳处,易导致大量结冰.娄振等
[14]采用液膜模型研究制动盘甩出的水滴在转向架表面及转向架舱的沉积,总结了水滴在转向架各区域的分布特性.除了数值模拟,风洞实验也被用于开展转向架积雪结冰领域的研究.Allain等
[15]在法国气候风洞开展了列车车底积雪实验研究,采用1∶2的列车缩比模型,探究了风速、气流温度和雪粒特性等不同参数对积雪结冰增长和分布的影响.Shon等
[16]和Ji等
[17]在维也纳的小型气候风洞开展了比例为 2∶3的转向架积雪实验,并以此验证了数值模拟方法的正确性.Liu等
[18]在风洞中模拟了转向架周围粒子运输过程,评估了转向架底部不同防护方案的效果.Wang等
[19]和Zhang等
[20]在中南大学开展了列车转向架积雪结冰实验,提出防结冰设计,并通过对比实验讨论了其效果;然而,这仅是为检验所设计方案的性能,并没有进行足够的研究来定量分析转向架的结冰特性.
综上所述,仅有少数学者开展了转向架结冰实验研究,高速列车转向架区域的结冰分布特性有待进一步探究,列车长时间运行条件下转向架的结冰速率尚不明确.因此,本文将开展高速列车轮对甩水致转向架结冰风洞实验,分析转向架区域的动态结冰过程,探明转向架区域结冰分布特性,建立转向架结冰速率公式,为转向架结冰质量快速预测提供参考依据.
1 转向架轮对甩水致结冰风洞实验
1.1 积雪结冰实验风洞系统
本研究在中南大学轨道车辆积雪结冰风洞开展
[1].如
图1所示,该风洞是专门为轨道车辆转向架积雪结冰实验所设计的闭路风洞,包括制冷系统、动力段、整流段和实验段.制冷系统可以将空气温度下降至0 ℃以下,并由动力段的风扇高速旋转将冷空气输送到实验段.本次实验在实验段进行,实验段的喷口尺寸为4 m×1 m,最大风速可达100 km/h.在制冷系统的作用下,整个风洞的最低温度可达-20 ℃.实验前对风洞实验段进行了流场标定,并评价了风洞流场品质,达到了《低速风洞和高速风洞流场品质要求》(GJB 1179A―2012)中规定的要求
[21].因此,该风洞系统可再现转向架区域的流场特性,满足本次研究的要求.
1.2 实验模型
本研究采用缩比尺寸为1∶3的某高速列车动力转向架,其基本结构如
图2(a)所示,除了部分线路、管路被简化外,转向架的主要结构均被保留,以准确再现转向架区域流场.此类转向架包含电机和齿轮箱等动力装置.列车在风雪环境下运行时,这些关键零部件区域容易受到风雪侵扰并产生结冰.为实 现转向架轮对的旋转,采用三相变频电机作为动力源,通过皮带轮带动转向架的轮对转动,动力装置如
图2(b)所示.
为了模拟列车运行时转向架周围车体边界,制作了车体模型,如
图3所示.本次实验主要对转向架以及车体底部(转向架舱)的结冰特性进行研究,车体的上部被适当简化
[19-20],但是保留了车体底部以及前、后端板等主要部分,以准确模拟转向架区域的真实流场特性.此外,根据实际轨道参数,还制作了相应的轨道模型.考虑到风洞实验的地面效应,将实验地板进行了一定程度的抬高.
1.3 实验设置
本次实验需在低温环境下开展。根据冬季东北地区高速列车运行条件,实验初始阶段通过制冷系统将整个风洞温度调整为-15 ℃左右,再由动力段风扇驱动系统高速旋转,产生持续恒定的低温气流.在冬季运行时,根据停站的需要,列车存在不同的运行工况,包括加速、匀速、制动、停止等阶段,共同组成一个运行周期.为探究转向架结冰分布特性,该实验对列车运行工况进行适当简化,在整个实验期间,列车轮对的转速保持恒定.风洞实验段喷口风速设置为27.78 m/s,对应100 km/h.调节电机频率值为 83 Hz,使轮对旋转速度与风速相匹配.
为模拟列车制动时制动夹钳周围融水产生的喷水现象以及轮对高速旋转时液滴四处飞溅的情况,设计了如
图4所示的喷水系统.该喷水系统共包含8个喷头,分别安装在4个制动夹钳的两侧.喷水总流量由一端的总阀控制,并通过电子流量计显示读数.根据Wang等
[19]和Zhang等
[20]喷水流量的选取依据,考虑到本文采用的是缩比模型,喷水总流量设计为1.3 L/min.为防止低温环境下喷水管道和喷头结冰导致的喷水流量下降,采用电加热带和保温材料包裹水管,并在整个实验阶段都开启加热功能.在冬季风雪环境下,列车匀速行驶阶段会采取间歇制动以保持制动能力
[19],因此,喷水系统在整个实验过程中都保持工作.
实验准备阶段,将转向架及车体模型安装在风洞实验段中间位置,以确保来流能够均匀地流过实验模型,整体布置如
图5所示.实验温度和实验风速达到实验条件时开始实验,设置电机参数,带动转向架轮对以特定速度旋转,同时启动喷水系统.实验设定运行时间结束后,收集整理实验结果.
1.4 实验数据处理
实验结果分为定性结果和定量结果两部分.定性实验结果包括实验过程中转向架舱和转向架表面动态结冰过程和最终结冰分布,定量实验结果包括转向架舱和转向架表面各区域的结冰质量和结冰速率.为了获得定性实验数据,实验时将两台分辨率为1 920×1 080的防雾高清摄像机布置在转向架两侧,从转向架左右侧面进行录像,确保可靠捕捉到转向架表面完整的结冰过程.每次实验结束后,通过相机拍照来记录转向架区域的结冰分布情况.为准确表征转向架表面的结冰分布特性,分别根据不同部位和不同部件
[20],将转向架舱和转向架表面分别划分为不同区域,统计每个区域结冰质量,如
图6所示.根据转向架各区域的结冰质量,可判断风雪环境运行时转向架结冰的严重程度,为防结冰方案的设计提供参考.
2 重复性实验
为验证实验的可靠性,进行三次重复性实验,每次实验结冰时间设定为40 min.实验结束后,根据
图6划分的区域分别收集并测量转向架舱表面和转向架表面的结冰质量,并进行数据分析.
转向架舱及转向架各个区域的结冰质量及总质量结果如
表1和
表2所示.引入相对实验误差
[19],其计算公式如下:
式中,D表示实验误差,M表示每次实验的结冰质量,Ma表示3次实验结冰质量的平均值.
从
表1和
表2可以看出,3次实验中转向架舱和转向架各个区域的结冰质量呈现出一致的分布趋势.且由于3次实验中总的喷水量相同,3次实验所收集的结冰总质量十分接近,转向架舱、转向架的结冰总质量最大误差分别为1.42%和2.88%.实验结果说明该实验重复性良好,满足本次实验的要求,同时表明了该实验系统可靠且实验设计等各方面具有合理性和可行性.
3 转向架结冰特性
3.1 转向架区域动态结冰过程
本次实验在一定程度上还原了高速列车转向架的真实动态结冰过程.
图7为结冰时间40 min工况下,每间隔8 min转向架区域的结冰特性.整体上,转向架区域的结冰和冰柱随结冰时间增加而迅速发展,不同部位的冰块展现出不同的形状、大小,直至整个转向架区域都被冰层覆盖.
对于转向架舱,在前、后端板处能够清晰观察到形成的块状冰和条带状冰,并且随时间增加结冰形态变得更加明显,如
图7绿色虚线区域所示.前者是由于前轮对甩出的水滴沿着前端板流出,在前方来流作用下形成,后者则是后端板处直接受到后轮对甩出水滴的冲击而形成,这与Wang等
[19]和Zhang等
[20]的研究结果一致.随着结冰时间的增加,紧贴车体底部的冰层逐渐变厚,如
图7蓝色轮廓线所示.以上分析表明列车运行时,轮对高速旋转甩出的水滴在紊乱流场作用下四处飞溅,导致车体底部大量结冰.
对于转向架局部区域,随结冰时间的增加,转向架的构架、轴箱、旋转臂表面覆有层冰.在第24 min以后,空气弹簧表面的结冰量基本不再改变,水滴开始沿着其表面滴落,在低温作用下形成细小的冰柱.如
图7红色虚线区域所示,由于制动夹钳靠近轮对,其结冰情况更为严重,造成大量不规则结冰堆积.
转向架区域的结冰形态在一定程度上反映了车体底部的流场特性。如
图7(d)~
图7(f)紫色虚线区域所示,冰柱向前倾斜,表明该区域形成一个回流,气流方向与来流方向相反.对于转向架底部以及前、后端板,冰柱向后倾斜,说明转向架底部的气流直接从前往后流动.这导致轮对甩出的水滴跟随气流流动,大部分聚集在转向架后部和转向架舱的后端板处.
3.2 转向架结冰分布特性
该小节对转向架区域的最终结冰分布特性,包括转向架舱、转向架各个区域的结冰分布和结冰质量进行详细分析.
图8给出了实验时间为40 min时,转向架舱的最终结冰分布.从
图8(a)可以看出,由于来流的作用,转向架舱前底板没有发生结冰,结冰主要集中在前底板以后的区域.
图8(b)展示了整个转向架舱有大量结冰产生.
图8(c)、
图8(d)显示前、后端板的结冰以块状为主,且后端板的结冰量大于前端板.转向架舱顶板表面呈水平方向,其结冰以柱状为主,水滴沿着来流方向在轮对正对位置两侧形成多个长短不一的冰柱.
图9展示了转向架舱前、后端板的结冰特性.由于正面轮对前、后端板形成了较厚的结冰.在来向气流的作用下,从轮对甩出的水滴沿着端板流出,在前端板底部形成体积较大的冰块,而后端板的水滴则流向后底板,形成多个向后倾斜的冰柱.列车在实际运行过程中,聚集在车底的结冰会逐渐增加,并且靠近轮对.这不仅增加了列车的重量,掉落的冰块还会威胁到车下设备,对列车运行环境造成不利影响
[22].
图10显示了转向架整体以及部分区域的结冰分布特性.构架、制动夹钳和空气弹簧是结冰较为严重的区域.由于构架靠近转向架舱顶板,从顶板流下的水滴不断形成冰柱,并填满构架与顶板之间的间隙,如
图10绿色线框区域所示.制动夹钳靠近电机、齿轮箱和横梁等部件,空间相对狭窄,水滴更容易粘附在其表面并迅速积聚,结冰几乎覆盖了整个制动夹钳,在其底部更是悬挂了多个向后倾斜的冰柱.部分水滴来不及结冰,沿着已经凝结的冰柱不断生长,甚至接触到地面,如
图10红色线框区域所示.当制动夹钳结冰时,其制动效率降低,导致制动失效,直接威胁列车的运行安全.空气弹簧表面覆盖了结冰以及悬挂的冰柱.值得注意的是,在空气弹簧下部及周边均有结冰,如
图10蓝色线框区域所示,这无疑会阻碍空气弹簧的有效行程.除此之外,在齿轮箱和电机的迎风面也观察到大量结冰.
为了定量分析转向架区域的结冰分布特性,收集了转向架舱和转向架各个区域的冰块并测量其质量.转向架舱各区域结冰质量如
图11所示.可以看出:整个转向架舱中,后端板的结冰质量较大,占结冰总质量的28%,这与之前的分析较为符合.其次是前端板的结冰质量,占结冰总质量的18%.转向架舱顶板的结冰质量之和占结冰总质量的37%,这是由于顶板所覆盖的面积最大,从轮对甩出的水滴和跟随流场运动的水滴易粘附到该区域.在顶板区域,顶板后部的结冰质量仍要大于顶板前部,这与转向架舱整体的结冰分布特性一致.
图12给出了转向架各个区域的结冰质量.结果显示:构架和制动夹钳区域的结冰质量最大,分别占转向架结冰总质量的34%和22%.一方面,构架的表面积占比大,容易粘附更多的水滴.另一方面,构架表面存在凹陷处,水滴易在这些区域聚集并不断形成冰块,从而导致构架结冰严重.对于制动夹钳,其正面轮对,从轮对甩出的水滴是其结冰的主要原因.齿轮箱和电机的结冰也较为严重,分别占结冰总质量16%和9%.由于空气弹簧、横梁和纵梁分布在转向架上半部分,水滴难以在这些区域聚集,结冰量相对于其它区域较少.转向架的整体结冰分布特性揭示了列车在运行过程中容易受到冰雪侵蚀的部位,为转向架有针对性防结冰和回库融冰提供了参考.
根据上述定性和定量分析得出如下结冰分布特性:转向架舱呈现出后部结冰大于前部的分布规律;对于转向架,由于水滴不断往底部流动聚集,在制动夹钳、电机和齿轮箱底部形成厚大的冰块和长短不一的冰柱,呈现出底部结冰量大、结冰形状复杂的特点.整个转向架区域的结冰呈现出底部多于上部,后部多于前部的分布特性.
3.3 转向架结冰速率模型
已有关于转向架结冰的实验研究主要集中在单一结冰时间下的结冰分布特性和防结冰装置设计
[19-20],并未考虑结冰速率.列车长时间运行过程中,转向架表面结冰不断增加
[22],可能引起结冰分布特性改变.因此,明确不同结冰时间下转向架表面结冰质量,提出结冰速率预测模型,可为列车转向架针对性防除冰提供参考依据.为此,除了开展结冰时间为40 min的实验外,还分别开展了结冰时间为 30 min、50 min和60 min的甩水实验,以更好地评估转向架的结冰速率特性.
图13给出了转向架各区域的结冰速率.由
图13可知,各个区域的结冰质量随结冰时间的增加而增加,但是结冰速率不同.构架、制动夹钳和横梁呈现出后端结冰质量大于前端结冰质量的特点.前、后端齿轮箱的结冰速率保持一致,且前端齿轮箱的结冰质量大于后端.随着结冰时间的增加,前、后端电机的结冰质量呈现出从基本相等到前端大于后端的规律.由于转向架结构的中心对称,左、右端空气弹簧的结冰质量和结冰速率并不相等.
图14给出了转向架舱、转向架的整体结冰速率模型,并拟合得出了结冰速率公式.转向架舱结冰总质量
y与结冰时间
x之间的拟合公式为
y=0.129 2
x,转向架结冰总质量
y与结冰时间
x之间的拟合公式为
y=0.130 7
x.上述两个公式的决定系数
R2分别为0.995 9和0.998 8,表明该预测曲线的估计值与实验数据之间的拟合程度较高.以上拟合公式表明转向架舱、转向架的结冰质量随结冰时间呈现线性增加的趋势,在本实验结冰时间范围内结冰速率保持恒定.
由于本文结冰质量预测模型是基于60 min内的结冰时间获得,该预测模型对于长运行时间下转向架结冰质量的预测,还存在不足.在后续的研究中,可以进一步开展更长时间下的转向架结冰风洞实验,拓宽本文结冰质量预测模型的适用范围.
4 结 论
本文基于中南大学轨道车辆积雪结冰风洞开展了高速列车轮对甩水实验,分析探明了转向架舱、转向架的结冰分布特性,并拟合了结冰速率曲线,得出如下结论:
1)轮对甩出的水滴在紊乱流场作用下扩散至转向架各个区域并产生结冰,转向架舱呈现出后端结冰质量多于前端的结冰特性,转向架呈现出底部结冰量大、结冰形状复杂的特性.
2)对于转向架舱,其后端板结冰严重,结冰质量占转向架舱结冰总质量的28%;对于转向架,其结冰主要分布在构架、制动夹钳、齿轮箱和电机底部,空气弹簧、横梁和纵梁等区域结冰较少,其中构架和制动夹钳结冰质量占比权重最高,分别为转向架结冰总质量的34%和22%.
3)建立了一定运行时间范围内转向架结冰质量快速预测模型,转向架舱、转向架的结冰总质量与结冰时间呈一次函数关系,其结冰质量随结冰时间线性增长,为转向架结冰质量预测提供了参考依据.
国家重点研发计划资助项目(2020YFA0710904-02); 中国国家铁路集团有限公司系统性重大项目(P2021J036)