列車(chē)空氣動(dòng)力學(xué)數值仿真研究

第23卷,第4期中國鐵道科學(xué)2002年8月CHINA RAILWAY SCIENCEAugust, 2002文章編號:1001-4632(2002)04013207博士論文摘要列車(chē)空氣動(dòng)力學(xué)數值仿真研究武青海鐵道科學(xué)研究院機輛所,北京100081)關(guān)鍵詞:列車(chē)空氣動(dòng)力學(xué);有限體積法;湍流模型;磁懸浮列車(chē)中圖分類(lèi)號:U260.111:V211.3文獻標識碼:A限差分法、有限單元法和有限體積法)1前言對于列車(chē)空氣動(dòng)力學(xué),與傳統的試驗研究方法相比數值仿真技術(shù)方法可以隨時(shí)隨意給定所需參數自從20世紀60年代,日本建成時(shí)速達210km進(jìn)行仿真計算,并利用計算結果選擇最佳設計,從的東海道新干線(xiàn)以來(lái),世界鐵路的高速化技術(shù)得到而節省人力、物力、財力,縮短產(chǎn)品設計周期。還迅速發(fā)展,使得鐵路運輸在經(jīng)濟發(fā)展中再次發(fā)揮了有工程中遇到的反問(wèn)題,諸如要求預先給定流動(dòng)結重要作用。2001年,法國建成了運營(yíng)速度達300km果,然后通過(guò)計算來(lái)確定物面形狀,而這些都是試h-的地中海線(xiàn),將世界鐵路的高速化進(jìn)程又推向驗所無(wú)法做到的?？梢哉f(shuō),對于解決工程問(wèn)題,仿個(gè)新的高潮。真計算可以幫助選型,減少設計費用,縮短設計周隨運行速度提高,空氣阻力等列車(chē)空氣動(dòng)力學(xué)期。并且對于某些特殊問(wèn)題,仿真計算還可以代替問(wèn)題日益受到重視。因此為實(shí)現高速鐵路的安全、試驗。因此本文采用了數值仿真方法來(lái)研究列車(chē)空舒適、降低能耗并滿(mǎn)足環(huán)境保護的要求,世界各國氣動(dòng)力學(xué)問(wèn)題。逐漸開(kāi)展了對列車(chē)空氣動(dòng)力學(xué)問(wèn)題的研究,并在此這其中有限體積法在列車(chē)空氣動(dòng)力學(xué)數值仿真基礎上形成了一個(gè)獨立的學(xué)科——列車(chē)空氣動(dòng)力研究中有著(zhù)其它計算方法無(wú)法比擬的優(yōu)勢,并且也學(xué)是計算流體力學(xué)中目前最為流行的計算方法,因此列車(chē)空氣動(dòng)力學(xué)問(wèn)題可以分為以下幾類(lèi):空氣本文采用了有限體積法進(jìn)行列車(chē)空氣動(dòng)力學(xué)數值仿阻力;側向風(fēng)的空氣動(dòng)力學(xué)特性;會(huì )車(chē)時(shí)作用在列真計算。車(chē)上的氣動(dòng)壓力;列車(chē)風(fēng);隧道內的壓力波動(dòng);隧道出口微氣壓波;車(chē)體表面湍流邊界層的壓力變3理論基礎及數值計算方法動(dòng)我國在此方面的研究起步較晚,90年代初才正列車(chē)是一個(gè)長(cháng)寬比很大的細長(cháng)物體,形狀比較式開(kāi)始,起初主要是進(jìn)行實(shí)車(chē)試驗及風(fēng)洞試驗;數復雜,而且近地運動(dòng),所以繞列車(chē)的空氣流動(dòng)是完值仿真技術(shù)盡管與世界先進(jìn)水平相差很大,但近幾全三維的流動(dòng)。在列車(chē)尾部會(huì )出現一對很強的尾年發(fā)展速度也很迅速渦。通過(guò)列車(chē)外部流場(chǎng)的基本特征以及流場(chǎng)雷諾數來(lái)判斷其流動(dòng)的形態(tài),可以知道列車(chē)外流場(chǎng)為湍流2列車(chē)空氣動(dòng)力學(xué)主要研究方法流場(chǎng)I凵中國煤化工論研究及實(shí)際應用中試驗研究方法(風(fēng)洞試驗、水洞試驗、電模擬都已CNMHG校滿(mǎn)意的結果特別是試驗實(shí)車(chē)試驗)和數值仿真研究方法(面元法、有對外部流動(dòng)的研究較多。因此本文選用kε雙方程湍收稿日期:200205作者簡(jiǎn)介方據"3-),男,河南安陽(yáng)人,博士第4期博士論文摘要133流模型來(lái)研究列車(chē)空氣動(dòng)力學(xué)問(wèn)題?！鱵23G,△xYH2k-ε雙方程湍流模型的求解方程組由以下基本Re△y2ayRe1△z方程組成連續方程變換后得到最終形式dU dF dg dHdFdG雷諾方程H△d2+(△xm△ym△zm…)0+pax,=p、少當a3,y分別取值小于或遠小于m時(shí)所采用的網(wǎng)格和差分格式才能比較正確地計入各方向的粘性貢獻。因此在α3,γ分別取值m時(shí)就可以得到dx: a湍動(dòng)能k方程臨界網(wǎng)格間距值Δx*、Δy*、△z*。Ig hIgErak, duiau'. dg△z=gRe湍流耗散率ε方程5計算域的選取2u dump在仿真計算中,一般采取有限計算域來(lái)代替無(wú)2限計算域,計算域的大小直接影響到計算結果的可dx dx: d信程度。但是對于不同類(lèi)型計算對象,其計算域的ut;/u’;du’;du大小不能一概而論,所以需要具體問(wèn)題具體對待。ax. ax adx: d本論文中的計算域的選取采取了試算分析的方法確au';o’oua-u'定。在對于計算域的選取中,試算方法不失為一種有效的方法數值計算方法采用有限體積法中常用的SIPE( Semi-Implicit Method for Pressure-Linked Equa6移動(dòng)網(wǎng)格技術(shù)ions)求解離散方程組。在列車(chē)交會(huì )和過(guò)隧道的瞬態(tài)計算中,存在計算4貼體網(wǎng)格間距的確定對象之間的相對運動(dòng),利用移動(dòng)網(wǎng)格來(lái)模擬物體之間的相對運動(dòng),可以很好地進(jìn)行處理移動(dòng)網(wǎng)格界面三維流動(dòng)無(wú)量綱化的NS方程之間的邊界條件。dU dF aG dH在移動(dòng)網(wǎng)格的一些實(shí)際應用中,計算域的大F dG aH在計算過(guò)程中有很大變化,此時(shí)在瞬態(tài)計算中可通Rell axt ay過(guò)單元移除或添加來(lái)解決上述問(wèn)題,以保證計算域階精度的差分格式求解后,方程成為的完整性。dU, aF(dFr aGat t ax tRe7以 STARCD為基礎的仿真程序dHd2/+△xm△ym△m灬…)中國煤化工CNMHG件作為核心進(jìn)行列車(chē)a grewIg re外部流功計算,世一個(gè)用軟件應用于具體專(zhuān)業(yè)時(shí)方程經(jīng)需要與其他程序一起使用以及進(jìn)行必要的二次開(kāi)發(fā)才能更好地為本專(zhuān)業(yè)服務(wù)aU3 F aG dH△xdF2o方方數據y+z=Re△xx+134國鐵道科學(xué)第23卷開(kāi)始型列車(chē)頭尾的壓力差,且前者的高壓區受力區域面給定列車(chē)模型主要外形尺寸參積要大于后者,因此鈍頭型列車(chē)的頭尾壓差阻力比流線(xiàn)型列車(chē)的頭尾壓差阻力大得多在CAD圖形軟件中建立列車(chē)模型的主要截面數據(4)鈍頭型列車(chē)總的運行空氣阻力要大于流線(xiàn)通過(guò)DF格式轉換,進(jìn)入UCAD程序建立三維列車(chē)模型型列車(chē)的空氣阻力。要降低列車(chē)的運行空氣阻力,輸出列車(chē)模型表面特征數據,生成格式轉換文件須從優(yōu)化設計列車(chē)頭尾形狀入手,好的流線(xiàn)型列車(chē)形狀是降低阻力的主要有效手段進(jìn)入 STARCD前處理程序將計算域合理分成若干子域,然后劃分計算網(wǎng)格10列車(chē)風(fēng)計算結果進(jìn)入 STARCD后處理程序采用數值仿真方法計算分析了流線(xiàn)型和鈍頭型云圖,編制程序從 STARCD中兩類(lèi)車(chē)體周?chē)熊?chē)風(fēng)基本特征(1)以神州號動(dòng)車(chē)組為代表的流線(xiàn)型車(chē)體周?chē)幹屏熊?chē)空氣動(dòng)力性能結果數據處理程序列車(chē)風(fēng)基本特征:列車(chē)風(fēng)合速度在車(chē)頭車(chē)尾各出現通過(guò)繪制圖表程序生成圖表個(gè)峰值,隨距地面高度增加,尾部波動(dòng)要大于頭結束部的波動(dòng);車(chē)頭附近壓力場(chǎng)出現一個(gè)由正到負的壓經(jīng)過(guò)上述這樣對ARCD軟件的修正以及二次力波,且最大正壓力出現在車(chē)頭的前方;車(chē)尾附近開(kāi)發(fā)后,可以讓工程設計人員很方便地應用該軟件有一個(gè)由負到正的壓力波,波峰出現在列車(chē)尾端的進(jìn)行產(chǎn)品的前期開(kāi)發(fā)和研究后方。(2)以Ss8列車(chē)編組為代表的鈍頭型車(chē)體周?chē)?計算模型選取與簡(jiǎn)化列車(chē)風(fēng)基本特征:列車(chē)風(fēng)合速度在頭部的波動(dòng)要遠大于尾部的波動(dòng);壓力場(chǎng)的分布規律也驗證了這為不失一般性,與實(shí)際情況相結合,選取了六點(diǎn),說(shuō)明頭部擾動(dòng)產(chǎn)生的危險程度要大于尾部。種車(chē)型:DF4D,DF11,SS8,神州號,藍箭號,270通過(guò)分析列車(chē)風(fēng)速度場(chǎng)研究分析,為確定我國kmh電動(dòng)車(chē)組鐵路沿線(xiàn)安全距離提供了理論和科學(xué)依據。計算模型建立應遵循兩條原則物理的真實(shí)性、數學(xué)計算的可行性11列車(chē)過(guò)隧道和交會(huì )時(shí)的數值傷真就以下幾方面對列車(chē)模型進(jìn)行簡(jiǎn)化:結果采用兩動(dòng)一拖模擬列車(chē)編組、簡(jiǎn)化列車(chē)表面、簡(jiǎn)化車(chē)體底部等。1)采用簡(jiǎn)化模型分析了列車(chē)進(jìn)入隧道到出隧道期間列車(chē)車(chē)身上的壓力變化和隧道內的壓力波9列車(chē)外部流場(chǎng)計算結果動(dòng)特征。若計算機工作能力許可,經(jīng)精確計算,可為隧道凈空面積和橫截面形狀設計時(shí)的空氣動(dòng)力學(xué)(1)對270kmh-電動(dòng)車(chē)組在橫風(fēng)工況下的題提供參考依據風(fēng)洞試驗結果與仿真計算結果進(jìn)行了比較,兩者吻(2)初步探討了列車(chē)交會(huì )時(shí),兩列車(chē)內側車(chē)身合很好,說(shuō)明采用 STARCD軟件能夠較好地描述和側壁承受了車(chē)頭車(chē)尾通過(guò)時(shí)的兩次壓力瞬變,但是研究鐵路列車(chē)空氣動(dòng)力學(xué)性能。分析了六個(gè)模型車(chē)對于兩車(chē)之間流場(chǎng)復雜的流動(dòng)杋理尚未有定論。身縱向對稱(chēng)面上的壓力分布規律,從而可以定量比較流線(xiàn)型列車(chē)與鈍頭型列車(chē)的區別。12(2)運營(yíng)速度在120kmh-1-160kmh-1范圍研究H中國煤化工動(dòng)力學(xué)數值仿真CNMHG內的列車(chē),若采用流線(xiàn)型車(chē)體,則與普通列車(chē)相比,可以使空氣阻力降低70%以上,從而達到降低〔1)對流場(chǎng)而言,不外包形狀的磁懸浮列車(chē)能耗的目的外部流場(chǎng)比外包形狀的要復雜,阻力系數也有所增3)鈍灰殞車(chē)車(chē)頭與車(chē)尾的壓力差大于流線(xiàn)大,所以流線(xiàn)型的外包形狀比較有利于列車(chē)周?chē)鷼獾?期博士論文摘要135體的流動(dòng)。但磁懸浮列車(chē)與線(xiàn)路是否外包對總旳運(3)高速鐵路隧道斷面形狀及岀入口形狀的研行空氣阻力的影響不大。究(2)在對列車(chē)模型簡(jiǎn)化的基礎上,磁懸浮列車(chē)4)隧道內部通風(fēng)以及散熱等問(wèn)題的運行空氣阻力計算結果比輪軌式列車(chē)大13%左5)列車(chē)交會(huì )時(shí)兩列車(chē)中間區域流場(chǎng)的空氣流右。因此,相同頭型和迎風(fēng)截面外形的輪軌式列車(chē)動(dòng)機理研究以及線(xiàn)間距的確定的運行空氣阻力要小于磁懸浮列車(chē)。(6)氣動(dòng)噪聲問(wèn)題13工作展望希望通過(guò)對列車(chē)空氣動(dòng)力學(xué)的數值仿真研究,能夠建立“數值風(fēng)洞”,使其在工程設計中得到廣(1)底部轉向架的近似模擬泛應用,加速發(fā)展我國的高速鐵路技術(shù)。(2)車(chē)內通風(fēng)問(wèn)題a Study on Numerical Simulation of Train Aerodynamicswu Qing-HaiChina Academy of Railway Sciences, Beijing 100081, Chinasure shift in tunnel the micro-pressure at tunnel exit, and1 Introductionthe pressure shift in the turbulence layer of train surfaceThis study in China is behind the world, startedAfter the first Shinkansen train in Japan run success- from actual train test and wind tunnel test at the early dayfully at 210 knt h-I in 1 960 s, the high-speed railway of 1990 s. Recently the numerical simulation techniquetechnique has been developed quickly all over the world has been developed rapidly although there is a big distanceand made railway transportation play a important role in e- between China and the worldconomic development again. In 2001, the Med Line wasbuilt in France so the high-speed technique of the world 2 Main study methods for train aerody-railway had reached a new steWith raising the train running speed, train aero-dynamIcsnamics problems have been regarded increasingly. Thestudy on train aerodynamics problems have been developedTest study method( wind tunnel test, water tunneltest electronic simulation test and actual train test ), andand comfortable of high-speed train reduce the powers numerical simulation method panel method, finite differconsume, and environmental protection. And based on encethese an independent subject, train aerodynamics has metheH中國煤化工 ethod and finite volCNMHGappearedFor train aerodynamics, compared with traditionalThe train aerodynamics can be classified as belotest study method, the simulation calculation can be madeAir drag side direction wind the air pressure on at any time when setting the required parameters. Thethe train wheh erosing, train induced airflow, the pres- best design will be selected according to the result, and it136國鐵道科學(xué)第23卷an cut down the consumption in manpower money andlL:antage of simulation method can be showed also throughTurbulence kinetic energy h equationthe anti-question in project. Such as the flow result areknown in advance then the surface shape of product willbe confirmed by calculation, while the target can not be2reached by test method. That is to say, for project probTurbulence dissipation rate E equationlems the shape can be selected, the design expense canbe reduced and the design period can be shortened alsojdx, dthrough simulation calculation And it can substitute fortest in some special problems. So the numerical simulationa∈)-2t'ujaxmethod is used to research the train aerodynamics proldu;/du’;o7u't0lems in this paperdCompared with other calculation methods the finitedu';du';Ou’volume method( FVM) has more advantage in train aerodxr dx: d.dynamics research. And it is one of the most popular calThe numerical calculation method adopts SIMPLEculation methods in CFD today. Therefore we adopt FVM( Semi-Implicit Method for Pressure-Linked Equations )ofto make the numerical simulation calculation for train FVM to solve the discredited equation setaerodynamIcs3 Theory foundation and numerical calcu- meshes near walllation method3D zero dimension N-S equationdU dF dG dHTrain has a slender body with large ratio of length towidth and complicated shape. And train runs close to thedF. dG. dHground. So the airflow around train is 3D completely. ARe, ox+ apair of strong tail vortex appears behind train. AccordinAfter m rank precision difference form solution,eto the basic characteristic of the airflow around train and quation becomesthe flow reynolds number we can know that the flow fieldau dF aG dH 1 dF aGaz Rearound train is a turbulence fieldIn theory research and application k-E turbulencea2)+△rAym△”灬)model has reached a mature degree and gotten many morelureIg Resatisfied results especially for exterior flow field. So weadopt the k-e turbulence model to study the train aerodSubstitute in the above equation and rearrangeU dF dg aH△xdF2The solution equation set of h-E turbulence modelRe1△xax+consists of the following aquatic△y2dG,△zH中國煤化工02+Continuity equationHCNMHGWe get the final expression1+3F+3+dFGa片數據a- op△z+O△xAym△znr…)第4期博士論文摘要137When a,B and y are less than m or more respectivelthe adopted mesh and difference form can include the viof traincosity contribution of each direction correctlyTherefore as a,B and are equal to m we canEstablish the main section data of train model in CAD softwareget the critical distance between meshes Ax,AyEstablish 3D train model in UGCAD by format transformOutput the surface data of train model,grelg reand create transform file withnto sTaRCd preprocedurelIgErColculation field divided into some sub-fields and5 Selection for the calculation fieldInto STARCDIn the simulation calculation people substitute theOutput flow field and pdistribution graphlimited calculation field for the infinite calculation field ingotten the result dataprogrammMinGgeneral. The size of the calculation field influences theDesign program to treat with the result data of train aerodynamicscreditable degree of calculation result directly. But for thedifferent type models their calculation field cannot beChart generated by prconsidered similarly. So specific problem should be treated especially. The calculation fields in thAfter above secondary development project plannerfirmed by the pilot calculation method. The pilot calculacan use STARCD software easily to make prophase develtionthod may be aneffective way for the calculation opment and research for productfield selection8 Selection and simplification of calcula-6 Sliding mesh techniquetion modelIn order to ensure generality six type trains are seIn the transient calculation of train crossing and trainlected: DF4D, DF11, S$8, Shenzhou, Lanjian 270 kmpassing a tunnel there is a relative movement between thehpower callculation objects. It can deal with the bounThe calculation model's establishment should followtion between the movement mesh interfaces nicely thoughIn some applications of sliding mesh, the size of cal- mathematical calculatio ity of physics and feasibility ofthe relative movement simulated by sliding meshculation field may change largely at the calculation proThe train models are simplified at the followincess. At this time we can use cell removal and addition topects according the principles: two power cars and onesolve the above problem, in order to ensure the integrality trail car used to simulate train set simplification of trainof calculation fieldsurface and the train body bottom9 The calculation result of train exterior7 Simulation program based onflow fieldSTARCD(1) After comparison, the calculation results of 270The calculation for the exterior flow around train is km.YH中國煤化工 condition of crosswind fitmade where the fluid software staRCd is considered as withC Gnel very well. It meansthe core. When general software is used to certain field, the software STARCD can describe and study the trainit should be used with other programs together or made aerodynamics betterecondary development in order to serve best for this spe-(2) For train which run at 120 kml h- -160 kth-, if a streamline body is adopted its air drag will be138中國鐵道科學(xué)第23卷less 70 than that of train with general body in order to passing tunnel and the pressure fluctuation in tunnel is anreach the target reducing power's consumeeIs. The study will provide ref-(3)The pressure differences of nose for blunt nose erence basement for the aerodynamics problems during detrain are larger than that of streamline train and the same signing tunnel area and section shape by accurate calculaas the end. The high-pressure size of the former is larger tion if the work ability of computer allowsthan that of the latter also. So the pressure differences be-(2) The sidewall of train body between two trainstween nose and end of blunt nose train is larger than that endures the pressure transient twice during trainf streamline trainBut the complicated flow mechanism between two trains is(4) The total running air drag of blunt nose train is not yet clear up to nowger than that of streamline train In order to reduce the12 Numerical simulation study of maglevgood streamline shape is main effective way to reduce airaerodynamIcs(1) For the flow field, no wrap maglev is compli-10 The calculation result of train inducedcated correspondingly compared with wrap maglev its airdrag coefficient is much larger thanof wraairflowSo the maglev with streamline wrap shape is propitious toThe basic characteristic of train induced airflow is the air flow around train but it has small influence foranalyzed with numerical simulation methodtotal air drag whether the shape between maglev and rail(1)The basic characteristic of train induced airflowaround the streamline train represented by Shenzhouhao(2)Based on the simplified model the calculationpower car set: the peak value of the resultant velocity Vxy result of maglev air drag is more than that of wheel trainappears twice at the nose and end respectively. The Vay about 13 percent. Therefore the air drag of wheel train isshift of the end is larger than that of the nose when the less than that of maglev with the same nose type and winddistances to grourise gradually. One press wave froward sectonpositive to negative appears near the train nose and themost pressure lies in front of train. Another press wave13 Work prospectfrom negative to positive appears near the train end and(1) Approximate simulation for the train bogiethe peak lies behind the train end(2) Ventilation problem inside train(2)The basic characteristic of train induced airflow(3) Study on the section shape entrance and exitaround the blunt train represented by SS8: The fluctuation shape of high-speed railway tunnelof the resultant velocity Iter than(4) Ventilation and heat emission in tunnelthat at the train end. The distribution law of pressure field(5) Study on the air flow mechanism duringalidates this point also. That is to say that the danger de- crossing and the interval between linesgree made by the nose fluctuation is larger than made by(6)Air noise problemthe train endset upBy analyzing the train-induced airflow this paper through studying train aerodynamics. Our high-speed railprovides theoretical and scientific basement for the railway way technique can be developed fast with the applicationsafety distance in China11 The simulation result of train passing KeyI Vpnd tunnel" widely in engineering design中國煤化工FVM TurbuleCNMHtunnel and crossingmoc責任編輯劉衛華)(1) The pressure change on train body during train