CFD数值仿真在高速列车设计中的应用Word文档下载推荐.docx
- 文档编号:18980505
- 上传时间:2023-01-02
- 格式:DOCX
- 页数:10
- 大小:716.43KB
CFD数值仿真在高速列车设计中的应用Word文档下载推荐.docx
《CFD数值仿真在高速列车设计中的应用Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《CFD数值仿真在高速列车设计中的应用Word文档下载推荐.docx(10页珍藏版)》请在冰豆网上搜索。
A
ApplicationofCFDNumericalSimulationinHighSpeedTrainDesign
WANGDongping,ZHAOWenzhong,MASiqun
(SchoolofMechanicalEngineering,DalianJiaotongUniversity,Dalian116028,China
Abstract:
Theaerodynamicperformanceofhighspeedtrain(HSTisrelatedtothepropulsiveefficiencyandhasinfluenceontheridingcomfortabilityandsecurity.CFDnumericalsimulationofHSTapplicationsisintroducedinthispaper:
theventsystemofpassengertrainsbyAirpak;
waterswayinginwaterdrumsofpassengertrains,outerflowfieldsofhighspeedtrainsandaerodynamicnoisesgeneratedbytheflowfluidaroundthe2dstreamlineHSTmodelbyFluent.Thecomputationalresultoftheventsystemagreeswiththeexperimentalresult.ThebucklinganalysisfromIDEASintegratingFluentagreeswiththefailurephenomenonofthewaterdrum.
Keywords:
CFDnumericalsimulation;
highspeedtrain;
viscosityresistance;
ventsystem;
waterswayinginawaterdrum;
aerodynamicnoise
近年来,随着计算机计算能力的不断提高,以及计算流体动力学(CFD本身理论和方法的不断改进,CFD已经被越来越多地应用到各个工程领域中。
通过计算机模拟获得某种流体在特定条件下的有关信息,是工程设计人员用于分析问题和解决问题的强有力工具。
高速列车的特点是庞大、细长、在地面轨道上运行,其空气动力学问题非常复杂。
过去,由于列车运行速度较低,空气阻力等空气动力学现象不显著,随着列车的高速化,空气动力学问题越发突出。
高速列车的设计提出了许多需要解决的问题。
诸如如何确定列车两侧人员的安全避让距离、列车尾流人的安全避让距离和如何确定列车的远场空气噪声收稿日期:
20061222;
修回日期:
20070602
基金项目:
国家高技术研究发展计划项目资助(2001AA411110作者简介:
王东屏(1962,女,陕西耀县人,教授,博士。
Email:
wdp1962@163.com等。
在列车启动和制动过程中,列车水箱中水的晃动对水箱产生的冲击力,也是设计者需要掌握的信息。
由于空气阻力与运行速度的平方成正比,所以高速列车运行时,其空气阻力约占列车全部阻力的80%,因而在气动阻力最小化的要求下,高速列车的外型设计不仅要求几何外形具有时代气质,还要求有良好的空气动力学性能。
而空气动力学性能是否优秀,在设计阶段只能借助于CFD数值仿真进行分析与评估。
日本、德国等高速列车发达国家都将高速列车空气动力学问题和列车头型设计作为关键技术进行研究,关于高速列车与空气相互作用的研究,己逐步形成了一个独立的学科高速列车空气动力学。
国内高速列车在空气动力性能数值仿真方面滞后于需求,因此,本文对如下所示高速列车在CFD领域内外流场的一些问
供技术储备。
列车头部
列车尾部
列车远场噪声
列车空调通风
列车水箱晃动
1算法简介
描述流场计算的控制方程是一组非线性偏微分方程组,为了能利用计算机进行求解,需要将这一组微分方程转化成每一个节点上的一组代数方程,该方程组中包含有该节点及邻近节点上所求函数的值,即离散方程。
建立离散方程的方法很多,本文采用的数值方法是有限体积法,即将守恒型的控制方程对区域离散后形成的控制容积作积分,对于节点间物理量的变化特性作出假设,从而把积分进行到底,得出节点间物理量间的代数方程式。
在将控制方程对控制容积作积分的过程中,需要对所求解的变量在两个节点之间的变化特性作出假设,而不同的假设就导致不同的离散格式。
本文中,不同的模型流场计算采用不同的离散格式,如表1所示。
表1各应用实例算法一览表
应用实例湍流模型算法
离散格式
对流项扩散项
空间维数时间因素空调通风系统零方程模型SIMPLE二阶迎风格式二阶精度中心差分三维定常气动阻力方程SIMPLE二阶迎风格式二阶精度中心差分三维定常水箱晃动方程PISO二阶迎风格式二阶精度中心差分三维非定常远场气动噪声大涡模拟PISO二阶迎风格式二阶精度中心差分二维非定常
湍流是一种高度复杂的非稳态三维流动,在湍流中流体的各种物理参数,如速度、压力等都随时间与空间发生随机的变化。
对于湍流的数值计算方法可以分为3类[1]:
直接模拟(DNS、大涡模拟(LES、雷诺时均方程法(RANS。
目前,以雷诺时均方程为基础的湍流模型模拟方法是处理工程实际问题更有效的、更常用的方法。
由于在讨论高速列车内外流场的气动特性时,感兴趣的是时均速度场、压力场,不需要知道湍流产生与发展的细节,因此在计算高速列车外流场、空调通风系统和水箱晃动中均选用雷诺时均方程为基础的湍流模型模拟方法。
但采用雷诺时均方程不能捕捉到宽带湍流噪声的信息,方程不能模拟由于耗散特性引起的压力振颤,所以在求解列车远场气动噪声时,采用大涡模拟方法。
LES是介于DNS和RANS方程之间的一种方法,它能直接求解流体流动中大于网格尺度的涡,因而宽带的气动噪声可以由大涡仿真求解,以便精确模拟流动中最主要的由于噪声产生的非定常特征。
本文讨论的4个流场均为不可压缩流动,空调通风系统和列车外流场计算采用算法SIMPLE,而水箱晃动和列车远场气动噪声采用算法PISO,如表1所示。
性能仿真的关键在于建模,正确建立模型需要有良好的理论基础和丰富的建模经验。
如何建立一个既的物理模型是第一关键。
由于工程设备结构的复杂性,在对发生于其间的流动和传热过程进行数值计算时,常常需要对结构做合理的简化,即在对实际物理流动现象深刻了解的基础上,把实际模型结构提炼成适用于空气动力流场计算的模型。
2商用软件简介
CFD通用软件包的出现并且商业化,对CFD技术在工程应用中的推广起了巨大的促进作用,本文使用美国Fluent公司的产品Fluent和Airpak软件包对列车内外流场进行数值研究。
Fluent是一个用于模拟和分析在复杂几何区域内的流体流动与热交换问题的专用CFD软件,它针对非结构化网格模型设计,是求解不可压缩流及可压缩流流场问题的CFD求解器,采用目前应用广泛又较成熟的有限体积法的数值计算方法,提供形象逼真、可视化的后处理结果。
Airpak是一个专用的空调通风系统设计软件,利用Fluent求解器进行流场计算,可计算超复杂空间模型的流动传热问题。
尽管有通用性很强的商业软件,但不存在可以解决所有问题的软件。
计算结果是否可靠与软件使用者对计算方法的理解以及对具体流动问题的理解很有关65
第5期CFD数值仿真在高速列车设计中的应用
3空调通风系统的数值仿真及验证
我国铁路客车空调通风起步较晚,与先进国家相比,仍然处于落后状态。
客车空调通风系统的设计基本上是按照我国铁路标准中规定的铁路客车热工计算方法及相关的客车技术条件要求做一些简单的计算,因而存在不能准确定量地进行风量分配,导致客车纵向断面和横向断面温度、风量分布不均匀现象,满足不了高水平的舒适性设计要求。
当前空调通风的大致估算、经验设计,成为制约准高速/高速客车整体设计水平的重要因素之一。
本文采用软件Airpak数值仿真客车空调通风系统,并将数值模拟结果与实验结果进行对比以考察CFD模型的可靠性。
3.1建模分析
由于空调领域内的流动问题多为低速流动,流体密度变化不大,故可将其看作不可压缩流动。
工程中常采用的湍流模型是两方程模型或其变形,但是模型对于等温流动情形能模拟得很好,对于空调通风系统内的非等温、混合对流流动的模拟却有较大的误差。
近年来,一些学者提出用零方程湍流模型对所关心的问题进行模拟。
由于湍流模型模拟是半经验的,故尽管零方程模型比较简单,但在专门的领域内却能获得比复杂模型更符合实际的结果[2],所以本文采用零方程湍流模型。
边界条件采用壁面函数法处理。
3.2计算结果分析
图1是摆式客车一等座车风道模型三维图,空调
系统中的气流从左往右流动。
数值仿真的优越性之一是方案对比,在计算过程中发现,当在风道静压腔前半程中设置隔板时,车厢前后温度均匀性较好,温度变化范围大约是18.5~21.7,最大温度差为3.2,基本符合设计要求,而在静压腔中全程设置隔板时,车厢前后温度分布均匀性差一些,同样位置温度变化范围大约是18~23,列车前半部分温度相对偏低,后半部分温度相对偏高。
原因是主风道的速度逐渐在减小,如果在主风道的后半部分再加隔板,主风道的气流再经过静压腔
中衰减,而从风口出风的风速将更加减小,导致车厢后半部分的温度相对较高。
通过CFD数值模拟,可以对客车空调通风系统模型进行多方面分析,反复修改模型,进而达到满意的结果。
经过多方案比较,摆式客车一等座车的空调通风系统设计得到明显的改善:
风道的设计更加合理,风口的出风量比较均匀,车厢内的温度分布也比较均匀,为工厂的设计提供量化指导。
图2是动车组头车空调通风系统几何模型。
由于风道结构复杂,计算区域的边界完全按照车厢内部形状建模,风道中的隔板与人体模型之间的间隙小,需要很密的网格才能满足计算精度,因此,用一个尺寸相当的块体代替人体模型,明显减少了网格数,同时经过试算也证明,
这样简化对结果的影响不大。
图3是空调机组端一位测风道出风口速度试验与数值模拟比较。
从图中可看到:
一位侧风口风速的试验与计算值具有良好的一致性,展示出相同的变化规律。
试验结果与数值计算结果的良好一致性证实本文采用的CFD
模型是可靠的。
4流线型高速列车气动阻力研究
列车外形关系着列车的空气阻力、列车风、列车交会压力波、列车在隧道内的微气压波、气动噪声等空气动力学问题。
列车头部流线型设计可以明显改善列车的空气动力学特性,有效地减小列车的运行阻力,使得相同速度下的列车空气阻力比传统列车减少60%以上,极大地推动了列车运行速度的提高,同时还可以降低隧道中出现的压力瞬变值,减弱列车尾流,降低列车的气动噪声[3]
。
本文应用Fluent软件对流线型高速车体外流场进行三维模拟计算。
66
铁道学报第29卷
4.1建模分析
计算模型取自流线型车头拱面,配一段车身。
由于结构基本对称,沿纵对称面取一半建模。
对任何CFD软件包,很难创造出三维流线型车头这样复杂的几何模型,所以首先在Ideas中建立模型,然后用STEP格式传到Fluent中进行仿真计算。
计算区域选为70m!
10m!
10m。
计算网格的划分采用非结构化网格,对车体近壁层的网格细化,远离车体的网格采用稀疏网格,以减少计算量和加快收敛速度。
列车表面划分三角形网格,空间采用四面体网格,空间体单元约为140万。
因高速列车速度为270km/h,故明线上列车运行绕流问题一般采用三维黏性不可压缩定常流动处理。
4.2计算结果分析
图4
是列车外流场空气流动轨迹。
在列车尾部,可以看到一个负压区域的尾流场。
在尾流的横截面上沿列车宽度方向对称地产生一对
涡,其涡心较低,在列车纵向对称面上涡强最大,沿宽度方向涡强减弱。
随着尾流发展,涡强减弱,并向地面靠近。
在距列车尾部7m处,涡流的影响减小。
因此,此车的安全退避距离为距列车尾部7m左右。
研究列车尾流场,可以为高速列车的安全退避距离提供依据。
计算得出此车的阻力系数Cx=0.262。
该值表明,车头拱面明显地减低了气动阻力,流线型拱面是合理的。
本次计算为后续的有限元计算提供了车头拱面压力分布。
经过以上分析,可得出如下结论:
(1此流线型车头的流线型拱面是合理的,明显地减低了气动阻力。
(2此列车尾部的安全退避距离大约是7m。
5列车水箱晃动问题的仿真计算
列车上的水箱在列车启动、刹车过程中,水的晃动问题显得越来越突出。
以往,人们通常是根据一些经验进行水箱的设计工作,曾发生过水箱失稳、起皱现
象。
确切地了解水箱在列车运行过程中所受水的冲击力,是一个重要问题。
采用Fluent软件可以对列车刹车过程中,水箱中的水晃动问题进行模拟仿真。
5.1建模分析
水箱模型如图5所示,外形是椭圆柱形,内设3个防波板把水箱分成4个容腔。
每个防波板上开有若干通水孔。
水箱的尺寸是1760mm!
1600mm!
576mm,装有2/3容积的水。
本文模拟了当列车以120km/h的速度行驶在实际线路上制动工况下,水箱中水的晃动情况。
坐标建立在水箱上,因而在本身作变速运动的非惯性坐标系下,质量场中包括附加的惯性力。
采用Fluent的UDF接口,用C语言编写一段程序定义动量方程中的源项。
水箱中包含水和空气两种流体,所以是两相流动问题。
定义空气为初项,水为第二项。
采用两相流VOF(VolumeofFluid模型,在不可压缩的湍流非定常情况下,通过求解动量方程、方程模型以及计算在区域内每一种流体的体积百分比,
可以得出水箱中两种流体的运动情况。
5.2计算结果分析
图6为水箱纵向对称面上水的变化情况。
列车匀速行驶,水自由表面基本水平,随列车的刹车,水朝着列车运动方向涌动,首先碰撞到隔板及右端的水箱壁上,然后又往回弹,水的自由表面在不断的振荡变化,
防波板有效地阻止了水的晃动。
分析水箱壁的压力分布,找出最具破坏性的压力波,并提取该压力波下水对水箱各个部位的作用力,为屈曲分析的真实性奠定基础。
屈曲仿真分析的结果与水箱实际破坏情况[4]
基本吻合。
通过与实际对比,说明对水箱晃动的流场计算是可靠的。
67
6高速列车远场气动噪声数值分析
多年来,发达国家一直进行高速列车的气动噪声研究工作。
其中,荷兰在风洞中对高速车模型进行了气动噪声测量。
分析表明:
从列车模型上有尖角和陡峭的部位产生的涡都是很重要的噪声发生源。
结果建议:
通过优化列车顶部区域的形状,可减少大约10dB的噪声[5]。
本文应用Fluent软件和PISO算法对270km/h的二维流线型列车头部模型的远场气动噪声进行预测和分析。
6.1噪声计算
气体流动或者物体和气体相互作用引起气体的扰动而辐射的噪声,称为气动噪声。
模型选自某流线型车头的中心平面图,对列车的外表面和车下走行部分都进行了简化。
来特希尔的气动噪声模拟分为两部分:
第一步是声音由流体在连续媒介中流动诱导产生。
第二步是声音在外部的脉动源作用下,在静止的声媒介中传播的过程。
首先需要精确计算高速列车表面的压力脉动情况,然后提取列车表面的瞬时压力脉动变化,计算列车外流场中某观察点的声强、声压、观察点声强随时间的变化情况等有关的声学参数。
控制方程的求解采用非耦合求解方法。
在进行声学计算之前,采用大涡模拟计算的模型外流场要达到动态的稳定状态,计算所选用的时间步长非常重要。
时间步长的选取既需要在大涡计算中合适地求解湍流涡,也要在声场计算中能捕捉到声音波动,两者中的小者即定为计算所需要的时间步长。
在此,时间步长取12s,采用四边形网格,因为在大涡模拟计算中,四边形网格比三角形网格的计算误差小;
在时间积分方案上,选择具有至少二阶精度的CrankNicolson半隐式方案;
在空间离散格式上,选用二阶迎风格式,壁面采用壁面函数法。
6.2计算结果与分析
采用大涡模拟,可以在列车尾流区域中计算出一个清楚的旋涡。
当网格数比较粗糙,并且采用湍流模型定常流动计算时,得不到明显的湍流涡。
本文在距离车尾不同距离的4个观察点上计算声压随时间的变化。
观察点距离列车尾部越远,到达观察点的时间越来越长,声压越来越小,如表2所示。
表2观察点位置及总声压级
观察点x/my/mz/mSPL/dB
1201.50128.88
2251.50126.31
3301.50126.25
4351.50125.33
7结论
通过上述分析,得到以下结论:
(1本文对高速列车的空调通风系统、列车的外流场、列车的水箱晃动及列车的气动噪声等四个方面的空气流动问题进行数值仿真;
并将空调通风系统的仿真结果与试验结果比较;
水箱晃动的仿真结果与实际水箱使用情况进行对比,实验和计算结果符合较好。
(2通过CFD数值模拟可以对所研究问题多方面分析,反复修改模型,选用最佳方案满足设计要求。
这是一个十分经济的技术路线。
(3对不同领域的工程应用,采用不同的湍流模型,不同的算法,能够取得较好的计算效果。
(4数值仿真可以解决高速列车空气动力流场方面的许多问题,并能达到较高的精度要求,在时间和经济上提高高速列车设计研发的效率。
参考文献:
[1]王福军.计算流体动力学分析[M].北京:
清华大学出版社,2004.
[2]赵彬,等.用零方程湍流模型模拟通风空调室内的空气流动[J].清华大学学报(自然科学版,2001,41(10:
109113.
ZHAOBin,etal.Simulationofindoorairflowinventilatedroombyzeroequationturbulencemodel[J].JounalofTsinghuaUniversity(ScienceandTechnology,2001,41(10:
[3]徐趁肖.高速列车头型共映射建模理论研究[D].北京:
北方交通大学,2002,7.
[4]闫雪冬,等.客车水箱起皱机理及CAE性能仿真集成[J].中国机械工程,2003,14(Supp.:
15.
YANXuedong,etal.FailureMechanismofCrumplingWaterBoxofPassengerTrainandItsCAESimulationIntegration[J].ChinaMechanicalEngineering,2003,14(Supp.:
[5]SjoerdW.Rienstra,etal.AeroacousticsResearchinEurope:
theCeasascReporton1997Highlights[EB/OL]//AnAnnualSummaryofAeroacousticsResearchinthe8CEASCountries.(1997[2004].http:
//www.win.tue.nl(责任编辑江峰
68铁道学报第29卷
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- CFD 数值 仿真 高速 列车 设计 中的 应用
![提示](https://static.bdocx.com/images/bang_tan.gif)