• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    空調(diào)安裝方式優(yōu)化對列車氣動特性影響研究

    2018-10-31 08:13:32洪琪琛楊明智丁暢
    關(guān)鍵詞:尾車頭車升力

    洪琪琛,楊明智,丁暢

    ?

    空調(diào)安裝方式優(yōu)化對列車氣動特性影響研究

    洪琪琛,楊明智,丁暢

    (中南大學(xué) 交通運(yùn)輸工程學(xué)院,軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410075)

    基于三維、不可壓和定常N-S方程的數(shù)值模擬計(jì)算方法,從空調(diào)機(jī)組的安裝位置和高度等2個方向著手,分析不同方案下列車運(yùn)行時的氣動特性。研究結(jié)果表明:當(dāng)頭車空調(diào)位置遠(yuǎn)離車頭、尾車空調(diào)位置靠近車尾時列車總阻力最大,較之其他方案大2.66%左右;空調(diào)阻力隨著車體高度增加而逐漸減小,當(dāng)車體增加高度達(dá)到350 mm時,3節(jié)車總阻力下降7.55%,但與此同時,列車在橫風(fēng)環(huán)境下運(yùn)行穩(wěn)定性變差。

    高速列車;數(shù)值計(jì)算;外形優(yōu)化;空調(diào)機(jī)組;氣動特性

    高速鐵路作為國內(nèi)外最重要的交通運(yùn)輸方式,具有運(yùn)輸量大、環(huán)保經(jīng)濟(jì)和安全性高等優(yōu)點(diǎn)。由于高速列車是在地面上高速運(yùn)行的細(xì)長體,導(dǎo)致一系列空氣動力學(xué)問題。當(dāng)列車運(yùn)行速度超過200km/h時,氣動阻力占總阻力70%左右,且隨著列車運(yùn)行速度的提高,其氣動阻力占比越高[1?3]。為了減少列車運(yùn)行時空氣阻力,國內(nèi)外學(xué)者進(jìn)行了大量關(guān)于列車整體外形優(yōu)化設(shè)計(jì)的研究并取得顯著成果[4?6],隨著列車整體外形優(yōu)化進(jìn)入瓶頸期,人們開始探究列車局部設(shè)備優(yōu)化來實(shí)現(xiàn)減阻。胡勇等[7]采用數(shù)值模擬和實(shí)車實(shí)驗(yàn)相結(jié)合的方法研究空調(diào)進(jìn)出氣口和風(fēng)道內(nèi)的氣動性能隨車速變化的規(guī)律,得出不同車速和各方案進(jìn)出氣口的風(fēng)量變化情況,對空調(diào)設(shè)備的選取和設(shè)計(jì)提供建議。肖云華等[8]采用數(shù)值計(jì)算的方式研究了8節(jié)車編組的動車組明線工況下運(yùn)行時空調(diào)表面的壓力沿車長和車寬方向的分布規(guī)律,并對空調(diào)裝置的安裝位置給出了意見。國內(nèi)對于列車空調(diào)裝置優(yōu)化的研究主要集中于空調(diào)進(jìn)、排風(fēng)道的外形設(shè)計(jì)和風(fēng)機(jī)外形優(yōu)化,主要研究目的是空調(diào)裝置對列車內(nèi)流場和車廂內(nèi)舒適度的影響[9?11],對空調(diào)裝置的優(yōu)化設(shè)計(jì)實(shí)現(xiàn)減阻和提升列車在橫風(fēng)環(huán)境下運(yùn)行穩(wěn)定性的研究較少。在安裝空調(diào)設(shè)備時,考慮其設(shè)備自身的重量,列車空調(diào)的冷凝器和散熱器一般都放置在列車的頂部,空調(diào)裝置作為列車頂部的凸起部件,在列車實(shí)際運(yùn)行中,容易在其安裝區(qū)域出現(xiàn)空氣分流現(xiàn)象,從而加大列車運(yùn)行時的空氣阻力。本文采用數(shù)值模擬方法,分析空調(diào)不同的安裝位置和安裝高度對列車運(yùn)行時阻力的影響,從而找到一個最合理的優(yōu)化方案,為列車空調(diào)機(jī)組的設(shè)計(jì)及安裝提供一定的參考依據(jù)。

    1 控制方程

    本文研究中的車速為350 km/h,馬赫數(shù)小于0.3,因而引起的空氣流動是不可壓縮的湍流流動,針對不可壓縮流的基本控制方程進(jìn)行求解。

    根據(jù)質(zhì)量守恒定律,單位時間內(nèi)流入流體微元的質(zhì)量和流出流體微元的質(zhì)量相等,可導(dǎo)出不可壓縮流體運(yùn)動的質(zhì)量守恒方程[12?13]:

    根據(jù)動量守恒定律:作用于流體微元上所有作用力的和與微元體中流體動量增長率相等,可得到3個方向的動量守恒方程:

    式中:u為流場速度,代表,和3個坐標(biāo)方向的速度分量;xx為3個坐標(biāo),代表,和3個坐標(biāo)分量;為流體密度,在不可壓縮流中為常數(shù);為壓力;為空氣動力黏度。

    -雙方程模型在列車周圍空氣流場流動的湍流數(shù)值模擬中應(yīng)用最為廣泛,控制方程組見文 獻(xiàn)[14?15]。

    2 計(jì)算模型

    2.1 動車模型

    我國高速列車通常采用8節(jié)車編組,由于列車中部截面保持不變,縮短后的模型不會改變列車流場結(jié)構(gòu)的基本特征,為了提高計(jì)算效率,本文的動車模型均采用3節(jié)車編組,即頭車+中間車+尾車,頭車及尾車長度為26.5 m,中間車長度為25 m,車高3.7 m,車寬3.4 m。除空調(diào)裝置及轉(zhuǎn)向架外,車體表面一些細(xì)小部件均不考慮。圖1為原始動車 模型。

    圖1 高速列車原始模型

    為了研究空調(diào)安裝位置對整車和其他設(shè)備的氣動性能的影響,采取4種安裝位置進(jìn)行研究:以位置1為基準(zhǔn),位置2為頭車車頂?shù)目照{(diào)向車頭方向移動了6 800 mm,位置3為頭車空調(diào)和尾車空調(diào)分別向車頭方向和車尾方向移動了6 800 mm;位置4為尾車空調(diào)向車尾方向移動了6 800 mm??照{(diào)安裝位置如圖2所示。

    (a) 位置1;(b) 位置2;(c) 位置3;(d) 位置4

    在保證車廂內(nèi)部的空間前提下,列車高度是有最低限度的,再加上頂部進(jìn)氣排氣通道和電路設(shè)備的存在,幾乎沒有了埋入空調(diào)裝置的空間,因此沒有辦法通過下沉空調(diào)高度以達(dá)到減阻目的,而只能考慮將車身加高,以達(dá)到空調(diào)高度降低的效果?;谝陨纤悸?,本文設(shè)計(jì)3種高度的列車模型:1) 車高為3 700 mm,此時空調(diào)導(dǎo)流罩完全露出車頂,車體增加的高度=0;2) 車高為3 860 mm,此時空調(diào)導(dǎo)流罩有一部分已經(jīng)埋入車頂,車體增加的高度=160;3) 車高為4 050 mm,此時空調(diào)導(dǎo)流罩已經(jīng)完全埋入車頂,車體增加的高度=350。具體模型見圖3。

    (a) h=0 mm;(b) h=160 mm;(c) h=350 mm

    2.2 計(jì)算區(qū)域、網(wǎng)格及邊界條件

    據(jù)歐洲標(biāo)準(zhǔn)BSEN 14067,計(jì)算域邊界不應(yīng)干擾車輛周圍流場,計(jì)算域應(yīng)保證流場充分發(fā)展:計(jì)算區(qū)域上游不小于8倍特征高度,計(jì)算域下游不小于16倍特征高度[16?17]。特征高度定義為列車頂面距離地面的距離,本文所有工況均為平地運(yùn)行,因此特征高度取車高=3.7 m。列車車輪底面距軌面高度為0.2 m,為便于設(shè)定入口及出口邊界條件,計(jì)算區(qū)域上游入口處和下游出口處都盡量遠(yuǎn)離車體,避免受來流或者尾流的影響。如圖4(a)所示為列車在明線運(yùn)行時的計(jì)算區(qū)域,區(qū)域前端面定義為速度入口,距離動車組鼻尖點(diǎn)為27,速度與列車運(yùn)行速度相等,均為97.22 m/s,方向相反;為消除地板附面層的影響,地面給定與速度入口一致的滑移邊界條件;區(qū)域尾端距離動車組尾部60,給定壓力出口邊界條件,靜壓為0 Pa,兩側(cè)面及頂面定義為對稱面邊界條件,車體表面為無滑移壁面邊界條件。圖4(b)為列車在橫風(fēng)環(huán)境下運(yùn)行時的計(jì)算區(qū)域,在橫風(fēng)環(huán)境下除了有與列車進(jìn)行相對運(yùn)動形成的流場之外,還有自然環(huán)境導(dǎo)致的橫風(fēng),風(fēng)速為20 m/s,方向沿軸正方向流動,入口面IJNM及面JKON均定義為速度入口邊界條件,面LKOP及面ILPM均定義為壓力出口邊界條件,靜壓為0 Pa,車體表面采用無滑移壁面邊界條件。

    (a) 明線環(huán)境下計(jì)算域;(b) 橫風(fēng)環(huán)境下計(jì)算域

    在網(wǎng)格離散過程中,對列車周圍的網(wǎng)格進(jìn)行了加密處理以保證計(jì)算的精確,同時在列車尾部后方流場由于尾渦的存在,該區(qū)域也進(jìn)行了加密處理,除轉(zhuǎn)向架和受電弓等復(fù)雜設(shè)備采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行離散外,車體表面的物面網(wǎng)格為四面體結(jié)構(gòu)化網(wǎng)格,具體網(wǎng)格圖見圖5。明線工況下網(wǎng)格總數(shù)為5 500萬,橫風(fēng)工況下網(wǎng)格總數(shù)為5 800萬。

    圖5 計(jì)算網(wǎng)格

    2.3 數(shù)值方法驗(yàn)證

    為了驗(yàn)證數(shù)值模擬計(jì)算方法的準(zhǔn)確性,在中國空氣動力研究與發(fā)展中心的8 m×6 m的風(fēng)洞進(jìn)行了風(fēng)洞試驗(yàn),用1:8的縮比模型的試驗(yàn)結(jié)果與數(shù)值模擬計(jì)算的結(jié)果進(jìn)行對比。試驗(yàn)?zāi)P统叽缗c計(jì)算模型尺寸一致,為3節(jié)車編組,為了便于分析,將各氣動力和力矩進(jìn)行無量綱化進(jìn)行對比。各氣動力系數(shù)和力矩系數(shù)分別定義為:

    =/(∞) (3)

    =/(∞) (4)

    =/(∞) (5)

    =/(∞) (6)

    式中:∞=0.5v2,為動壓;為氣流密度,通常取1.225 kg/m3;為列車橫截面面積,在本試驗(yàn)中=0.19 m2;為阻力;為側(cè)向力;為升力;為傾覆力矩;為橫向參考長度,在試驗(yàn)中取0.203 m,傾覆力矩的取矩點(diǎn)為列車背風(fēng)側(cè)輪軌接觸點(diǎn);為阻力系數(shù);為側(cè)向力系數(shù);為升力系數(shù);為傾覆力矩系數(shù)。

    數(shù)值計(jì)算模型完全依照風(fēng)洞試驗(yàn)?zāi)P徒?,地面邊界條件設(shè)置為固定壁面,風(fēng)洞試驗(yàn)結(jié)果與數(shù)值計(jì)算結(jié)果對比如表1所示,數(shù)值計(jì)算與風(fēng)洞試驗(yàn)整車阻力相差4.21%,數(shù)值計(jì)算準(zhǔn)確可靠,滿足工程應(yīng)用的要求。

    表1 風(fēng)洞試驗(yàn)與數(shù)值計(jì)算結(jié)果對比

    3 計(jì)算結(jié)果與分析

    3.1 列車表面壓力分布

    圖6為列車明線工況下以350 km/h運(yùn)行,4種不同空調(diào)機(jī)組安裝位置,列車鼻尖點(diǎn)頂部區(qū)域縱截面車體表面壓力系數(shù)沿車長方向分布情況。觀察圖中壓力分布曲線,沿車長方向,頭車空調(diào)壓力變化要明顯大于中車和尾車空調(diào)。位置2和位置3在頭車空調(diào)迎風(fēng)側(cè)和背風(fēng)側(cè)曲面過渡處的負(fù)壓絕對值均大于位置1和位置4;位置3和位置4在尾車空調(diào)曲面過渡處的負(fù)壓絕對值要大于位置1和位置2。因此可以得出,靠近車頭和車尾的空調(diào)裝置的最大負(fù)壓絕對值要大于遠(yuǎn)離車頭和車尾的空調(diào)裝置的最大負(fù)壓絕對值。

    圖6 列車縱截面壓力系數(shù)

    圖7為不同車高方案下空調(diào)表面的壓力云圖,隨著車高的增加,空調(diào)的高度減小,空調(diào)裝置對流場的阻滯減弱,空調(diào)迎風(fēng)面的正壓減小,前面的研究中已經(jīng)發(fā)現(xiàn)空調(diào)的阻力主要受迎風(fēng)面正壓的影響,故空調(diào)的阻力減小??照{(diào)高度的降低也使得空調(diào)導(dǎo)流曲面對流場的加速效應(yīng)減弱,流場速度減小,使空調(diào)頂部曲面過渡處的最大負(fù)壓絕對值減小。繞過空調(diào)頂部的流場一部分會沿著空調(diào)背面下沖,直接作用于風(fēng)擋,空調(diào)高度降低后,下沖作用于風(fēng)擋的氣流減少,故風(fēng)擋的正壓減小。

    當(dāng)車體增加的高度達(dá)到350 mm時,空調(diào)已經(jīng)完全埋入車頂,車頂?shù)牧鲌龅姆植际志鶆?,但由于列車高度的增加,如圖7所示,車頭的流線型拉長,過渡曲面處拐點(diǎn)后移,對流場的加速效應(yīng)增強(qiáng),故車頭曲面過渡處的負(fù)壓面積會增大,頭車升力也會因此改變,且最大負(fù)壓的位置較=0和=160時更靠后。

    (a) h=0 mm;(b) h=160 mm;(c) h=350 mm

    3.2 氣動力結(jié)果分析

    表2為明線運(yùn)行時不同安裝位置的空調(diào)列車各部位阻力系數(shù)和升力系數(shù)。頭車空調(diào)在位置2和位置3時頭車的阻力系數(shù)大于在位置1和位置4時,說明頭車空調(diào)在靠近車頭位置時,頭車的阻力系數(shù)增大;尾車空調(diào)在位置1和位置2時尾車的阻力系數(shù)要大于在位置3和位置4時,說明尾車空調(diào)在靠近車尾位置時,尾車的阻力系數(shù)會減小。對整車來說,位置1的阻力系數(shù)稍大于位置2和位置3,位置4的阻力系數(shù)最大,比位置3的阻力系數(shù)大2.66%。頭車和中車的升力為負(fù),尾車升力為正,說明頭車的空調(diào)在靠近車頭時,升力系數(shù)會增大;尾車升力系數(shù)在位置3和位置4時更小,說明尾車升力系數(shù)在空調(diào)靠近車尾位置時會減小。

    表3為明線運(yùn)行時不同車高的空調(diào)列車各部位氣動力系數(shù)。隨著增加到350 mm,頭車、中車和尾車的阻力系數(shù)分別減少1.41%,11.30%和15.38%,列車的總阻力系數(shù)減少7.55%,由此可見,采取增加車高的方式對尾車的減阻效果最好,中間車次之,對頭車的減阻效果不明顯。

    表2 不同安裝位置下列車氣動力系數(shù)

    表3 不同車高下列車氣動力系數(shù)

    升力方面,從大小來看,頭車的升力系數(shù)變化了26.32% ,中車升力系數(shù)變化了26.53%,尾車的升力系數(shù)變化了27.08%,從方向來看,頭車和中車的升力方向?yàn)樨?fù),尾車方向?yàn)檎?。可見,列車高度的變化對列車升力的影響非常大,這可能也與明線運(yùn)行時升力基數(shù)很小有關(guān)。

    由于是通過列車車身高度增加,這樣增加了橫風(fēng)情況下的迎風(fēng)面積,必須分析其在橫風(fēng)環(huán)境下的側(cè)向力和傾覆力矩。表4為橫風(fēng)工況下列車各部位的氣動力系數(shù)隨車身高度變化規(guī)律。在同一車高的工況中,列車各部位的側(cè)向力系數(shù)和傾覆力矩系數(shù)大小關(guān)系為:頭車最大,尾車最小。隨著車身高度增加,列車的側(cè)向力系數(shù)和傾覆力矩系數(shù)都變大。頭車、中車和尾車的側(cè)向力系數(shù)分別增加了1.91%,9.07%和10.42%,受電弓的側(cè)向力系數(shù)變化不大。頭車、中車和尾車的傾覆力矩系數(shù)分別增加了4.44%,8.39%和5.99%,受電弓的傾覆力矩系數(shù)增加了1.54%。由此可知,隨著列車高度的增加,列車在橫風(fēng)下的穩(wěn)定性會變差。

    表4 橫風(fēng)下不同車高下列車氣動力系數(shù)

    3.3 流線分布

    圖8為頭車頂部的空調(diào)裝置縱剖面流線圖,列車周圍區(qū)域用速度進(jìn)行渲染,將車體顏色設(shè)置為白色以便區(qū)別于周圍速度流場。隨著列車高度的增加,空調(diào)高度則會降低,流場流線變得更加平順,流場在空調(diào)曲面過渡處的流速減小,風(fēng)擋上方區(qū)域流場流速加快。當(dāng)=0時,由于空調(diào)裝置的迎風(fēng)面對來流空氣的阻擋作用,相當(dāng)一部分的空氣被滯止,在滯止區(qū)域形成一個高壓區(qū),其余空氣沿著空調(diào)迎風(fēng)面向上,在空調(diào)迎風(fēng)面和水平面之間存在一個曲率較大的過渡區(qū)域,此處空氣速度達(dá)到最大,形成一個低壓區(qū)??諝獯撕笱刂照{(diào)機(jī)組水平面向后發(fā)展,在內(nèi)風(fēng)擋和外風(fēng)擋之間的區(qū)域產(chǎn)生了2個對稱的旋渦,由于該區(qū)域幾何尺度狹窄,氣流流向改變較小。

    隨著列車高度增加,風(fēng)擋空腔靠迎風(fēng)一側(cè)開始受流,旋渦消失,僅剩背風(fēng)一側(cè)的旋渦,且該旋渦的幾何尺寸隨著空調(diào)高度的減小而變大。

    (a) h=0 mm;(b) h=160 mm;(c) h=350 mm

    通過增加列車車身高度以達(dá)到降低空調(diào)機(jī)組高度的效果,可以較大地減小空調(diào)自身的氣動阻力,但是該方案增加了車體的高度,高速列車速度高、自重較輕、抗傾覆能力低,尤其是在強(qiáng)橫風(fēng)的作用下,列車氣動性能惡化,列車所受的空氣阻力和升力迅速增加,影響了列車的橫向穩(wěn)定性,嚴(yán)重時甚至將導(dǎo)致列車發(fā)生脫軌等事故,社會影響極大。因此必須對該方案的橫風(fēng)氣動特性進(jìn)行評估。

    圖9為頭車空調(diào)位置列車橫截面的流線圖,橫風(fēng)從列車上方繞過車頂后在列車背風(fēng)側(cè)形成了一個旋渦,可以明顯看出不同車高下漩渦的幾何尺度存在差異,由于車高的增加,空調(diào)機(jī)組的高度降低,空調(diào)導(dǎo)流罩對流場的加速效應(yīng)減弱,空調(diào)頂部流場流速減小。

    (a) h=0 mm;(b) h=160 mm;(c) h=350 mm

    4 結(jié)論

    1) 沿車長方向,頭車空調(diào)區(qū)域的壓力變化要明顯大于中車及尾車空調(diào)區(qū)域的壓力變化;靠近車頭和車尾的空調(diào)裝置的最大負(fù)壓絕對值要大于遠(yuǎn)離車頭和車尾的空調(diào)裝置的最大負(fù)壓絕對值。

    2) 當(dāng)頭車和尾車車頂?shù)目照{(diào)分別向車頭和車尾方向移動時,頭車的阻力系數(shù)會增大,尾車的阻力系數(shù)會減小,二者變化的程度不一樣,對整車來說,列車的阻力系數(shù)在位置4時最大,比其他方案的阻力系數(shù)大2.66%左右。

    3) 隨著車體高度的增加,頭車、中車和尾車的阻力系數(shù)均在減小,尾車受影響的程度最大,當(dāng)從0 mm增加到350 mm時,整車的阻力系數(shù)減少7.55%。

    4) 隨著車體高度的增加,列車在橫風(fēng)下的傾覆力矩系數(shù)增加,因此,要根據(jù)具體的運(yùn)行條件選取更優(yōu)方案,若列車在常年有大風(fēng)的環(huán)境下運(yùn)行,如新疆地區(qū),此時列車所受的氣動力會因受到自然風(fēng)的影響而出現(xiàn)明顯的波動,則選取=0 mm的方案,即不應(yīng)該對車身進(jìn)行加高;若列車在環(huán)境風(fēng)很小的地區(qū)運(yùn)行,則選取=350 mm的方案更優(yōu)。

    [1] 劉玉標(biāo), 張營營, 刑云林. 定常橫風(fēng)作用下高速列車的安全性分析[J]. 科學(xué)技術(shù)與工程, 2014, 14(8): 75?82. LIU Yubiao, ZHANG Yingying, XING Yunlin. Safety analysis of the high speed-train running under the steady aerodynamics forces of cross wind[J]. Science Technology and Engineering, 2014, 14(8): 75?82.

    [2] 梁習(xí)鋒, 田紅旗, 鄒建軍. 動力車縱向氣動力風(fēng)洞試驗(yàn)及數(shù)值模擬[J]. 國防科技大學(xué)學(xué)報(bào), 2003, 25(6): 101?105. LIANG Xifeng, TIAN Hongqi, ZOU Jianjun. The wind tunnel test and numerical simulation of longitudinal aerodynamic force of the traction car[J]. Journal of National University of Defense Technology, 2003, 25(6): 101?105.

    [3] 蔡國華. 高速客車模型氣動特性實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué), 2007, 21(4): 27?31. CAI Guohua. An experimental research of aerodynamic characteristics of the high speed passenger train model[J]. Journal of Experiments in Fluid Mechanics, 2007, 21(4): 27?31.

    [4] 馮志鵬. 高速列車氣動性能與外形設(shè)計(jì)[D]. 成都: 西南交通大學(xué), 2011: 1?88. FENG Zhipeng. Aerodynamic performance of high speed train with the train shape design[D]. Chengdu: Southwest Jiaotong University, 2011: 1?88.

    [5] 周丹,田紅旗,魯寨軍. 國產(chǎn)磁浮列車氣動外形的優(yōu)化[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2006, 37(3): 613?617. ZHOU dan, TIAN Hongq, LU Zhaijun. Optimization of aerodynamic shape for domestic maglev vehicle[J]. Journal of Central South University (Science and Technology), 2006, 37(3): 613?617.

    [6] 張?jiān)谥? 周丹. 不同頭部外形高速列車氣動性能風(fēng)洞試驗(yàn)研究[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 44(6): 2603?2608. ZHANG Zaizhong, ZHOU Dan. Wind tunnel experiment on aerodynamic characteristics of streamline head of high speed train with different head shapes[J]. Journal of Central South University (Science and Technology), 2013, 44(6): 2603?2608.

    [7] 胡勇. 高速列車橫向冷凝進(jìn)排風(fēng)道的氣動性能研究[D]. 成都: 西南交通大學(xué), 2013: 1?73. HU Yong. The study of the aerodynamic characteristics of the horizontal condensing air tunnel in high-speed train[D]. Chengdu: Southwest Jiaotong University, 2013: 1?73.

    [8] 肖云華, 李行, 熊小慧. 明線運(yùn)行時動車組空調(diào)裝置表面壓力分布數(shù)值分析[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2013, 10(3): 108?111. XIAO Yunhua, LI Hang, XIONG Xiaohui. Numerical simulation of surface pressure distribution of air conditioning equipment when EMU running in open air[J]. Journal of Railway Science and Engineering, 2013, 10(3): 108?111.

    [9] 歐陽仲志. 鐵路空調(diào)客車車內(nèi)設(shè)置溫度和控制分析[J]. 鐵道車輛, 2017, 55(1): 33?37. OUYANG Zhongzhi. Analysis of temperature setting and control inside air conditioned railway passenger cars[J]. Rolling Stock, 2017, 55(1): 33?37.

    [10] 張寶霞, 何忠韜. 鐵道車輛制冷與空氣調(diào)節(jié)[M]. 北京:中國鐵道出版社, 2005. ZHANG Baoxia, HE Zhongtao. Railway vehicle refrigeration and air conditioning[M]. Beijing: China Railway Press, 2005.

    [11] 吳俊云. 空氣客車廂體結(jié)構(gòu)傳熱數(shù)值分析[J]. 上海理工大學(xué)學(xué)報(bào), 2005, 27(4): 287?290. WU Junyun. Numerical heat transfer analysis of the air-conditioning bus envelope[J]. University of Shanghai for Science and Technology, 2005, 27(4): 287?290.

    [12] 田紅旗. 列車空氣動力學(xué)[M]. 北京: 中國鐵道出版社, 2007: 114?115. TIAN Hongqi. Train aerodynamics[M]. Beijing: China Railway Press, 2007: 114?115.

    [13] 陶文銓. 數(shù)值傳熱學(xué)[M]. 2版.西安: 西安交通大學(xué)出版社, 2013: 2?4. TAO Wenquan. Numerical heat transfer[M]. 2nd ed. Xi’an: Xi’an Jiaotong University Press, 2013: 2?4.

    [14] 吳望一. 流體力學(xué)[M]. 北京:北京大學(xué)出版社, 1998: 164?166. WU Wangyi. Fluent mechanics[M]. Beijing: Beijing University Press, 1998: 164?166.

    [15] 武青海, 周虹偉, 朱勇更. 高速列車湍流流場數(shù)值仿真計(jì)算探討[J]. 鐵道學(xué)報(bào), 2002, 24(3): 99?106. WU Qinghai, ZHOU Hongwei, ZHU Yonggeng. Research on numerical simulation calculation for high speed train in turbulence field[J]. Journal of China Railway Society, 2002, 24(3): 99?106.

    [16] Baker C J. The determination of topographical exposure factors for railway embankments[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1985, 21: 89?99.

    [17] GE/RT2542. Recommendations for determination of aerodynamic rolling moment coefficient[S].

    (編輯 陽麗霞)

    Study of the optimization of air conditioner installation mode’s effect on aerodynamic performance of high speed train

    HONG Qichen, YANG Mingzhi, DING Chang

    (Key Laboratory of Traffic Safety on Track, Ministry of Education, School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China)

    Based on steady N-S equation of three-dimensional and uncompressible viscous fluid, the aerodynamic characteristics of the train under different schemes are simulated from two directions: the installation position and the height of the air conditioner. The results indicate that when the position of head car’s air conditioner is far from the front while the air conditioner position of the tail car keep close to the end of the train, the total resistance turns out to be the maximum, which is 2.6% higher than other programs; The resistance of the air conditioner decreases as the height of the train increases, when the increase of train height reaches 350 mm, the total resistance of three cars decreased by 7.55%, but at the same time the stability of the train running in the cross-wind will be worse.

    high speed train; numerical calculation; shape optimization; air conditioner; aerodynamic performance

    10.19713/j.cnki.43?1423/u.2018.10.024

    U270.38+3

    A

    1672 ? 7029(2018)10 ? 2631 ? 08

    2017?08?13

    國家自然科學(xué)基金資助項(xiàng)目(U11372360);國家科技支撐計(jì)劃資助項(xiàng)目(2015BAG12B01)

    楊明智(1972?),男,湖南望城人,副教授,從事空氣動力學(xué)研究;E?mail:yqyymz@126.com

    猜你喜歡
    尾車頭車升力
    混勻堆料機(jī)尾車調(diào)偏研究
    高速列車車頂–升力翼組合體氣動特性
    單堆單取斗輪機(jī)尾車防脫鉤保護(hù)裝置及連鎖
    無人機(jī)升力測試裝置設(shè)計(jì)及誤差因素分析
    丁輝:閱兵坦克“頭車”駕駛員
    晚晴(2020年8期)2020-12-03 13:49:57
    頭車半自動鉤緩裝置傾斜問題研究
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    斗輪堆取料機(jī)尾車改造研究
    Numerical simulation of Gurney flap on SFYT15thick airfoil
    升力式再入飛行器體襟翼姿態(tài)控制方法
    亚洲国产毛片av蜜桃av| 啦啦啦在线免费观看视频4| 国产男女内射视频| 久久综合国产亚洲精品| 我要看黄色一级片免费的| 丝袜脚勾引网站| 女人爽到高潮嗷嗷叫在线视频| 超碰成人久久| 亚洲熟女毛片儿| 午夜免费成人在线视频| 国产一区有黄有色的免费视频| 午夜免费观看性视频| 久久久欧美国产精品| 人妻一区二区av| 日韩熟女老妇一区二区性免费视频| 男女之事视频高清在线观看| 免费在线观看黄色视频的| 老司机午夜十八禁免费视频| 国产片内射在线| 91精品国产国语对白视频| 美女国产高潮福利片在线看| 久久精品国产综合久久久| 色老头精品视频在线观看| 国产日韩欧美视频二区| 巨乳人妻的诱惑在线观看| 亚洲成人国产一区在线观看| 精品福利永久在线观看| 成人国产av品久久久| 少妇人妻久久综合中文| 9色porny在线观看| 黄色视频在线播放观看不卡| 蜜桃国产av成人99| 国产伦理片在线播放av一区| 久久久久久久久久久久大奶| 美女主播在线视频| 人人妻人人添人人爽欧美一区卜| 欧美精品亚洲一区二区| 欧美精品高潮呻吟av久久| 久久久久精品人妻al黑| 在线看a的网站| 18禁裸乳无遮挡动漫免费视频| www.熟女人妻精品国产| 国产人伦9x9x在线观看| 王馨瑶露胸无遮挡在线观看| 欧美日韩亚洲高清精品| av线在线观看网站| 免费少妇av软件| 黄色视频不卡| 亚洲激情五月婷婷啪啪| 欧美另类亚洲清纯唯美| 久久精品国产综合久久久| 成人国产一区最新在线观看| 亚洲男人天堂网一区| 国产福利在线免费观看视频| 日本wwww免费看| 国产99久久九九免费精品| 色老头精品视频在线观看| av线在线观看网站| 人人妻人人澡人人爽人人夜夜| 汤姆久久久久久久影院中文字幕| 狂野欧美激情性xxxx| 久久人人爽av亚洲精品天堂| 欧美人与性动交α欧美软件| 午夜免费成人在线视频| 女性生殖器流出的白浆| tocl精华| 国产精品一区二区在线不卡| 老熟妇乱子伦视频在线观看 | 国产91精品成人一区二区三区 | 最近最新中文字幕大全免费视频| 性少妇av在线| 丁香六月欧美| 欧美日韩黄片免| 久久狼人影院| 亚洲精品中文字幕一二三四区 | 国产精品欧美亚洲77777| 国产免费视频播放在线视频| 亚洲精品在线美女| 80岁老熟妇乱子伦牲交| 亚洲av欧美aⅴ国产| av天堂在线播放| 成人亚洲精品一区在线观看| 日本一区二区免费在线视频| 青春草亚洲视频在线观看| 亚洲欧美清纯卡通| 久久久国产欧美日韩av| 欧美日韩亚洲综合一区二区三区_| 午夜福利在线观看吧| 菩萨蛮人人尽说江南好唐韦庄| 国产色视频综合| 国产成人系列免费观看| 老熟妇仑乱视频hdxx| 丝袜脚勾引网站| 在线观看免费高清a一片| 日日爽夜夜爽网站| 国产免费视频播放在线视频| 日本av手机在线免费观看| 国产区一区二久久| 亚洲精品乱久久久久久| 女人精品久久久久毛片| 99久久国产精品久久久| 亚洲伊人久久精品综合| 老司机午夜十八禁免费视频| 国产麻豆69| 午夜激情久久久久久久| 一级毛片电影观看| 久久久久久久精品精品| 久久久精品区二区三区| 国产在线免费精品| 亚洲国产av影院在线观看| 啦啦啦在线免费观看视频4| 久久久久久免费高清国产稀缺| 国产精品久久久人人做人人爽| 狠狠狠狠99中文字幕| 亚洲国产中文字幕在线视频| 午夜激情av网站| 青春草视频在线免费观看| 精品国产一区二区三区久久久樱花| 99精品久久久久人妻精品| 各种免费的搞黄视频| 国产成人免费无遮挡视频| 国产精品一二三区在线看| 欧美变态另类bdsm刘玥| 电影成人av| 欧美一级毛片孕妇| 亚洲自偷自拍图片 自拍| 午夜免费成人在线视频| 久久国产精品大桥未久av| 欧美日韩福利视频一区二区| 亚洲色图综合在线观看| 老司机靠b影院| 欧美黑人欧美精品刺激| 别揉我奶头~嗯~啊~动态视频 | 纯流量卡能插随身wifi吗| 国产91精品成人一区二区三区 | 十八禁人妻一区二区| 黑人欧美特级aaaaaa片| 妹子高潮喷水视频| 精品人妻熟女毛片av久久网站| 亚洲人成电影观看| 国产成人精品久久二区二区免费| 亚洲伊人色综图| 日本黄色日本黄色录像| 精品一区二区三区四区五区乱码| 黄色毛片三级朝国网站| 男人舔女人的私密视频| 俄罗斯特黄特色一大片| 桃花免费在线播放| 欧美 亚洲 国产 日韩一| 精品一区二区三卡| 麻豆av在线久日| 国产精品免费大片| 悠悠久久av| 一二三四社区在线视频社区8| 中文欧美无线码| 亚洲欧洲精品一区二区精品久久久| 韩国高清视频一区二区三区| 日本黄色日本黄色录像| 性少妇av在线| a 毛片基地| 人妻 亚洲 视频| 欧美精品一区二区大全| 国精品久久久久久国模美| 12—13女人毛片做爰片一| 久久久久精品人妻al黑| e午夜精品久久久久久久| 精品卡一卡二卡四卡免费| 亚洲少妇的诱惑av| 国产成人精品久久二区二区91| 色综合欧美亚洲国产小说| 成年人午夜在线观看视频| 巨乳人妻的诱惑在线观看| 91av网站免费观看| 亚洲第一欧美日韩一区二区三区 | 精品国产乱子伦一区二区三区 | 色婷婷av一区二区三区视频| 国产成+人综合+亚洲专区| 国产成人系列免费观看| 久久精品亚洲av国产电影网| 亚洲中文日韩欧美视频| 午夜福利视频在线观看免费| 电影成人av| 亚洲欧洲精品一区二区精品久久久| 欧美成人午夜精品| 欧美黑人欧美精品刺激| 少妇人妻久久综合中文| 成年人午夜在线观看视频| 欧美变态另类bdsm刘玥| 亚洲美女黄色视频免费看| 脱女人内裤的视频| 12—13女人毛片做爰片一| 日本黄色日本黄色录像| 99热国产这里只有精品6| 日韩免费高清中文字幕av| 新久久久久国产一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品中文字幕在线视频| 国产免费现黄频在线看| bbb黄色大片| 又紧又爽又黄一区二区| 久久精品熟女亚洲av麻豆精品| 天天操日日干夜夜撸| 久久香蕉激情| 久久综合国产亚洲精品| 好男人电影高清在线观看| www.av在线官网国产| 久久久久久久久久久久大奶| 天堂中文最新版在线下载| 国产在线观看jvid| 丝袜在线中文字幕| 欧美97在线视频| 国产在线观看jvid| 国产欧美日韩一区二区三 | 大型av网站在线播放| 十八禁人妻一区二区| 国产成人影院久久av| 国产成人免费观看mmmm| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美日韩高清在线视频 | 亚洲情色 制服丝袜| 啦啦啦中文免费视频观看日本| 激情视频va一区二区三区| 91字幕亚洲| 欧美午夜高清在线| 国产精品av久久久久免费| 免费看十八禁软件| 欧美中文综合在线视频| 久久99热这里只频精品6学生| 久久性视频一级片| 在线观看免费午夜福利视频| 青春草亚洲视频在线观看| 国产精品久久久久成人av| 久久久久久人人人人人| 亚洲国产精品999| 超色免费av| 久久ye,这里只有精品| 五月天丁香电影| 777久久人妻少妇嫩草av网站| 超碰97精品在线观看| 少妇精品久久久久久久| 成在线人永久免费视频| 亚洲熟女精品中文字幕| 久久久精品国产亚洲av高清涩受| 熟女少妇亚洲综合色aaa.| 9色porny在线观看| 岛国在线观看网站| 丝袜喷水一区| 少妇裸体淫交视频免费看高清 | 亚洲男人天堂网一区| 一级,二级,三级黄色视频| 电影成人av| 女警被强在线播放| 成人影院久久| 久热爱精品视频在线9| 成年av动漫网址| 国产精品久久久av美女十八| 成人亚洲精品一区在线观看| 国产精品偷伦视频观看了| av网站在线播放免费| 亚洲情色 制服丝袜| 美女福利国产在线| 九色亚洲精品在线播放| 亚洲成人免费av在线播放| 岛国毛片在线播放| 99国产精品一区二区蜜桃av | 国产极品粉嫩免费观看在线| 我的亚洲天堂| 色婷婷久久久亚洲欧美| 午夜精品国产一区二区电影| 国产av一区二区精品久久| 肉色欧美久久久久久久蜜桃| 久9热在线精品视频| 爱豆传媒免费全集在线观看| 国产成人a∨麻豆精品| 99热全是精品| 十八禁高潮呻吟视频| 国产又色又爽无遮挡免| 亚洲欧美一区二区三区黑人| 亚洲第一欧美日韩一区二区三区 | 国产熟女午夜一区二区三区| 亚洲精品国产av成人精品| 午夜免费观看性视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美激情久久久久久爽电影 | 国产精品av久久久久免费| 婷婷丁香在线五月| 欧美人与性动交α欧美软件| 日本五十路高清| 12—13女人毛片做爰片一| 动漫黄色视频在线观看| 日韩欧美一区二区三区在线观看 | 久久久久久久国产电影| 亚洲成人国产一区在线观看| 免费女性裸体啪啪无遮挡网站| 成年女人毛片免费观看观看9 | 国产成人啪精品午夜网站| 久久国产亚洲av麻豆专区| 国产成人欧美| 久久久久久久国产电影| 亚洲精品粉嫩美女一区| 欧美日韩一级在线毛片| 99国产精品99久久久久| 国产精品影院久久| 女人久久www免费人成看片| 国产极品粉嫩免费观看在线| 高潮久久久久久久久久久不卡| 精品国产超薄肉色丝袜足j| 亚洲第一青青草原| 免费看十八禁软件| 成人三级做爰电影| 不卡av一区二区三区| 91字幕亚洲| 国产精品国产av在线观看| 亚洲av欧美aⅴ国产| 777米奇影视久久| 国产精品久久久久成人av| 好男人电影高清在线观看| 日本精品一区二区三区蜜桃| 永久免费av网站大全| 99久久精品国产亚洲精品| 精品国产乱子伦一区二区三区 | 亚洲精品美女久久久久99蜜臀| 十八禁人妻一区二区| 男女高潮啪啪啪动态图| 色婷婷久久久亚洲欧美| 涩涩av久久男人的天堂| 成人免费观看视频高清| 香蕉国产在线看| 久久国产精品人妻蜜桃| 在线观看舔阴道视频| 国产在线一区二区三区精| 国产黄色免费在线视频| 老熟妇乱子伦视频在线观看 | 亚洲精品乱久久久久久| 成年人午夜在线观看视频| av国产精品久久久久影院| 午夜福利视频精品| 国产精品久久久久久精品电影小说| 无遮挡黄片免费观看| av有码第一页| 最黄视频免费看| av不卡在线播放| 国产伦理片在线播放av一区| 国产一区二区三区av在线| 高清在线国产一区| 欧美精品亚洲一区二区| 免费在线观看完整版高清| 国产精品久久久久久人妻精品电影 | 淫妇啪啪啪对白视频 | 久久人人97超碰香蕉20202| 视频区欧美日本亚洲| 成人黄色视频免费在线看| 咕卡用的链子| 日韩电影二区| 男女下面插进去视频免费观看| 黑人猛操日本美女一级片| 建设人人有责人人尽责人人享有的| 国产伦理片在线播放av一区| 亚洲人成电影免费在线| 亚洲精品久久午夜乱码| 97在线人人人人妻| av网站免费在线观看视频| 国产高清videossex| 欧美精品啪啪一区二区三区 | 飞空精品影院首页| 久久女婷五月综合色啪小说| 老司机午夜十八禁免费视频| 热99国产精品久久久久久7| 在线永久观看黄色视频| 成人亚洲精品一区在线观看| 欧美一级毛片孕妇| 亚洲国产欧美日韩在线播放| 成在线人永久免费视频| 国产精品一区二区在线不卡| 巨乳人妻的诱惑在线观看| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 日韩视频一区二区在线观看| 国产成人免费观看mmmm| 久久久久久久国产电影| 又大又爽又粗| 91九色精品人成在线观看| 久久久久网色| 日韩视频一区二区在线观看| 精品久久蜜臀av无| 91大片在线观看| 欧美 日韩 精品 国产| 欧美激情高清一区二区三区| 老司机午夜福利在线观看视频 | 巨乳人妻的诱惑在线观看| 久久99一区二区三区| 久久99热这里只频精品6学生| 亚洲少妇的诱惑av| 成人亚洲精品一区在线观看| 最新的欧美精品一区二区| 久久久久久久精品精品| 美女福利国产在线| 国产在线一区二区三区精| 国产伦理片在线播放av一区| 女人精品久久久久毛片| 国产老妇伦熟女老妇高清| 天天影视国产精品| 久久久精品免费免费高清| 视频在线观看一区二区三区| 亚洲欧美成人综合另类久久久| √禁漫天堂资源中文www| 亚洲三区欧美一区| 婷婷成人精品国产| 91精品三级在线观看| 日韩欧美免费精品| 18禁裸乳无遮挡动漫免费视频| 国产人伦9x9x在线观看| 一本综合久久免费| 老汉色av国产亚洲站长工具| 黑人巨大精品欧美一区二区蜜桃| 久久久精品区二区三区| 欧美 日韩 精品 国产| 狠狠婷婷综合久久久久久88av| 法律面前人人平等表现在哪些方面 | 久久精品国产a三级三级三级| 亚洲久久久国产精品| 日韩人妻精品一区2区三区| 黑人猛操日本美女一级片| 天天躁日日躁夜夜躁夜夜| 久久久精品免费免费高清| 国产男人的电影天堂91| 人人妻人人添人人爽欧美一区卜| avwww免费| 各种免费的搞黄视频| 王馨瑶露胸无遮挡在线观看| 男女床上黄色一级片免费看| 国产精品.久久久| 日韩精品免费视频一区二区三区| 国产一区二区三区av在线| 欧美精品亚洲一区二区| 免费观看人在逋| 国产一区二区 视频在线| 在线观看一区二区三区激情| 老鸭窝网址在线观看| 精品第一国产精品| 两个人免费观看高清视频| 亚洲美女黄色视频免费看| 欧美性长视频在线观看| 一二三四社区在线视频社区8| 夫妻午夜视频| 亚洲成人免费电影在线观看| 亚洲精品中文字幕一二三四区 | 欧美97在线视频| 18在线观看网站| 国产男女内射视频| 美女视频免费永久观看网站| 精品国产一区二区三区久久久樱花| 国产成人啪精品午夜网站| 侵犯人妻中文字幕一二三四区| 手机成人av网站| 天天躁狠狠躁夜夜躁狠狠躁| 18在线观看网站| 中文字幕高清在线视频| 久久久久久久大尺度免费视频| 久久久国产精品麻豆| 欧美一级毛片孕妇| 热re99久久精品国产66热6| 免费在线观看完整版高清| 各种免费的搞黄视频| 欧美日韩黄片免| 亚洲视频免费观看视频| 国产区一区二久久| 一本综合久久免费| 老汉色av国产亚洲站长工具| av网站在线播放免费| 蜜桃国产av成人99| 老汉色∧v一级毛片| 成人av一区二区三区在线看 | 青草久久国产| 久久99热这里只频精品6学生| 亚洲三区欧美一区| 99国产综合亚洲精品| e午夜精品久久久久久久| 欧美日韩av久久| 性高湖久久久久久久久免费观看| 伊人亚洲综合成人网| 真人做人爱边吃奶动态| 久热这里只有精品99| 一级毛片电影观看| 午夜两性在线视频| 免费高清在线观看视频在线观看| 一本综合久久免费| 亚洲国产精品一区二区三区在线| 中文字幕高清在线视频| 中文字幕人妻丝袜制服| 免费在线观看视频国产中文字幕亚洲 | 亚洲av电影在线进入| 黑人操中国人逼视频| 色94色欧美一区二区| 夫妻午夜视频| 久久av网站| 亚洲精品国产色婷婷电影| 无遮挡黄片免费观看| 91麻豆av在线| 日日夜夜操网爽| 国产日韩欧美亚洲二区| 国产精品成人在线| 不卡一级毛片| 一区二区三区乱码不卡18| 一个人免费看片子| 2018国产大陆天天弄谢| 日韩欧美国产一区二区入口| 欧美xxⅹ黑人| 国产老妇伦熟女老妇高清| 天天影视国产精品| 午夜激情久久久久久久| av线在线观看网站| 亚洲欧洲日产国产| 欧美黑人精品巨大| 久久久精品94久久精品| 国产精品影院久久| 免费在线观看黄色视频的| 99国产极品粉嫩在线观看| 精品一区二区三区四区五区乱码| 国产在线一区二区三区精| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美一区二区三区久久| 亚洲av日韩在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 美女视频免费永久观看网站| 日本a在线网址| 在线 av 中文字幕| 国产一区二区三区在线臀色熟女 | 大陆偷拍与自拍| 一本色道久久久久久精品综合| 欧美国产精品va在线观看不卡| 亚洲人成77777在线视频| 午夜福利影视在线免费观看| 视频区欧美日本亚洲| 午夜福利视频精品| 久久精品人人爽人人爽视色| 青青草视频在线视频观看| 日韩视频在线欧美| 18禁裸乳无遮挡动漫免费视频| 久久人人97超碰香蕉20202| 日韩大片免费观看网站| 老熟妇乱子伦视频在线观看 | 国产一区二区激情短视频 | 亚洲中文日韩欧美视频| 香蕉丝袜av| 老司机亚洲免费影院| 18在线观看网站| 99九九在线精品视频| 最近最新免费中文字幕在线| 国产高清国产精品国产三级| 欧美黑人欧美精品刺激| 精品国内亚洲2022精品成人 | 日韩制服骚丝袜av| 精品亚洲成a人片在线观看| 黄色片一级片一级黄色片| 亚洲 国产 在线| 美女高潮到喷水免费观看| 精品少妇内射三级| 99国产精品免费福利视频| 大片电影免费在线观看免费| 国产欧美日韩一区二区三区在线| 国产精品秋霞免费鲁丝片| 国产精品久久久久久精品古装| 久久这里只有精品19| 日本撒尿小便嘘嘘汇集6| 亚洲成人免费电影在线观看| 制服诱惑二区| 国产深夜福利视频在线观看| 女人精品久久久久毛片| 天天躁狠狠躁夜夜躁狠狠躁| www.熟女人妻精品国产| 男人添女人高潮全过程视频| 国产激情久久老熟女| 侵犯人妻中文字幕一二三四区| 成人手机av| 久久久久久久大尺度免费视频| 色婷婷av一区二区三区视频| 午夜福利乱码中文字幕| 黄片播放在线免费| 亚洲精品粉嫩美女一区| 久久国产精品影院| 日韩欧美免费精品| 黄色视频在线播放观看不卡| 男女边摸边吃奶| 亚洲国产欧美日韩在线播放| 超色免费av| 婷婷成人精品国产| 母亲3免费完整高清在线观看| 国产99久久九九免费精品| 精品福利永久在线观看| 亚洲中文日韩欧美视频| 国产成人精品久久二区二区91| 9色porny在线观看| 精品久久久久久久毛片微露脸 | 99精品欧美一区二区三区四区| 9色porny在线观看| 久久国产精品大桥未久av| 99香蕉大伊视频| 亚洲国产看品久久| 大片电影免费在线观看免费| 黄片小视频在线播放| xxxhd国产人妻xxx| 免费观看人在逋| 777久久人妻少妇嫩草av网站| 午夜精品久久久久久毛片777| 性色av一级| 免费不卡黄色视频| 热re99久久精品国产66热6| 久久久国产成人免费| 国产在线观看jvid| 免费久久久久久久精品成人欧美视频| 人人妻人人爽人人添夜夜欢视频| 精品国产一区二区久久| 亚洲精品中文字幕在线视频| 女警被强在线播放|