低真空管道交通系统列车气动阻力研究.pdf
《低真空管道交通系统列车气动阻力研究.pdf》由会员分享,可在线阅读,更多相关《低真空管道交通系统列车气动阻力研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、技 术 装 备2023/07CHINA RAILWAY低真空管道交通系统列车气动阻力研究王友彪,宣言,张俊博(中国铁道科学研究院集团有限公司 铁道科学技术研究发展中心,北京 100081)摘要:采用计算流体动力学软件,结合动网格技术,建立全长200 m的8节全尺寸列车模型,基于三维瞬态可压缩黏性流体数值模拟,分析低真空管道交通系统列车的气动阻力分布规律,以及真空度、阻塞比、管道截面形状、列车运行速度等因素与列车气动阻力的关系。研究表明,列车气动阻力主要由头车和尾车产生,中间各节车气动阻力不同但均占比不大,气动阻力成分主要为压差阻力,且远大于摩擦阻力;列车气动阻力随真空度、阻塞比增大而增大,且均
2、近似成线性关系;在同阻塞比情况下,方形截面管道内列车气动阻力大于半圆形截面管道内列车气动阻力;在低真空管道环境下,列车气动阻力与速度近似为二次关系。关键词:低真空管道;气动阻力;真空度;阻塞比;仿真分析中图分类号:U260.17 文献标识码:A 文章编号:1001-683X(2023)07-0053-07DOI:10.19549/j.issn.1001-683x.2023.02.13.0010 引言 列车在开放空间的大气环境下运行,头部长度为10 m的流线型列车在运行速度300 km/h时,气动阻力占总阻力的75%1;当运行速度进一步提高时,气动阻力占比更高。气动阻力成为列车运行速度进一步提升
3、的关键制约因素之一。在低真空环境下,列车气动阻力大幅减小,可实现列车更高速运行2。因此,低真空管道高速列车作为新一代运输方式被学者提出3。近年来,随着Hyperloop真空管道概念的提出4,以及我国西南交通大学、中国中车股份有限公司、中国航天科工集团第三研究院等单位的大力推进,低真空管道交通系统发展迅猛5。空气动力学是低真空管道交通系统的一个重要研究领域6,国内外对一系列相关问题开展了研究,仿真研究是其中的重要研究手段之一7。周晓等8基于二维可压缩黏性流体,研究了管道气压 103 Pa、基金项目:中国国家铁路集团有限公司科技研究开发计划项目(N2021J033);中国铁道科学研究院集团有限公司
4、科研开发基金项目(2021YJ026)第一作者:王友彪(1988),男,助理研究员。E-mail:-53低真空管道交通系统列车气动阻力研究 王友彪 等技 术 装 备2023/07CHINA RAILWAY阻塞比0.160.78、列车运行速度180720 km/h时列车气动阻力的变化规律;Kim 等9采用静止列车模型系统地研究了阻塞比、真空度与列车运行速度的关系,并 计 算 出 列 车 极 限 速 度 和 激 波 临 界 速 度;刘加利等10-12基于三维可压缩计算流体力学数值模拟,建立“头车+中间车+尾车”的简化 3 节车模型,分析了列车气动阻力与管道真空度、阻塞比、列车运行速度的关系;Llu
5、esma-Rodrguez 等13研究指出,真空管道内气动阻力系数随阻塞比增大而增大,但几乎不随真空度变化;Le等14进一步分析了不同真空度、阻塞比等参数条件下,摩擦阻力和压差阻力的变化规律。目前,针对低真空管道环境下列车空气动力学三维全尺寸模型的数值研究较少。管道直径越大、真空度越高,建设成本越大;管道直径越小、真空度越低,列车运行阻力越大。真空度、阻塞比等指标成为平衡建设成本与运行成本的关键指标。因此,有必要对低真空管道高速列车空气动力学进行系统性研究,探索三者关系,为低真空管道轨道交通的发展提供理论支撑。以低真空管道交通系统为研究对象,利用计算流体力学对管道内超高速列车的空气动力学特性进
6、行三维全尺寸数值模拟研究,系统地分析真空度、阻塞比、管道截面形状、列车运行速度等对气动阻力的影响规律。1 仿真计算方法 1.1基本控制方程目前,高速列车商业运行速度普遍在 350 km/h(马赫数 0.286)以下,可将空气简化为不可压缩流体,即假定密度为常数;而对于低真空管道交通系统,预计其运行速度可达 1 000 km/h,考虑该研究的速度范围为 4001 000 km/h,马赫数均0.300,仿真模型应考虑空气的压缩效应。因此,采用三维瞬态可压缩 Navier-Stokes 方程为控制方程,表示如下:t+xi(ui)=0 ,(1)t(ui)+xi(uiuj)=-pxi+xj (uixj+
7、ujxi-23ijuixi),(2)t(E)+xiui(E+p)=xj kTxj+ui(uixj+ujxi)-23uiuixiij ,(3)E=T+12(u2i+u2j+u2k),(4)p=RT ,(5)式中:为流体密度;t为时间;u为流体速度;p为流体压强;E为内能;为流体动力黏度;T为流体温度;k为导热率;为流体比热容;R为气体常数。1.2湍流模型与列车相关的流动现象绝大多数是湍流流动。湍流的数值模拟方法主要有直接数值模拟湍流、大涡模拟和湍流模型(包括涡黏性模型、雷诺应力模型)。涡黏性模型和雷诺应力模型中的-两方程模型在列车周围空气流场流动的湍流数值模拟中应用最广泛,因此采用-湍流模型进行
8、湍流模拟。2 模型建立 2.1列车模型计算车型选用标准动车组作为低真空管道交通系统的运行列车。根据某动车组实际尺寸,建立动车组三维几何模型(见图1)。列车全长200 m,包括1节头车、6节车身和1节尾车。列车保留对气动阻力起主要作用的车身,可在保证达到研究目标的前提下提高计算效率。2.2低真空管道模型仿真模拟基于“动模型”方法,列车在低真空管道内运行。为保证流动充分发展,经过试算低真空管道长度选为2 453 m,作为流体计算外域,列车在管道内运行一段时间后,流动及气动阻力趋于稳定。低真空管道模型见图2。图1动车组三维几何模型-54低真空管道交通系统列车气动阻力研究 王友彪 等技 术 装 备20
9、23/07CHINA RAILWAY如图2所示,低真空管道截面为半圆形,阻塞比为0.180.70,动车组放置于管道中间位置。不同阻塞比条件下的低真空管道计算外域模型形状见图3。同时,为考察管道截面形状对列车气动阻力的影响,建立方形截面低真空管道作为对比。2.3计算网络列车高速运行属于非定常问题,为模拟列车与隧道之间的相对运动,计算区域采用分区对接网格技术,地面、列车均采用六面体结构化网格离散。各分区之间的数据交换通过公共滑移界面进行。研究中,列车运行速度高达1 000 km/h,对流体动力学仿真计算方法的准确性要求更高。因此,需要对比不同网格数量、边界条件、湍流方程对计算结果的影响规律,确定更
10、准确的流体动力学计算方法。对比分析500万、750万、1 000万等不同粗细网格条件下的仿真计算结果。结果表明,网格数量达到一定数量后,结果变化不大。考虑细网格所需计算资源可以接受,因此,后续研究中采用网格数量均约为1 000万。2.4初始及边界条件设置不同初始条件及边界条件,以实现不同仿真工况:(1)设置不同的管道环境压力,模拟不同真空度:分别为500、1 000、2 000、3 000、4 000 Pa;(2)设置不同的滑移区域运动速度,模拟列车运行分别为:400、500、600、700、800、900、1 000 km/h。(3)设置出入口为压力出口边界,列车、轨道及隧道等结构表面设置为
11、固定壁面边界。3 结果分析 3.1气动阻力分布规律各工况条件下列车气动阻力规律基本相同,选择管道压力 500 Pa,列车运行速度 1 000 km/h,阻塞比0.18(车体横截面积 12.8 m2、半圆形隧道横截面积72.0 m2)为典型工况。将长度200 m列车分为8节,对各节车厢气动阻力进行单独监控(见表1)。由表1可知,每节车所受气动阻力不同,尾车气动阻力最大,头车其次,其他6节中间车气动阻力差别不大,均40 N,远低于头车和尾车的气动力。上述气动阻力分布规律与列车在常压、开放环境下运行基本一致。气动阻力包括压差阻力和摩擦阻力,头车、尾车产生的气动阻力主要为压差阻力,而中间各节车主要为摩
12、擦阻力。另外,中间各节车形状完全一致,但边界层沿列车方向发展,因此各节车的气动阻力也有所区别。考虑列车气动优化重点关注头车,后续分析中气动阻力均为列车头车的气动阻力。表1各节车厢气动阻力车厢编号头车2车3车4车5车6车7车尾车阻力/N3 301.00.217.120.022.926.333.75 471.0(a)阻塞比0.18(d)阻塞比0.40(e)阻塞比0.50(b)阻塞比0.20(f)阻塞比0.60(c)阻塞比0.30(g)阻塞比0.70图3低真空管道计算外域模型形状图2低真空管道模型-55低真空管道交通系统列车气动阻力研究 王友彪 等技 术 装 备2023/07CHINA RAILWA
13、Y3.2阻塞比对气动阻力的影响设置列车运行速度1 000 km/h,真空度500 Pa,计算不同阻塞比条件下的列车和低真空管道的横截面积(见表2),分析不同阻塞比对列车气动阻力的影响。对列车气动阻力进行监控,不同阻塞比条件下的列车气动阻力时程曲线见图4。由图4可知,低真空管道内列车运行各工况的气动阻力曲线形状几乎一致,在开始阶段不断升高,较短时间(约0.5 s)后气动阻力趋于稳定,说明列车周围的流场已处于稳定状态。列车气动阻力随阻塞比的增大而增大。稳定后,不同阻塞比条件下的列车气动阻力见表3。根据表3数据,绘制不同阻塞比条件下的列车气动阻力拟合曲线(见图5)。选择拟合曲线形式为幂函数形式,即y
14、=axb,式中:a=16 434;b=0.941。由图5可知,列车气动阻力与阻塞比近似成线性关系。3.3真空度对气动阻力的影响设置列车运行速度1 000 km/h,阻塞比0.18(车体横截面积12.8 m2、隧道横截面积72.0 m2),计算不同真空度条件下的列车气动阻力(见表4),分析不同真空度对列车气动阻力的影响规律。管道压力分别设置为500、1 000、2 000、3 000、4 000 Pa。不同真空度条件下,列车气动阻力时程曲线见图6。由图6可知,不同真空度条件下的气动阻力随时间的变化规律几乎一致,在较短时间后趋于稳定,且随着真空度的增大,列车气动阻力不断增大。根据表4数据,绘制不同
15、真空度条件下的列车气动阻力拟合曲线(见图7)。拟合曲线公式为y=axb,式中:a=6.93;b=0.991,可知列车气动阻力与真空度近似成线性关系。图6不同真空度条件下的列车气动阻力时程曲线图5不同阻塞比条件下的列车气动阻力拟合曲线表3不同阻塞比条件下的列车气动阻力N阻塞比0.180.200.300.400.500.600.70最大值3 3013 5815 2267 0148 50010 20011 721图4不同阻塞比条件下的列车气动阻力时程曲线表2不同阻塞比条件下的列车、低真空管道横截面积m2阻塞比0.180.200.300.400.500.600.70横截面积列车12.812.812.8
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 真空 管道 交通 系统 列车 气动 阻力 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。