Impact of air flow characteristics underneath carbody on snow accumulation in bogie region of high-speed train
-
摘要: 针对高速列车转向架区域的积雪问题, 建立了包含精细化转向架的列车空气动力学模型; 采用分离涡模拟方法, 对运行速度为350 km·h-1的高速列车周围空气流场进行了模拟, 分析了空气流场特性对车底与转向架区域雪粒输运的影响; 提取了涡核线, 研究了转向架区域的涡流特征与雪粒输运的关系。研究结果表明: 车底气流主要由前后轮对后部向上翻转进入转向架区域, 绕轮轴形成旋转气流; 转向架底部区域涡量大于1 000 s-1, 涡流基本为纵向; 转向架顶部区域涡量小于200 s-1, 涡流基本为纵向; 转向架轮对与前后端墙的空隙处涡流多为竖向, 且后部轮对处的涡量较前部轮对处大5倍以上; 转向架内部区域涡量小于200 s-1, 涡流走向杂乱; 涡流的尺度、强度与走向特性反映出进入转向架区域的气流具有较强的挟带雪粒的能力, 而流出转向架的气流挟带雪粒的能力较弱; 头车下部区域负压较大, 车底与裙板两侧存在强度较大的涡流, 易卷起轨道积雪形成雪烟; 除头车外, 车底与转向架表面绝大部分区域壁面剪切应力小于1 Pa, 对应的摩擦风速小于0.9 m·s-1, 沉积的雪粒不易被内部气流剪切走。Abstract: Aiming at the problem of snow accumulation around the bogie region of high-speed train, the aerodynamics model of the train including refined bogie was established. The detached-eddy simulation was used to simulate the air flow field around the high-speed train at the running speed of 350 km·h-1. The influence of air flow field characteristics on the snow particle transport in the vehicle bottom and bogie region was analyzed. The vortex core lines were extracted, and the relationship between the vortex characteristics of the bogie region and the snow particle transport was studied. Analysis result shows that the underneath air flow is mainly turned upside down from the front and rear wheelsets into the bogie region, forming a swirling air flow around axle. The vorticity in the bottom of bogie region is large than 1 000 s-1, and the vortexes are basically longitudinal. The vorticity in the top of bogie region is less than 200 s-1, and the vortexes are basically longitudinal. The vortexes in the gap between the bogie wheelset and the front/rear end walls are mostly vertical. The vorticity around the back wheelset is more than 5 times larger than that around the front wheelset. The vorticity in the interior region of bogie is less than 200 s-1, and the vortexes are chaotic. The scale, strength and direction characteristics of the vortex reflect that the air flow entering the bogie region has a relative strong ability to carry snow particles, while the outflow of bogie has weaker ability to carry snow particles. The negative pressure in the lower part of head car is large, and there are strong vortexes on both sides of the vehicle bottom and the skirt, which makes it easy to roll up the snow on the track to form snow smoke. Except for head car, the shear stress on the vehicle bottom and the bogie is less than 1 Pa in most areas, and the corresponding frictional wind speed is less than 0.9 m·s-1. The deposited snow particles are difficulty sheared by the internal air flow.
-
0. 引言
近年来,随着中国交通建设的飞速发展,出现了大量的黄土公路隧道。但由于黄土垂直节理发育,强度较低,具有明显的结构性、水敏性和区域差异性,使得隧道开挖引起的扰动变形和后期蠕变变形比较大,并且公路隧道多具有断面大、扁坦状、地质复杂等特点,所以对黄土隧道支护体系的力学特性研究成了亟待解决的问题。尽管中国黄土隧道工程实践已有半个多世纪,但设计理论研究仍然滞后于工程实践,黄土隧道一直是工程界和学术界关注的热点之一[1-2]。
在此背景下,许多学者对黄土的物理性质与黄土隧道的力学特性进行了多年的研究,取得了一定成果。Kruse等结合多个工程案例研究了黄土结构对黄土变形及力学特性的影响[3];Li等对黄土进行了一系列室内土工试验,分析了水含量和环境压力与黄土结构的应力应变关系,并对黄土结构坍塌进行试验评估[4-5];赵占厂等在浅埋黄土隧道衬砌受力测试的基础上,研究了围岩压力的分布形式以及随时间的变化规律,对深浅埋的界定标准也作了相应的论述[6];来弘鹏等以青土岘隧道为依托进行了现场测试,对围岩压力、格栅拱架钢筋轴力、初支和二衬接触压力等变化规律及分布特性进行了研究[7];陈建勋等基于刘家坪2#隧道对黄土隧道洞口段支护的受力状况进行施工监测,并采用有限元法进行了分析[8];杨建民等针对郑州至西安客运专线大断面黄土隧道设计中存在的问题,对深浅埋分界、深浅埋围岩压力进行了研究[9];李鹏飞等在对44座隧道91个监测断面围岩压力统计的基础上,研究了隧道围岩压力的分布特征及其与隧道埋深、围岩级别的关系[10]。
以上针对不同区域与条件下的黄土隧道的力学特性进行了相关研究,为黄土隧道的设计与施工提供了一定技术支撑,但由于黄土的区域差异性明显,多个已建成隧道出现了各种病害,并且在建隧道施工过程中也多次出现塌方现象[11-12],加之黄土自身的特殊性和复杂性,形成时期不同的湿陷性黄土的力学性质也相差较大,因此,有必要对其做进一步的研究。大有山隧道所在西宁市无黄土隧道施工先例,地质条件复杂多变,洞身同时处在湿陷性和非湿陷性黄土地层,地表有杂填土和新近堆积黄土等,在进洞过程中频遇塌方。鉴于此,本文以大有山隧道为依托,采用钢弦式传感器对软弱黄土隧道各典型断面的支护力学特性进行系统测试和分析,为软弱黄土隧道优化设计提供重要基础数据,以期提高黄土隧道的修筑技术水平。
1. 工程概况
大有山隧道是丹东至拉萨西宁过境高速公路的关键性控制工程[13],所属地貌单元主要为黄土塬梁地貌,范围内地层为第四系全新统的黄土、上更新统风积黄土及人工填土,具有压缩性高的特点,湿陷深度达10~20 m。隧道进口段为北川河西岸Ⅲ级阶地地貌,出口段为黄土塬梁深切沟谷斜坡地貌,海拔高程为2314.8~2444.2 m,隧址区山体冲沟发育。
该隧道为分离式双洞隧道,右线起、讫桩号分别为K2+655、K5+190,全长2535 m,左线起、讫桩号分别为ZK2+655、ZK5+215,全长2560 m。围岩以第四系黄土为主,且多具有湿陷性,级别均为Ⅴ级,根据隧道埋深以及上覆土层的不同将结构分为冲沟回填加强段(杂填土段)、深埋段与浅埋段。隧道衬砌设计遵循“强初支,弱二衬”的思想,其主要支护参数见表 1,L为支护长度。
表 1 隧道支护参数Table 1. Parameters of tunnel liner and support2. 测试内容与方案
2.1 测试内容
为研究支护体系的力学特性,现场测试内容包括围岩和初支间的接触压力、初支和二衬间的接触压力、二衬和仰拱混凝土的应变。对现场测试仪器的布设应高度重视,本次受力测试全部采用适宜于土质条件的沥青囊式压力盒,埋设时将其放在预先挖好的深约3 cm的圆坑中,并使沥青囊上面与土面相持平,固定之后对压力盒施加预压力。同时,为了使测试数据更加准确[2],埋设过程中用纱布将压力盒周边绕实,防止二衬混凝土砂块嵌入。二衬混凝土应变计在拱圈内用环向U形支架固定,使其与相应位置二衬外表面切线相平行。总之,在测试仪器埋设过程中,要避免因埋设方法不当使量测“误差”变成“错误”。
2.2 测试方案
为了使测试结果更有研究价值,根据隧道的埋深和上覆土层种类及其厚度的不同,分别在杂填土浅埋段、黄土深埋段、黄土浅埋段各埋设2个量测断面,包括:K3+227、ZK3+160(Ⅴ级杂填土浅埋段,简称227、160断面),ZK3+915、ZK3+243(Ⅴ级深埋段,简称915、243断面),ZK4+770、K4+450(Ⅴ级浅埋段,简称770、450断面),各个断面土层厚度及埋深见表 2。其中450断面埋设于初支与围岩之间的测试元件由于塌方致使损坏,所以浅埋段围岩与初支之间压力规律都由770断面得出。断面布设见图 1,数字为测试元件编号。
表 2 测试断面上覆土层厚度Table 2. Thicknesses of covering soils at test sections3. 初支压力分析
3.1 空间分布规律
初支压力分布见图 2,数字1~13为测点编号,其余数字为初支压力值,单位为kPa。在227、915、770断面中,10、11号压力盒数值都较大,边墙底部承受了大部分的垂直压力。衬砌边墙位置处的压力值也比较大,这是由于拱顶上方土体向两侧挤压所致。227、160、770断面仰拱所受的压力较小,由于这3个断面位于浅埋段,且上覆杂填土或者湿陷性黄土较非湿陷性黄土厚度较大,土体变形速度较快,完成变形所需时间也较短,并且由于施工的滞后导致仰拱没能及时闭合,所以仰拱处的压力值较小。915、243断面位于深埋段,仰拱处的压力值相比浅埋段较大,但增加的幅度不大,由于隧道上覆土层为一定厚度孔隙率相对较小的非湿陷性黄土,易形成自稳性较强的拱圈,比浅埋隧道的稳定性较好,土体承担的释放荷载较大[14]。
整体来看,围岩压力普遍较大,洞周的不均匀性受力比较明显,主要由地层初始应力状态不均匀以及地层的复杂性造成的。最大围岩压力值出现在拱顶位置处,说明拱顶是隧道拱圈最危险的部位。大有山隧道监测数据显示拱顶沉降普遍较大,应力得到了有效释放,由于土体强度较低,自稳能力较差,围岩压力仍然较大[15]。另外,拱肩处的围岩压力也较大,容易发生局部坍塌,在实际埋设压力盒的过程中也验证了这一点,所以在施工过程中应注意该位置的支护工作。
图 3为初支压力沿纵向分布,数字1~13为测点编号。除左拱肩外,915断面的初支压力普遍大于243断面,这是由湿陷性黄土的孔隙率较大引起的[16],而243断面左拱肩位置压力较大则可能由偏压导致。综合浅埋段曲线变化可以看出,由于上覆土体类型不同,随着隧道埋深的增加,部分位置的初支压力随之减小,说明上覆土层的类型对初支压力大小有较大的影响。对比浅埋段和深埋段的初支压力分布,由于浅埋段的初期支护设计更强,有效抑制了土层的变形,从而使得接触压力较大。
表 3列出了利用不同方法计算得到的围岩压力值,其中实测值推算结果是按静力等效近似计算所得,表中θ为滑面的摩擦角,φc为围岩计算摩擦角。对于杂填土断面,无论是按照谢家烋公式[17]、Terza g hi公式[18]还是《公路隧道设计规范》(JTG D70—2004)相应公式来计算,都与实测结果相差较大,其中160断面的实测结果与Terzaghi公式计算结果相对较接近,考虑到杂填土成分复杂,物理力学参数选取不准确,所以227断面计算值均小于实测值。深埋段243断面的实测结果为411.78 kPa,与普氏理论计算结果390.63 kPa比较吻合,而915断面的实测结果大于计算结果且相差较大,说明现有计算方法在计算黄土深埋断面时仍存在一些问题。浅埋770断面的实测结果为508.59 kPa,与Terza g hi公式和《公路隧道设计规范》(JTG D70—2004)计算结果相差较大,而与谢家烋公式所得结果532.25 kPa比较接近。公路隧道设计规范深埋隧道围岩压力计算是基于247座岩质隧道塌方统计得到的经验公式,所以与土质隧道差异大,《公路隧道设计规范》(JTG D70—2004)计算结果与实测结果相比偏小40%~50%。
表 3 竖向围岩压力Table 3. Pressures of vertical surround ling rocks3.2 时间分布规律
图 4为围岩与初支接触压力随时间的变化情况,数字1~9为测点编号。160断面的拱顶、右拱肩和左拱脚处的压力较大,并且时间效应显著,随监测时间增长不断增大,其他测点的压力则较小,只在前30 d略有增长,其后便趋于稳定。结合227、160断面的时程曲线来看,杂填土断面拱顶位置处的压力随时间变化幅度大,并且达到稳定所需的时间较长。深埋段915、243断面压力均比较小,且稳定时间较快,特别是243断面,除拱顶和左拱肩测点外,其余测点压力均不到200 kPa,并且在20 d之内便趋于稳定;拱顶和左拱肩处压力较大,拱顶处最大为1049 kPa,左拱肩处最大值为791 kPa,在监测50 d之后逐渐趋于稳定。770断面位于浅埋段,各个测点的接触压力最大值出现在30 d左右,之后随监测时间增长变化幅度不大,其中,拱腰处的压力随时间的增长幅度比较小,而边墙的压力较大。
图 5为仰拱压力随时间的变化情况,数字10~13为测点编号。仰拱压力盒读数规律与各个断面的拱墙压力规律基本相同,数值普遍较小,说明隧道底部土体产生底鼓的趋势较小,相比之下,杂填土断面仰拱压力值较大且达到稳定所需时间较长,符合杂填土压缩性高、承载力低的特点。
总的来说,在监测2个月以后,围岩压力基本趋于稳定,其中深埋断面稳定所需的时间最短,浅埋段面次之,杂填土断面稳定所需时间最长,这为二衬施作时机的选择提供了依据。通过时程曲线与施工进度对比可发现,各个断面在仰拱施作时,围岩压力呈现减小的趋势,说明仰拱对于提高初支承担围岩压力发挥了一定的作用。
4. 二衬压力分析
4.1 空间分布规律
隧道初支与二衬的接触压力分布见图 6,数字1~9为测点编号,其余数字为接触压力值,单位为kPa。由图 6(a)、(b)可知,227、160断面接触压力的最大值出现在拱顶或拱肩处,并且呈现出由拱顶至拱墙下部逐渐减小的规律。其中,160断面左侧压力大于右侧,这是由隧道洞身两侧地质条件差异以及浇筑混凝土时模板台车自身结构的原因造成的。表 4为杂填土断面初支与二衬分担的荷载比例,P1为初支与围岩接触压力,P2为初支与二衬接触压力,P为总围岩压力。227断面二衬分担比例为17.88%,160断面二衬分担比例为36.17%,两者相差较大,一方面是由于初支还没充分变形情况下便开始施作二衬,另一方面也与杂填土成分复杂、规律性差的特点有关。
表 4 杂填土段断面初支与二衬接触压力Table 4. Contact pressures of primary liner and secondary liner at miscellaneous fill sections915、243断面位于深埋段,由图 6(c)、(d)可知,二衬受力较小,特别是243断面,压力对称分布,说明围岩的自稳能力强,浇筑二衬时,初支压力已经基本稳定,二衬所受压力较小。915断面的压力最大值出现在左拱肩处,最大值为289.30 kPa,该位置处的混凝土浇筑密实度较低,左拱肩处出现应力集中现象。由表 5可知,深埋段二衬的荷载分担比较为接近,均在18%左右,说明在深埋断面中,大部分围岩压力都由初支承担。
表 5 深埋段断面初支与二衬接触压力Table 5. Contact pressures of primary liner and secondary liner at deep depth sections由图 6(e)、(f)可知,770断面压力由拱肩至拱脚呈现逐渐增大的规律。450断面最大压力出现在拱顶位置,最大值为200.64 kPa,拱顶初支的变形较大,在实际施工过程中,当初支施作完成后该断面发生了塌方事故,证实了测试结果的准确性。表 6为770断面初支与二衬接触压力,该断面二衬压力分担比为11.15%,这与新庄岭隧道测试结果10.50%比较接近[19],虽然新庄岭隧道为传统的两层模筑混凝土设计,其计算方法与“新奥法”二衬设计有较大区别,但由于大有山隧道所采取的初期支护比较强,与基于荷载结构原理的两层衬砌设计效果类似,此时二衬主要作为安全储备。
表 6 浅埋段断面初支与二衬接触压力Table 6. Pressures of primary liner and secondary liner at shallow depth sections图 7为初支与二衬接触压力的纵向分布,数字1~9为测点编号。初支与二衬接触压力最大值一般出现在拱顶和拱肩位置,并且分布不均匀,主要由土体分布不均匀以及二衬混凝土浇筑密实度的不同造成的。由于初期支护承受了较大压力,使得浅埋段和深埋段初支与二衬之间接触压力对比并不明显。相对而言,杂填土段160断面的压力分布较为凌乱,这是由杂填土的特点所决定的。由表 4~6可知,大有山隧道二衬荷载承担比在11%~36%范围内,将此结果与以往黄土隧道研究[19-20]进行对比分析,见表 7。大有山隧道的二衬荷载承担比11%~17% 远远小于西安地铁2号线隧道标准断面的测试结果60%[20],而与赵占厂在谗口至兰州高速公路3座黄土隧道(新庄岭隧道、白虎山隧道、青土岘隧道)的测试结果比较接近[19],由此可见,其主要作用仍为安全储备。基于目前黄土公路隧道“新奥法”二衬设计理念,黄土隧道的二衬荷载分担比例比较大,对二衬施工质量的要求也较高;通过大有山黄土隧道的长期监测数据可知,“强初支,弱二衬”的设计可充分发挥初支对围岩变形的控制作用,既达到了传统黄土公路隧道两层衬砌设计的支护效果,又可降低二衬的施作力度,经济效益可观。
表 7 黄土隧道二衬承受的荷载比Table 7. Load ratios of secondlary liners of loess tunnels4.2 时间分布规律
图 8为初支与二衬接触压力随时间的变化情况,数字1~9为测点编号。227、160断面为杂填土,初支与二衬间接触压力整体比较大,特别是上拱圈压力比较大。其中227断面的压力最大值580.12 kPa出现在右拱肩处,160断面的压力最大值483.30 kPa出现在拱顶处。在拱墙上部测点,拱顶和左右拱肩3个测点的变化最明显,接触压力随着监测时间持续增大。竖向围岩压力主要由隧道拱圈承担,在岩质隧道中,拱圈顶部形成一个不受围岩约束的脱离区,而在黄土隧道中,特别是杂填土,土层与拱圈顶部密贴使拱圈顶部形成一个抗力区,但由于杂填土较松散,土体变形受到拱圈约束,使得接触压力不断增大,变形稳定时间也较长。915、243断面位于深埋段,由图 8(b)可知,二衬受力较小,承载作用不明显。初支与二衬间接触压力40 d之后便基本趋于稳定。770、450断面位于浅埋段,拱顶和拱脚处二衬受力较大,拱腰较小,但从整体来看,各个测点受力较小,稳定较快。其中450断面的拱顶和右拱脚处的压力比较大,并且随着监测时间不断增大,考虑到偏压的影响,作用的围岩应力不断释放的缘故,建议隧道支护结构采取不对称的设计参数。
随着监测时间的增长,这几个断面在40 d内均呈现先增大后减小再增大的变化规律,而压力开始减小的时刻正是仰拱开始施作时。在仰拱开挖或撤除临时支护后,隧道边墙底部由三向应力状态变为单向受力状态,衬砌结构产生向洞内方向的变形,使得接触压力逐渐减小;然后随着仰拱混凝土强度的增大,应力状态发生改变,接触压力又缓慢增大。另外,二衬压力由开始监测到最后趋于稳定所需要的时间基本在二衬施作后的40 d左右。
5. 二衬应变分析
5.1 空间分布规律
二衬与仰拱混凝土应变分布见图 9,数字1~13为测点编号,其余数字为混凝土应变,单位为10-6。227断面最大应变出现在左右拱脚处,并且左侧应变比右侧应变大,这与压力盒所测数据基本一致。整体来看,这2个断面不同位置处混凝土的应变差别较大,其中227断面的拱顶处和160断面的仰拱右侧为压应变,其余部分为拉应变。深埋段915、243断面二衬应变都比较小,并且不同位置应变的差别相对较小,说明该断面处结构变形稳定较快,其中915断面的最大应变为181.13×10-6,位于仰拱中心偏左处,作用在仰拱底部的土层,承担了较大的荷载。浅埋段断面(770、450)二衬应变最大值出现在拱脚和拱顶附近,尤其是450断面的拱脚处,混凝土应变最大值为436.25×10-6。在二衬施作后,拱脚和拱顶处初支的残余变形较大,导致这2个位置处的二衬应变值较大,这也说明对于不同的围岩,其侧压力系数存在着差异性,不能一概而论。
5.2 时间分布规律
图 10为二衬混凝土应变随时间的变化情况,数字1~9为测点编号。杂填土段227与160断面在初测阶段除左拱腰外,其他测点处均为压应变,应变先增大后逐渐减小并演变成拉应变,而拱顶处则一直为压应变,并且随时间不断增大;深埋断面各个测点处的应变普遍都比较小,除243断面的拱脚测点,多数都在(-40~40)×10-6范围内,说明了隧道结构变形稳定的时间较快;浅埋段770与450断面各个测点应变值随时间增长规律比较接近,开始不断增大,随后在监测60 d左右基本稳定。图 11为仰拱混凝土应变随时间的变化情况,数字10~13为测点编号。仰拱测点以拉应变为主,应变值在监测10~ 20 d急剧增大。随后逐渐趋于稳定,但稳定后仍有微小的波动,这是由于施工中行车等因素造成的。
整体来看,二衬混凝土在初测阶段均为压应变并不断增大,到达峰值之后逐渐减小,之后变成拉应变,并逐渐增大,最后趋于稳定,这个规律在227与243断面更为明显。
6. 结语
(1)在浅埋杂填土段,拱顶和拱肩位置处的围岩压力较大且波动大,基本呈现出由拱顶至边墙逐渐减小的规律;在深埋段,围岩压力最大值出现在拱顶和左拱肩处,分别为1049、791 kPa;在浅埋段,边墙围岩压力较大,约为1050 kPa,而仰拱受力较小,为使仰拱得到充分利用,在施工阶段应及时闭合成环。
(2)围岩压力在2个月后基本趋于稳定,其中深埋断面稳定所需的时间最短,浅埋断面次之,杂填土断面时间最长。
(3)拱顶、拱肩处的压力值较大,是黄土隧道最危险的部位,容易发生坍塌,设计和施工过程中应注意超前支护参数的选取。
(4)围岩压力计算值与实测值差距较大,仅深埋段断面的实测结果与普氏理论的计算结果相吻合,说明现有方法在计算黄土断面围岩压力时存在较大的误差,本文实测值可为相似地域工程提供一定参考。
(5)软弱黄土隧道设计施工时,在“新奥法”原理指导下,可进一步遵循“强初支,弱二衬”的设计理念,充分发挥初支对围岩变形的抑制作用,一定程度上既达到了传统黄土隧道两层衬砌设计的支护效果,又可降低二衬的施作力度,使二衬作为安全储备。
-
表 1 转向架周围不同风速区间的空间体积
Table 1. Spatial volumes in different wind speed ranges around bogies
m3 转向架 0~5 m·s-1 5~10 m·s-1 10~15 m·s-1 15~20 m·s-1 b1 1.056 1.999 1.730 1.068 b2 1.761 3.041 1.610 0.854 b3 1.545 2.563 1.609 1.033 b4 2.584 2.540 1.140 0.907 b5 2.051 2.375 1.704 1.301 b6 2.216 3.018 1.345 0.711 -
[1] KLOOW L. High-speed train operation in winter climate[R]. Stockholm: KTH Railway Group and Transrail, 2011. [2] SHISHIDO M, NAKADE K, IDO A, et al. Development of deflector to decrease snow-accretion to bogie of vehicle[J]. Rtri Report, 2009, 23 (3): 29-34. [3] 韩运动, 姚松, 陈大伟, 等. 高速列车转向架舱内流场实车测试和数值模拟[J]. 交通运输工程学报, 2015, 15 (6): 51-60. doi: 10.3969/j.issn.1671-1637.2015.06.007HAN Yun-dong, YAO Song, CHEN Da-wei, et al. Real vehicle test and numerical simulation of flow field in high-speed train bogie cabin[J]. Journal of Traffic and Transportation Engineering, 2015, 15 (6): 51-60. (in Chinese). doi: 10.3969/j.issn.1671-1637.2015.06.007 [4] 丁叁叁, 田爱琴, 董天韵, 等. 端面下斜导流板对高速列车转向架防积雪性能的影响[J]. 中南大学学报(自然科学版), 2016, 47 (4): 1400-1405. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201604041.htmDING San-san, TIAN Ai-qin, DONG Tian-yun, et al. Influence of inclined guiding plate on anti-snow performance of high-speed train bogie[J]. Journal of Central South University (Science and Technology), 2016, 47 (4): 1400-1405. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201604041.htm [5] 何德华, 王刚义, 陈厚嫦, 等. 基于风雪两相流的高寒动车组转向架防冰雪绕流技术研究[J]. 铁道机车车辆, 2016, 36 (4): 38-42. doi: 10.3969/j.issn.1008-7842.2016.04.09HE De-hua, WANG Gang-yi, CHEN Hou-chang, et al. Analysis on two-phase flowing anti-ice/snow of low-temperature EMU bogie[J]. Railway Locomotive and Car, 2016, 36 (4): 38-42. (in Chinese). doi: 10.3969/j.issn.1008-7842.2016.04.09 [6] 冯永华, 黄照伟, 张琰, 等. 高寒动车组转向架区域积雪结冰数值仿真研究[J]. 铁道科学与工程学报, 2017, 14 (3): 437-444. doi: 10.3969/j.issn.1672-7029.2017.03.002FENG Yong-hua, HUANG Zhao-wei, ZHANG Yan, et al. Research of numerical simulation of the snow icy phenomenon of the high-speed train bogie area[J]. Journal of Railway Science and Engineering, 2017, 14 (3): 437-444. (in Chinese). doi: 10.3969/j.issn.1672-7029.2017.03.002 [7] RAGHUNATHAN R S, KIM H D, SETOGUCHI T. Aerodynamics of high-speed railway train[J]. Progress in Aerospace Sciences, 2002, 38 (6/7): 469-514. [8] SCHULTE-WERNING B. Research of European railway operators to reduce the environmental impact of high-speed trains[J]. Journal of Rail and Rapid Transit, 2003, 217 (4): 249-257. doi: 10.1243/095440903322712856 [9] PAZ C, SUÁREZ E, GIL C, et al. Effect of realistic ballasted track in the underbody flow of a high-speed train via CFD simulations[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2019, 184: 1-9. doi: 10.1016/j.jweia.2018.11.007 [10] PAZ C, SUÁREZ E, GIL C. Numerical methodology for evaluating the effect of sleepers in the underbody flow of a high-speed train[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2017, 167: 140-147. doi: 10.1016/j.jweia.2017.04.017 [11] 陈羽, 杨志刚, 高喆, 等. 底部结构对高速列车流场及气动优化规律的影响[J]. 同济大学学报(自然科学版), 2016, 44 (6): 930-936. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201606016.htmCHEN Yu, YANG Zhi-gang, GAO Zhe, et al. Influences of underbody structures on flow field and aerodynamic optimization laws of high speed train[J]. Journal of Tongji University (Natural Science), 2016, 44 (6): 930-936. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201606016.htm [12] XIA Cao, WANG Han-feng, SHAN Xi-zhuang, et al. Effects of ground configurations on the slipstream and near wake of a high-speed train[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2017, 168: 177-189. doi: 10.1016/j.jweia.2017.06.005 [13] 郗艳红, 毛军, 高亮, 等. 横风作用下高速列车转向架非定常空气动力特性[J]. 中南大学学报(自然科学版), 2014, 45 (5): 1705-1714. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201405045.htmXI Yan-hong, MAO Jun, GAO Liang, et al. Aerodynamic force/moment for high-speed train bogie in crosswind field[J]. Journal of Central South University (Science and Technology), 2014, 45 (5): 1705-1714. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201405045.htm [14] ZHANG Jie, LI Jing-juan, TIAN Hong-qi, et al. Impact of ground and wheel boundary conditions on numerical simulation of the high-speed train aerodynamic performance[J]. Journal of Fluids and Structures, 2016, 61: 249-261. doi: 10.1016/j.jfluidstructs.2015.10.006 [15] SHUR M L, SPALART P R, STRELETS M K, et al. A hybrid RANS-LES approach with delayed-DES and wall-modelled LES capabilities[J]. International Journal of Heat and Fluid Flow, 2008, 29 (6): 1638-1649. doi: 10.1016/j.ijheatfluidflow.2008.07.001 [16] GRITSKEVICH M S, GARBARUK A V, SCHUTZE J, et al. Development of DDES and IDDES formulations for the k-ω shear stress transport model[J]. Flow, Turbulence and Combustion, 2012, 88 (3): 431-449. doi: 10.1007/s10494-011-9378-4 [17] FAVRE T, EFRAIMSSON G. An assessment of detached-eddy simulations of unsteady crosswind aerodynamics of road vehicles[J]. Flow, Turbulence and Combustion, 2011, 87 (1): 133-163. doi: 10.1007/s10494-011-9333-4 [18] MORDEN J A, HEMIDA H, BAKER C J. Comparison of RANS and detached-eddy simulation results to wind-tunnel data for the surface pressures upon a class 43 high-speed train[J]. Journal of Fluids Engineering, 2015, 137 (4): 1-9. [19] 苗秀娟, 高广军. 基于DES的车辆横风性能模拟[J]. 中南大学学报(自然科学版), 2012, 43 (7): 2855-2860. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201207058.htmMIAO Xiu-juan, GAO Guang-jun. Aerodynamic performance of train under cross-wind based on DES[J]. Journal of Central South University (Science and Technology), 2012, 43 (7): 2855-2860. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201207058.htm [20] 张亮, 张继业, 李田, 等. 高速列车不同位置受电弓非定常气动特性研究[J]. 机械工程学报, 2017, 53 (12): 147-155. https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201712018.htmZHANG Liang, ZHANG Ji-ye, LI Tian, et al. Research on unsteady aerodynamic characteristics of pantographs in different positions of high-speed trains[J]. Journal of Mechanical Engineering, 2017, 53 (12): 147-155. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201712018.htm [21] 潘永琛, 姚建伟, 梁策, 等. 高速列车近尾流区涡旋结构的湍流特性分析[J]. 中国铁道科学, 2017, 38 (2): 83-88. doi: 10.3969/j.issn.1001-4632.2017.02.13PAN Yong-chen, YAO Jian-wei, LIANG Ce, et al. Analysis on turbulence characteristics of vortex structure in near wake of high speed train[J]. China Railway Science, 2017, 38 (2): 83-88. (in Chinese). doi: 10.3969/j.issn.1001-4632.2017.02.13 [22] YAO Shuan-bao, SUN Zhen-xu, GUO Di-long, et al. Numerical study on wake characteristics of high-speed trains[J]. Acta Mechanica Sinica, 2013, 29 (6): 811-822. doi: 10.1007/s10409-013-0077-3 [23] MULD T W, EFRAIMSSON G, HENNINGSON D S. Wake characteristics of high-speed trains with different lengths[J]. Journal of Rail and Rapid Transit, 2014, 228 (4): 333-342. doi: 10.1177/0954409712473922 [24] 余以正, 姜旭东, 孙健. 不同排障器导流罩对高速列车阻力及升力的影响[J]. 大连交通大学学报, 2017, 38 (6): 30-34. https://www.cnki.com.cn/Article/CJFDTOTAL-DLTD201706008.htmYU Yi-zheng, JIANG Xu-dong, SUN Jian. Drag force and lift force performance of a high speed train affected by different cowcatcher air-deflector[J]. Journal of Dalian Jiaotong University, 2017, 38 (6): 30-34. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-DLTD201706008.htm [25] SUJUDI D, HAIMES R. Identification of swirling flow in 3-D vector fields[C]//AIAA. 12th Computational Fluid Dynamics Conference. Reston: AIAA, 1995: 792-799. [26] 李相鹏, 周昊, 岑可法. 涡结构对小颗粒在圆管背风面碰撞分布的影响[J]. 浙江大学学报(工学版), 2006, 40 (4): 605-609. doi: 10.3785/j.issn.1008-973X.2006.04.013LI Xiang-peng, ZHOU Hao, CEN Ke-fa. Influences of vortices on impaction and distribution of small ash particles on rear side surface of boiler tube[J]. Journal of Zhejiang University (Engineering Science), 2006, 40 (4): 605-609. (in Chinese). doi: 10.3785/j.issn.1008-973X.2006.04.013 [27] LYU Xiao-hui, HUANG Ning, TONG Ding. Wind tunnel experiments on natural snow drift[J]. Science China Technological Sciences, 2012, 55 (4): 927-938. doi: 10.1007/s11431-011-4731-3 [28] 周晅毅, 顾明, 朱忠义, 等. 首都国际机场3号航站楼屋面雪载荷分布研究[J]. 同济大学学报(自然科学版), 2007, 35 (9): 1193-1196. doi: 10.3321/j.issn:0253-374X.2007.09.008ZHOU Xuan-yi, GU Ming, ZHU Zhong-yi, et al. Study on snow loads on roof of terminal 3 of Beijing Capital International Airport[J]. Journal of Tongji University (Natural Science), 2007, 35 (9): 1193-1196. (in Chinese). doi: 10.3321/j.issn:0253-374X.2007.09.008 [29] KIND R J, MURRAY S B. Saltation flow measurements relating to modeling of snowdrifting[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1982, 10: 89-102. doi: 10.1016/0167-6105(82)90056-3 [30] SATO T, KOSUGI K, MOCHIZUKI S, et al. Wind speed dependences of fracture and accumulation of snowflakes on snow surface[J]. Cold Regions Science and Technology, 2008, 51: 229-239. doi: 10.1016/j.coldregions.2007.05.004 -