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

    動(dòng)車組橫風(fēng)環(huán)境下的交會(huì)氣動(dòng)效應(yīng)

    2016-06-24 07:15:15劉峰姚松張潔張娜中南大學(xué)交通運(yùn)輸工程學(xué)院軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室湖南長沙410075
    關(guān)鍵詞:橫風(fēng)交會(huì)動(dòng)車組

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

    ?

    動(dòng)車組橫風(fēng)環(huán)境下的交會(huì)氣動(dòng)效應(yīng)

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

    摘要:采用三維、可壓縮、非定常N?S方程的數(shù)值計(jì)算方法,對8輛編組的動(dòng)車組在20 m/s橫風(fēng)下以250 km/h速度交會(huì)時(shí)列車表面瞬變壓力和車體所受氣動(dòng)力及力矩進(jìn)行分析,并采用間接驗(yàn)證方法,將風(fēng)洞實(shí)驗(yàn)、動(dòng)模型實(shí)驗(yàn)得到的結(jié)果分別與數(shù)值模擬結(jié)果進(jìn)行對比。研究結(jié)果表明:間接驗(yàn)證方法下所得氣動(dòng)效應(yīng)實(shí)驗(yàn)結(jié)果和數(shù)值模擬結(jié)果變化規(guī)律一致,壓力幅值相對誤差在 5%以內(nèi);動(dòng)車組橫風(fēng)下交會(huì)時(shí),車體頭、尾處測點(diǎn)壓力差別較大,中部位于同側(cè)測點(diǎn)壓力差異較小,同一高度、不同縱向測點(diǎn)的壓力變化波形及幅值基本一致,車體頂部測點(diǎn)壓力始終為負(fù);對于車體所受橫向氣動(dòng)力及傾覆力矩,頭車比中間車和尾車的大,背風(fēng)車比迎風(fēng)車的大;隨著橫風(fēng)風(fēng)速的增加,列車所受橫向氣動(dòng)力及傾覆力矩峰值也迅速增加,嚴(yán)重威脅著動(dòng)車組的安全運(yùn)行。

    關(guān)鍵詞:動(dòng)車組;橫風(fēng);交會(huì);風(fēng)洞實(shí)驗(yàn);動(dòng)模型實(shí)驗(yàn);氣動(dòng)效應(yīng)

    動(dòng)車組高速交會(huì)時(shí),交會(huì)側(cè)的空氣壓力發(fā)生突變,產(chǎn)生瞬態(tài)壓力沖擊,會(huì)對車體鋼結(jié)構(gòu)、側(cè)窗和車體橫向穩(wěn)定性帶來不利影響[1?3]。交會(huì)過程中若遭遇較大的橫風(fēng),可使交會(huì)壓力波幅值增加,列車將承受更大的橫向沖擊載荷[4?5]。我國地形及氣象條件比較復(fù)雜,不少行車區(qū)段已處于強(qiáng)風(fēng)區(qū),而隨著高速鐵路的不斷發(fā)展,強(qiáng)風(fēng)區(qū)內(nèi)動(dòng)車組交會(huì)不可避免[6]:因此,研究動(dòng)車組橫風(fēng)下交會(huì)時(shí)的氣動(dòng)效應(yīng)具有很強(qiáng)的現(xiàn)實(shí)意義。對于橫風(fēng)效應(yīng)和列車交會(huì)問題,國內(nèi)外學(xué)者開展了諸多數(shù)值模擬研究[7?14],但將二者進(jìn)行耦合分析即同時(shí)考慮列車橫風(fēng)環(huán)境下發(fā)生交會(huì)的研究較少。為此,本文作者以國內(nèi)處于風(fēng)區(qū)中的某客運(yùn)專線高速鐵路為研究對象,采用數(shù)值模擬方法對動(dòng)車組在橫風(fēng)下以250 km/h 速度交會(huì)時(shí)的車體表面測點(diǎn)壓力、車體所受橫向氣動(dòng)力及傾覆力矩變化進(jìn)行分析,為動(dòng)車組在風(fēng)區(qū)內(nèi)交會(huì)安全性提供氣動(dòng)依據(jù)。另外,對于列車橫風(fēng)下交會(huì)這種復(fù)雜工況,采用的數(shù)值模擬方法很難通過實(shí)驗(yàn)手段進(jìn)行直接驗(yàn)證。據(jù)文獻(xiàn)[15],在對物理假說進(jìn)行驗(yàn)證時(shí),對于受實(shí)驗(yàn)水平的限制暫時(shí)達(dá)不到要求的情況,可以采用間接驗(yàn)證的方法?;谶@一思路,采用風(fēng)洞模型實(shí)驗(yàn)和動(dòng)模型實(shí)驗(yàn),將動(dòng)車組單列橫風(fēng)下運(yùn)行、2 列無風(fēng)下交會(huì)運(yùn)行的動(dòng)車數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果分別進(jìn)行對比,以間接驗(yàn)證本文所采用數(shù)值計(jì)算方法的可行性。

    1數(shù)值計(jì)算模型及計(jì)算方法

    依照所選列車及線路的實(shí)際幾何特征,建立列車橫風(fēng)交會(huì)數(shù)值計(jì)算模型。交會(huì)列車選用國內(nèi)高鐵線路上較常見的2款動(dòng)車組,為了與實(shí)際運(yùn)營盡量一致,列車模型均采用8輛編組,即頭車、6 節(jié)中間車和尾車,2款動(dòng)車組全長分別為201.4 m和205.2 m。根據(jù)“鐵路200~250 km/h 既有線技術(shù)管理辦法”,動(dòng)車組在環(huán)境風(fēng)速不大于20 m/s時(shí)可以正常速度運(yùn)行,本文選擇動(dòng)車組交會(huì)速度為250 km/h,環(huán)境風(fēng)速為20 m/s,風(fēng)向角為90°。交會(huì)開始前2車頭相距100 m,線間距為5m??拷鼨M風(fēng)入口的列車為迎風(fēng)車,離橫風(fēng)入口較遠(yuǎn)的列車為背風(fēng)車。采用混合網(wǎng)格對計(jì)算區(qū)域進(jìn)行離散,即在幾何外形復(fù)雜的列車周圍區(qū)域采用非結(jié)構(gòu)網(wǎng)格,其他區(qū)域采用結(jié)構(gòu)化網(wǎng)格。同時(shí)加密列車周圍網(wǎng)格密度,遠(yuǎn)離車體網(wǎng)格采用稀疏網(wǎng)格,密網(wǎng)格和稀疏網(wǎng)格之間以一定的增長因子均勻過渡。計(jì)算區(qū)域和計(jì)算網(wǎng)格分別如圖1和圖2所示。

    圖1 計(jì)算區(qū)域Fig.1 Computational zone

    圖2計(jì)算網(wǎng)格Fig.2Computational mesh

    基于三維、非定常、可壓縮、黏性流場對動(dòng)車組橫風(fēng)下交會(huì)進(jìn)行流場數(shù)值分析,采用工程上應(yīng)用廣泛的k?ε湍流模型,時(shí)間步長為0.005 s,詳細(xì)流場控制方程及湍流模型見文獻(xiàn)[16]。為了得到穩(wěn)定初始流場,先讓列車靜止,加載橫風(fēng),待風(fēng)場充分發(fā)展后再啟用滑移網(wǎng)格使列車運(yùn)動(dòng)。采用分區(qū)對接滑移網(wǎng)格技術(shù),即在有相對運(yùn)動(dòng)的不同網(wǎng)格間設(shè)置公共滑移界面,從而實(shí)現(xiàn)彼此的數(shù)據(jù)交換。為了對車體表面瞬變壓力進(jìn)行分析,在動(dòng)車組表面布置多個(gè)測點(diǎn)。迎風(fēng)車與背風(fēng)車布點(diǎn)方式一致,單列車共66個(gè);頭車和尾車測點(diǎn)布置位置相同,且數(shù)量相對較多;車體中部壓力變化較小,所以,只在每節(jié)中間車的同一截面布置5個(gè)測點(diǎn)。具體布點(diǎn)方式如圖3所示。

    圖3 動(dòng)車組模型測點(diǎn)布置圖Fig.3 Points arrangement of EMU model

    2 橫風(fēng)下交會(huì)空氣動(dòng)力效應(yīng)分析

    2.1動(dòng)車組表面壓力變化

    當(dāng)動(dòng)車組無風(fēng)交會(huì)、車體表面測點(diǎn)壓力在對方頭車鼻尖到達(dá)時(shí),會(huì)產(chǎn)生1個(gè)正、負(fù)脈沖,即頭波;在最大負(fù)脈沖出現(xiàn)后開始等幅波動(dòng),直到對方尾車鼻尖通過時(shí),則會(huì)產(chǎn)生1個(gè)負(fù)、正脈沖,即尾波。由于本文選用動(dòng)車組編組較長,故尾波幅度遠(yuǎn)小于頭波幅度[17]。

    當(dāng)動(dòng)車組交會(huì)發(fā)生在橫風(fēng)下時(shí),流場則會(huì)同時(shí)具有橫風(fēng)特性和交會(huì)特性。一方面,橫風(fēng)會(huì)改變動(dòng)車組表面尤其是頭尾區(qū)域原有的壓力分布,致使交會(huì)壓力波也發(fā)生相應(yīng)改變。以迎風(fēng)車車頭為例,橫風(fēng)作用會(huì)使車頭表面正壓區(qū)朝逆風(fēng)向偏移,負(fù)壓區(qū)朝順風(fēng)向偏移,因此,當(dāng)背風(fēng)車車身測點(diǎn)與其發(fā)生交會(huì)時(shí),單由交會(huì)所產(chǎn)生的頭波正脈沖被削減,負(fù)脈沖被加強(qiáng)。另一方面,2 列動(dòng)車組在橫風(fēng)下交會(huì)時(shí),由于彼此之間具有擋風(fēng)作用,使的車體表面壓力幅值會(huì)產(chǎn)生進(jìn)一步改變。圖 4所示為動(dòng)車組橫風(fēng)下交會(huì)時(shí)車體表面的壓力分布。

    圖4動(dòng)車組橫風(fēng)下交會(huì)車體表面壓力云圖Fig.4PressureContour of EMU in passing each other underCrosswind

    圖5和圖 6所示分別為迎風(fēng)車、背風(fēng)車車身中部測點(diǎn)(21號測點(diǎn))的壓力變化曲線。從圖5和圖6可以看出:由于橫風(fēng)和交會(huì)壓力波的雙重作用,動(dòng)車組橫風(fēng)交會(huì)時(shí)車體表面測點(diǎn)壓力與無風(fēng)交會(huì)時(shí)有明顯不同;對于迎風(fēng)車車體交會(huì)側(cè)測點(diǎn),與無風(fēng)交會(huì)相比,壓力頭波正波幅值較大,負(fù)波幅值較?。欢鴫毫ξ膊ㄕǚ递^小,負(fù)波幅值較大;頭、尾波之間的區(qū)間段壓力波動(dòng)差別不大;對于背風(fēng)車車體交會(huì)側(cè)測點(diǎn),幅值較大的為頭波負(fù)波和尾波正波;非交會(huì)側(cè)壓力主要受橫風(fēng)影響,有風(fēng)與無風(fēng)時(shí),測點(diǎn)壓力幅值有一定的變化,而交會(huì)頭波和尾波帶來的幅值變化較小。近(4號測點(diǎn)),分別為?1.76 kPa和2.37 kPa。尾車測點(diǎn)壓力幅值相對較小,鼻尖壓力(66 號測點(diǎn))最大值及峰峰值分別為0.30 kPa和0.54 kPa。

    圖5 迎風(fēng)車車身中部測點(diǎn)壓力變化曲線Fig.5 Curves of pressure for measuring points in the middle of upwind vehicle

    圖6 背風(fēng)車車身中部測點(diǎn)壓力變化曲線Fig.6 Curves of pressure for measuring points in the middle of leeward vehicle

    2.2動(dòng)車組壓力變化幅值

    圖7所示為動(dòng)車組頭車對稱面上不同位置測點(diǎn)(見圖3)的壓力變化曲線。從圖7可以看出:在同一時(shí)刻,頭車鼻尖與車身過渡處壓力相差4.00 kPa左右;動(dòng)車組橫風(fēng)下交會(huì)時(shí),除車體交會(huì)側(cè)與非交會(huì)側(cè)壓力幅值差異較大外,頭、尾車由于曲率變化較大,其測點(diǎn)壓力幅值差異也比較大。以迎風(fēng)車車體表面測點(diǎn)壓力為例,交會(huì)過程中車體表面最大正壓發(fā)生在車頭鼻尖附近(1號測點(diǎn)),為3.55 kPa,最大負(fù)壓及最大壓力峰峰值(最大值?最小值)發(fā)生在車頭鼻尖偏交會(huì)側(cè)附

    動(dòng)車組中部幾何結(jié)構(gòu)相對簡單,位于車體同側(cè)的測點(diǎn)壓力變化較小。圖 8所示為車體中部同一斷面不同測點(diǎn)壓力變化曲線。從圖8可以看出:對于動(dòng)車組同側(cè)的測點(diǎn),隨著車體垂向高度不同,測點(diǎn)壓力變化不大,壓力幅值及壓力峰峰值僅相差0.05 kPa左右;車體頂部測點(diǎn)壓力始終為負(fù)值,其最大正壓及峰峰值分別為?0.16 kPa和0.41kPa。沿車身不同縱向位置的測點(diǎn),由于發(fā)生交會(huì)時(shí)刻不同,使得其壓力變化也有所不同。圖9所示為車體交會(huì)側(cè)不同縱向位置測點(diǎn)壓力變化曲線。從圖9可以看出:除了相位有所不同外,各測點(diǎn)壓力變化波形及幅值基本一致。背風(fēng)車的各測點(diǎn)壓力變化波形及幅值與迎風(fēng)車的類似。

    圖7 車頭不同位置測點(diǎn)壓力變化曲線Fig.7 Curves of pressure for measuring points of different locations on head surface

    圖8 車體中部同一截面不同測點(diǎn)壓力變化曲線Fig.8 Curves of pressure for different measuring points in the same section of the middle of EMU

    圖9 車體中部不同縱向位置測點(diǎn)壓力變化曲線Fig.9 Curves of pressure for measuring points of different longitudinal locations on the middle part of train

    圖10 迎風(fēng)車所受橫向氣動(dòng)力變化曲線Fig.10 Curves of lateral force acting on upwind vehicle

    2.3動(dòng)車組所受氣動(dòng)力及力矩

    圖10~13所示分別為動(dòng)車組橫風(fēng)下交會(huì)時(shí),車體所受橫向氣動(dòng)力及傾覆力矩變化曲線。由于中間車所受力及力矩變化幅值相差較小,因此,只選擇第4節(jié)車進(jìn)行分析。從圖10~13可以看出:受交會(huì)壓力波的排斥、吸引作用,迎風(fēng)車所受橫向氣動(dòng)力及傾覆力矩在交會(huì)開始時(shí)先減小后增大,而在交會(huì)結(jié)束時(shí)先增大后減??;而背風(fēng)車在交會(huì)開始時(shí)先增大后減小,在交會(huì)結(jié)束時(shí)先減小后增大;與中間車和尾車相比,頭車所受橫向氣動(dòng)力及傾覆力矩較大;當(dāng)頭車在另一列車區(qū)段內(nèi)行駛時(shí),由于對方動(dòng)車組的擋風(fēng)作用,其所受橫向氣動(dòng)力及傾覆力矩較非交會(huì)時(shí)有所減小。

    綜合比較各節(jié)車輛所受氣動(dòng)載荷,背風(fēng)車頭車所受橫向氣動(dòng)力及傾覆力矩相對最大,這是由于交會(huì)初始時(shí)產(chǎn)生的壓力頭波正脈沖與橫風(fēng)作用產(chǎn)生了疊加。表1所示為不同橫風(fēng)風(fēng)速下,動(dòng)車組以250 km/h交會(huì),背風(fēng)車頭車所受橫向氣動(dòng)力及傾覆力矩峰值。從表1可以看出:隨著橫風(fēng)風(fēng)速增加,車體所受橫向氣動(dòng)力及傾覆力矩峰值也迅速增加,當(dāng)風(fēng)速為40 m/s時(shí),背風(fēng)車頭車所受橫向氣動(dòng)力正峰值達(dá)到213.1kN,傾覆力矩正峰值達(dá)到356.2 kN?m,嚴(yán)重影響了動(dòng)車組的安全運(yùn)行穩(wěn)定性。

    圖11 背風(fēng)車所受橫向氣動(dòng)力變化曲線Fig.11 Curves of lateral force acting on leeward vehicle

    圖12迎風(fēng)車所受傾覆力矩變化曲線Fig.12Curves of overturning moment acting on upwind vehicle

    圖13 背風(fēng)車所受傾覆力矩變化曲線Fig.13 Curves of lateral force overturning moment acting on leeward vehicle

    表1 不同橫風(fēng)風(fēng)速下背風(fēng)車頭車所受橫向氣動(dòng)力及傾覆力矩峰值Table1 Lateral force and overturning moment acting on headCar of leeward vehicle underCrosswind with different speeds

    3 模型實(shí)驗(yàn)間接驗(yàn)證

    3.1風(fēng)洞實(shí)驗(yàn)

    風(fēng)洞實(shí)驗(yàn)在長×寬×高為15 m×8 m×6 m的大型低速風(fēng)洞進(jìn)行[18]。選擇3輛編組(頭車+中車+尾車)的迎風(fēng)車動(dòng)車組為研究對象。模型縮比為1:15,合成風(fēng)速為60 m/s,側(cè)滑角分別為0°,30°,60°,70°,80° 和 90°。正式實(shí)驗(yàn)前,通過調(diào)整風(fēng)速,確保實(shí)驗(yàn)雷諾數(shù)滿足自模擬區(qū)域,即氣動(dòng)力系數(shù)、力矩系數(shù)不再隨來流速度的變化而改變。對于單車的數(shù)值模擬部分,模型縮比為1:15,合成速度為60 m/s,馬赫數(shù)Ma小于0.3,可認(rèn)為是不可壓縮流動(dòng)問題,但為了保證與模擬雙車交會(huì)時(shí)計(jì)算方法的統(tǒng)一性,按照可壓縮問題進(jìn)行處理。圖14所示為在相同來流條件下,相同外形的動(dòng)車組在平地上所受橫向氣動(dòng)力及傾覆力矩的數(shù)值計(jì)算結(jié)果與風(fēng)洞實(shí)驗(yàn)結(jié)果對比。從圖14可以看出:數(shù)值計(jì)算所得橫向氣動(dòng)力和傾覆力矩與實(shí)驗(yàn)結(jié)果變化規(guī)律一致,除 90°側(cè)滑角外,其他工況計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果相對誤差在 5%以內(nèi)。計(jì)算橫風(fēng)交會(huì)工況的側(cè)滑角并沒有涉及 90°等大角度,這說明本文模擬橫風(fēng)所使用的計(jì)算方法是合理的。

    圖14 風(fēng)洞實(shí)驗(yàn)結(jié)果與計(jì)算結(jié)果比較Fig.14Comparison between wind tunnel test results andCalculated results

    3.2動(dòng)模型實(shí)驗(yàn)

    采用中南大學(xué)軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室的列車氣動(dòng)性能動(dòng)模型實(shí)驗(yàn)裝置,驗(yàn)證本文對于動(dòng)車組無風(fēng)交會(huì)時(shí)氣動(dòng)效應(yīng)模擬的有效性。該實(shí)驗(yàn)裝置可以真實(shí)再現(xiàn)動(dòng)車組高速交會(huì)時(shí)的空氣三維非定常流動(dòng)現(xiàn)象,是模擬列車交會(huì)氣動(dòng)問題的專用大模型實(shí)驗(yàn)設(shè)備[19]。實(shí)驗(yàn)所用2列動(dòng)車組模型與數(shù)值計(jì)算模型完全一致,編組同為8 輛,模型縮比均為1:31。根據(jù)流動(dòng)相似原理,動(dòng)模型實(shí)驗(yàn)中模擬雷諾數(shù)是重要的相似準(zhǔn)則,而當(dāng)實(shí)驗(yàn)雷諾數(shù)大于臨界雷諾數(shù)3.6×10 5(由實(shí)驗(yàn)測得)時(shí),氣動(dòng)特性就不再隨雷諾數(shù)變化而改變[19]。本次動(dòng)模型實(shí)驗(yàn)列車速度為250 km/h,特征長度取模型車高0.119 m,在大氣壓為1×105Pa下,20℃時(shí)的空氣運(yùn)動(dòng)黏性系數(shù)為15.0×10?6m2/s,得出實(shí)驗(yàn)雷諾數(shù)為 5.5×105,大于臨界雷諾數(shù),說明本次的實(shí)驗(yàn)結(jié)果滿足相似性要求。

    兩動(dòng)車組在無風(fēng)下以250 km/h 速度交會(huì)時(shí),數(shù)值計(jì)算與動(dòng)模型實(shí)驗(yàn)中位于相同位置的測點(diǎn)壓力變化進(jìn)行對比。圖15所示為動(dòng)車組中部交會(huì)側(cè)測點(diǎn)(21號測點(diǎn))隨時(shí)間的壓力變化曲線對比。從圖15可以看出:通過數(shù)值計(jì)算與動(dòng)模型實(shí)驗(yàn)得到的車體表面壓力曲線變化規(guī)律基本一致,幅值相對誤差在 5%以內(nèi),說明本文采用的數(shù)值計(jì)算方法可以較好地模擬動(dòng)車組交會(huì)時(shí)所引發(fā)的空氣動(dòng)力效應(yīng)問題。需注意的是:在用模型實(shí)驗(yàn)間接驗(yàn)證數(shù)值計(jì)算方法的可行性時(shí),只針對模型縮比影響進(jìn)行分析,對于列車模型編組差異帶來的影響將有待進(jìn)一步研究。

    圖15 動(dòng)模型實(shí)驗(yàn)結(jié)果與計(jì)算結(jié)果比較Fig.15 Comparison between moving model experiment results andCalculation results

    4 結(jié)論

    1)風(fēng)洞實(shí)驗(yàn)、 動(dòng)模型實(shí)驗(yàn)得到的結(jié)果分別與數(shù)值模擬進(jìn)行對比,所得到氣動(dòng)效應(yīng)變化規(guī)律基本一致,動(dòng)車組表面測點(diǎn)壓力幅值相對誤差在 5%以內(nèi),從而間接證明了本文采用的數(shù)值計(jì)算方法能夠較好地模擬動(dòng)車組橫風(fēng)下交會(huì)時(shí)所誘發(fā)的空氣動(dòng)力效應(yīng)問題。

    2)動(dòng)車組在20 m/s橫風(fēng)下以250 km/h交會(huì)時(shí),車體表面最大正壓發(fā)生在車頭鼻尖附近,最大負(fù)壓及最大壓力峰峰值發(fā)生在車頭鼻尖偏交會(huì)側(cè)附近,且在車體頭、尾部不同位置測點(diǎn)壓力差別較大,頭車鼻尖處壓力比車身過渡處大4.00 kPa左右;車體頂部測點(diǎn)壓力始終為負(fù)值;車體中部處于同側(cè)的不同測點(diǎn)壓力差別較小,不同高度側(cè)點(diǎn)壓力幅值相差0.05 kPa左右,同一高度不同縱向測點(diǎn)壓力變化波形及幅值基本一致。

    3)迎風(fēng)車所受橫向氣動(dòng)力及傾覆力矩在交會(huì)開始時(shí)先減小后增大,而在交會(huì)結(jié)束時(shí)先增大后減小,背風(fēng)車與之相反;與迎風(fēng)車相比,背風(fēng)車所受橫向氣動(dòng)力及傾覆力矩較大;與中間車和尾車相比,頭車所受橫向氣動(dòng)力及傾覆力矩較大;隨著橫風(fēng)風(fēng)速增大,車體所受橫向氣動(dòng)力及傾覆力矩峰值也迅速增大,嚴(yán)重威脅著動(dòng)車組的安全運(yùn)行。

    參考文獻(xiàn):

    [1]李人憲,趙晶,劉杰,等.高速列車會(huì)車壓力波對側(cè)窗的影響[J].機(jī)械工程學(xué)報(bào),2010,46(4): 87?92.LI Renxian,ZHAO Jing,LIU Jie,et al.Influence of air pressure pulse on side windows of high-speed trains passing each other[J].Journal of Mechanical Engineering,2010,46(4): 87?92.

    [2]熊小慧,梁習(xí)鋒.CRH2型動(dòng)車組列車交會(huì)空氣壓力波試驗(yàn)分析[J].鐵道學(xué)報(bào),2009,31(6):15?20.XIONG Xiaohui,LIANG Xifeng.Analysis of air pressure pulses in meeting ofCRH2EMU Trains[J].Journal of theChina Railway Society,2009,31(6):15?20.

    [3]田紅旗,姚松,姚曙光.列車交會(huì)壓力波對車體和側(cè)窗的影響[J].中國鐵道科學(xué),2000,21(4): 6?12.TIAN Hongqi,YAO Song,YAO Shuguang.Influence of the air pressure pulse onCar-body and side-windows of two meeting trains[J].China Railway Science,2000,21(4): 6?12.

    [4]陶澤平,楊志剛,陳羽.側(cè)風(fēng)風(fēng)場特征對高速列車交會(huì)的影響研究[J].力學(xué)與實(shí)踐,2013,35(2): 22?28.TAO Zeping,YANG Zhigang,CHEN Yu.The influence ofCrosswind fieldCharacteristics on operating performance of high speed trains in passing each other[J].Mechanics in Engineering,2013,35(2): 22?28.

    [5]梁習(xí)鋒,沈嫻雅.環(huán)境風(fēng)與列車交會(huì)耦合作用下磁浮列車橫向氣動(dòng)性能[J].中南大學(xué)學(xué)報(bào)學(xué)報(bào)(自然科學(xué)版),2007,38(4): 751?757.LIANG Xifeng,SHEN Xianya.Lateral aerodynamic performances of maglev train when two trains meet with wind blowing[J].Journal ofCentral South University(Science and Technology),2007,38(4): 751?757.

    [6]崔濤,張衛(wèi)華,王琰.側(cè)風(fēng)環(huán)境下列車高速交會(huì)流固耦合振動(dòng)安全性分析[J].振動(dòng)與沖擊,2013,32(11): 75?79.CUI Tao,ZHANG Weihua,WANG Yan.Fluid-solidCoupled vibration safety analysis for two trains passing by each other at high-speed in side wind[J].Journal of Vibration and Shock,2013,32(11): 75?79.

    [7]FUJII T,MAEDA T,ISHIDA H.Wind-induced accidents of train/vehicles and their measure in Japan[J].Quarterly Report of Railway Technical Research Institute,1999,40(1): 50?55.

    [8]肖京平,黃志祥,陳立.高速列車空氣動(dòng)力學(xué)研究技術(shù)綜述[J].力學(xué)與實(shí)踐,2013,35(2):1?10.XIAO Jingping,HUANG Zhixiang,CHEN Li.Review of aerodynamic investigations for high speed train[J].Mechanics in Engineering,2013,35(2):1?10.

    [9]李雪冰,侯傳倫,張曙光,等.高速列車交會(huì)時(shí)的風(fēng)致振動(dòng)研究[J].振動(dòng)與沖擊,2009,28(7): 81?84.LI Xuebing,HOUChuanlun,ZHANG Shuguang,et al.Flow-induced vibration of high-speed train in passing events[J].Journal of Vibration and Shock,2009,28(7): 81?84.

    [10]BAKERC J,JONES J,LOPEZ-CALLEJA F,et al.Measurements of theCross wind forces on trains[J].Journal of Wind Engineering and Industrial Aerodynamics,2004,92(7): 547?563.

    [11]ANDERSSON E,HAGGSTROM J,SIMA M.Assessment of train overturning risk due to strongCross-winds[J].Rail and Rapid Transit,2004,218(3): 213?223.

    [12]毛軍,郗艷紅,楊國偉.側(cè)風(fēng)風(fēng)場特征對高速列車氣動(dòng)性能作用的研究[J].鐵道學(xué)報(bào),2011,33(4): 22?30.MAO Jun,XI Yanhong,YANG Guowei.Research on influence ofCharacteristics ofCross wind field on aerodynamic performance of a high-speed train[J].Journal of theChina Railway Society,2011,33(4): 22?30.

    [13]劉加利,于夢閣,張繼業(yè),等.基于大渦模擬的高速列車橫風(fēng)運(yùn)行安全性研究[J].鐵道學(xué)報(bào),2011,33(4):14?20.LIU Jiali,YU Mengge,ZHANG Jiye,et al.Study on running safety of high-speed train underCrosswind by large eddy simulation[J].Journal of theChina Railway Society,2011,33(4):14?20.

    [14]高廣軍,苗秀娟.強(qiáng)橫風(fēng)下青藏線客車在不同高度橋梁上的氣動(dòng)性能分析[J].中南大學(xué)學(xué)報(bào)學(xué)報(bào)(自然科學(xué)版),2010,41(1): 376?380.GAO Guangjun,MIAO Xiujuan.Aerodynamic performance of passenger train on different heights of bridge of Qinghai—Tibet railway line under strongCross wind[J].Journal ofCentral South University(Science and Technology),2010,41(1): 376?380.

    [15]李忠.只能間接驗(yàn)證的假說的物理特征[J].物理通報(bào),1998,28(8):10?12.LI Zhong.PhysicalCharacteristics of hypothesis validated indirectly only[J].Physics Bulletin,1998,28(8):10?12.

    [16]王福軍.計(jì)算流體動(dòng)力學(xué)分析[M].北京: 清華大學(xué)出版社,2004: 210?215.WANG Fujun.Analysis ofComputational fluid dynamics[M].Beijing: Tsinghua University Press,2004: 210?215.

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

    [18]梁習(xí)鋒,楊明智,周丹,等.蘭新第二雙線路基防風(fēng)風(fēng)洞實(shí)驗(yàn)研究報(bào)告[R].長沙: 中南大學(xué)軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,2010:11?14.LIANG Xifeng,YANG Mingzhi,ZHOU Dan,et al.Research report of wind tunnel experiment of the windbreak of the second Lanzhou—Xinjiang Double Line Railway[R].Changsha:Central South University.Key Laboratory of Traffic Safety on Track of Ministry of Education,2010:11?14.

    [19]田紅旗.中國列車空氣動(dòng)力學(xué)研究進(jìn)展[J].交通運(yùn)輸工程學(xué)報(bào),2006,6(1):1?9.TIAN Hongqi.Study evolvement of train aerodynamics inChina[J].Journal of Traffic and Transportation Engineering,2006,6(1):1?9.

    (編輯 陳燦華)

    Aerodynamic effect of EMU passing by each other underCrosswind

    LIU Feng,YAO Song,ZHANG Jie,ZHANG Na
    (Key Laboratory of Traffic Safety on Track of Ministry of Education,School of Traffic & Transportation Engineering,Central South University,Changsha 410075,China)

    Abstract:Based on unsteady N?S equation of three-dimensional andCompressible viscous fluid,the transient pressure on theCar surface,aerodynamic pressure and its moment were analyzed on theCondition that the EMU(electric multiple units)Consisting of8carriages interact with each other at the speed of 250 km/h under aCrosswind speed of 20 m/s.The results show that byContrasting the results of wind tunnel test and moving model experiment with numerical simulation using indirect verification method,theChange law of aerodynamic effect shows agreeable accordance with an error under 5%.When passing by each other underCrosswind,the pressure difference of measuring points on body and tail is obvious,the measuring points on the same side of middle part shows little difference,the waveform and amplitude of pressure variation on the measure points of different longitude directions with one height are the same generally,and the measurement point on the top of body presents negative pressure.The lateral force and overturning moment of the headCar are greater than those of the middle and the tail,and the later force and overturning moment of the leeward vehicle are greater than those of the upward one.With the increase ofCrosswind speed,the aerodynamic lateral force and overturning moment increase by the power law,which threatens the safe running of EMU greatly.

    Key words:EMU(electric multiple units)?Crosswind? train-crossing? wind tunnel test? moving model experiment? aerodynamic effect

    中圖分類號:U271.91;U25

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1672?7207(2016)01?0307?07

    DOI:10.11817/j.issn.1672-7207.2016.01.042

    收稿日期:2015?01?10;修回日期:2015?03?25

    基金項(xiàng)目(Foundation item):國家自然科學(xué)基金資助項(xiàng)目(11372360,U1134203,51205418);中南大學(xué)自由探索基金資助項(xiàng)目(2014zzts038);中國鐵路總公司科技研究開發(fā)計(jì)劃項(xiàng)目(2014T001-A);湖南省杰出青年基金資助項(xiàng)目(14JJ103)(Projects(11372360,U1134203,51205418)supported by the National Natural Science Foundation ofChina? Project(2014zzts038)supported by the Exploration and Innovation Funds for Graduate Students ofCentral South University? Project(2014T001-A)supported by Science and Technology Research and Development Program ofChina RailwayCorporation? Project(14JJ103)supported by Outstanding Youth Fund of Hunan Province)

    通信作者:姚松,副教授,從事軌道交通安全研究;E-mail: song_yao@csu.edu.cn

    猜你喜歡
    橫風(fēng)交會(huì)動(dòng)車組
    “織交會(huì)”打造大朗樣板
    “95后”動(dòng)車組女司機(jī)的首個(gè)春運(yùn)
    海峽姐妹(2020年2期)2020-03-03 13:36:34
    橫風(fēng)下高速列車駛?cè)胨淼罆r(shí)瞬態(tài)氣動(dòng)性能研究
    2019年,水交會(huì)來了!
    動(dòng)車組BTM帶內(nèi)干擾的排查與整治
    橫風(fēng)和軌道不平順聯(lián)合作用下的車輛—軌道系統(tǒng)隨機(jī)分析模型
    CRH3型動(dòng)車組輪對壓裝曲線研究
    擋風(fēng)墻后動(dòng)車組與普速列車交會(huì)氣動(dòng)特性
    立方體星交會(huì)對接和空間飛行演示
    太空探索(2016年9期)2016-07-12 09:59:53
    高速鐵路動(dòng)車組站內(nèi)對標(biāo)停車難的研究
    人体艺术视频欧美日本| 最新的欧美精品一区二区| 精品一区二区三区四区五区乱码 | 中国三级夫妇交换| 国产成人a∨麻豆精品| 久久久久久免费高清国产稀缺| 另类亚洲欧美激情| 免费女性裸体啪啪无遮挡网站| 欧美精品一区二区免费开放| 搡老乐熟女国产| 亚洲国产最新在线播放| 新久久久久国产一级毛片| 九草在线视频观看| 国产精品成人在线| 激情五月婷婷亚洲| 国产亚洲一区二区精品| 蜜桃在线观看..| 电影成人av| 色综合欧美亚洲国产小说| 在线观看免费高清a一片| 婷婷色av中文字幕| 日韩 欧美 亚洲 中文字幕| 色播在线永久视频| 18禁动态无遮挡网站| 亚洲男人天堂网一区| 老鸭窝网址在线观看| 日韩制服骚丝袜av| 亚洲欧美清纯卡通| 青春草视频在线免费观看| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 777久久人妻少妇嫩草av网站| 精品一区二区三卡| 久久精品人人爽人人爽视色| 一级毛片 在线播放| 亚洲成人手机| 国产亚洲精品第一综合不卡| 制服丝袜香蕉在线| 免费av中文字幕在线| 亚洲欧美中文字幕日韩二区| 亚洲视频免费观看视频| 亚洲av中文av极速乱| 美女扒开内裤让男人捅视频| 97人妻天天添夜夜摸| 一区福利在线观看| 永久免费av网站大全| 成人国产麻豆网| 中文字幕制服av| 久久久久国产精品人妻一区二区| 夫妻午夜视频| 各种免费的搞黄视频| 久久久精品国产亚洲av高清涩受| 亚洲欧洲国产日韩| 一本色道久久久久久精品综合| 久久久国产一区二区| 久久ye,这里只有精品| xxxhd国产人妻xxx| 亚洲欧美一区二区三区久久| 女性被躁到高潮视频| 一区二区三区激情视频| 韩国高清视频一区二区三区| 纯流量卡能插随身wifi吗| 天天影视国产精品| 国产在线视频一区二区| 美女大奶头黄色视频| 成人黄色视频免费在线看| 在线观看三级黄色| 丝袜美足系列| 黄网站色视频无遮挡免费观看| 亚洲精品国产av蜜桃| 91精品三级在线观看| 亚洲熟女精品中文字幕| 日本av手机在线免费观看| 如何舔出高潮| 国产熟女欧美一区二区| 亚洲欧洲日产国产| 久久精品久久久久久久性| 尾随美女入室| 精品国产一区二区久久| 在现免费观看毛片| 国产午夜精品一二区理论片| 观看av在线不卡| 一区二区三区四区激情视频| 日本91视频免费播放| 国产有黄有色有爽视频| 国产成人啪精品午夜网站| 国产精品免费视频内射| 国产成人一区二区在线| 国产 一区精品| 免费观看性生交大片5| 国产精品蜜桃在线观看| 欧美日韩亚洲综合一区二区三区_| 中国三级夫妇交换| 国产亚洲av高清不卡| 亚洲精华国产精华液的使用体验| 亚洲色图 男人天堂 中文字幕| xxxhd国产人妻xxx| 亚洲伊人色综图| av.在线天堂| 日韩av免费高清视频| 亚洲国产精品成人久久小说| 亚洲,欧美,日韩| 两性夫妻黄色片| 久久99热这里只频精品6学生| 午夜精品国产一区二区电影| 黄频高清免费视频| 男女高潮啪啪啪动态图| 18禁观看日本| 最近最新中文字幕大全免费视频 | av在线播放精品| 亚洲专区中文字幕在线 | 久久久久精品性色| 亚洲在久久综合| 亚洲av成人精品一二三区| 日韩欧美精品免费久久| 蜜桃在线观看..| 男人操女人黄网站| 亚洲欧美清纯卡通| 国产成人一区二区在线| 国产麻豆69| 无限看片的www在线观看| 日韩电影二区| 欧美日韩一级在线毛片| 成年人午夜在线观看视频| 国产精品久久久人人做人人爽| 超色免费av| 日韩成人av中文字幕在线观看| 中文字幕制服av| 肉色欧美久久久久久久蜜桃| 久久热在线av| 久久久久久久大尺度免费视频| 国产精品一区二区在线观看99| 中国三级夫妇交换| 亚洲精品国产区一区二| av.在线天堂| 国产精品 欧美亚洲| 青春草国产在线视频| 捣出白浆h1v1| 中文精品一卡2卡3卡4更新| 日韩 亚洲 欧美在线| 国产毛片在线视频| 日日撸夜夜添| 国产亚洲一区二区精品| 亚洲欧美精品自产自拍| 久久人人爽av亚洲精品天堂| 在线观看一区二区三区激情| av在线观看视频网站免费| 中文字幕av电影在线播放| 丁香六月天网| 香蕉国产在线看| 亚洲,一卡二卡三卡| 丁香六月天网| 纵有疾风起免费观看全集完整版| 久久综合国产亚洲精品| 国产深夜福利视频在线观看| 少妇人妻久久综合中文| 一边摸一边做爽爽视频免费| 欧美最新免费一区二区三区| 国产日韩一区二区三区精品不卡| 涩涩av久久男人的天堂| 一级片免费观看大全| 欧美日韩视频高清一区二区三区二| 亚洲美女搞黄在线观看| 亚洲av欧美aⅴ国产| 久久久久久久久免费视频了| 日韩视频在线欧美| 少妇人妻 视频| 老司机影院成人| 亚洲成人国产一区在线观看 | 一二三四中文在线观看免费高清| 久久国产精品男人的天堂亚洲| 婷婷色麻豆天堂久久| 日日摸夜夜添夜夜爱| 国产色婷婷99| 午夜福利网站1000一区二区三区| 久久久久精品久久久久真实原创| 日韩av免费高清视频| 超碰成人久久| 久久久久久人妻| 成人午夜精彩视频在线观看| 国产精品av久久久久免费| 亚洲色图 男人天堂 中文字幕| 国产亚洲最大av| 中文乱码字字幕精品一区二区三区| 亚洲av成人不卡在线观看播放网 | 欧美 亚洲 国产 日韩一| 妹子高潮喷水视频| 久久久精品免费免费高清| 久久精品国产亚洲av高清一级| 精品国产一区二区久久| 免费观看性生交大片5| 丰满饥渴人妻一区二区三| a级片在线免费高清观看视频| 亚洲精品视频女| 青青草视频在线视频观看| 黄色 视频免费看| 亚洲中文av在线| 久久久久精品久久久久真实原创| 久久影院123| 美女视频免费永久观看网站| 欧美av亚洲av综合av国产av | 国产av一区二区精品久久| 日韩av不卡免费在线播放| 操出白浆在线播放| 国产成人啪精品午夜网站| 在线免费观看不下载黄p国产| 午夜免费观看性视频| 日本色播在线视频| 啦啦啦视频在线资源免费观看| 另类精品久久| 久久亚洲国产成人精品v| 日韩一本色道免费dvd| 精品午夜福利在线看| av电影中文网址| 国产一级毛片在线| 亚洲,欧美精品.| 国产熟女欧美一区二区| a级毛片在线看网站| 国产av国产精品国产| 天天躁夜夜躁狠狠躁躁| 性高湖久久久久久久久免费观看| 99热国产这里只有精品6| 国精品久久久久久国模美| 在线观看www视频免费| 国产一卡二卡三卡精品 | 少妇人妻久久综合中文| 久久国产亚洲av麻豆专区| 日韩 亚洲 欧美在线| 另类亚洲欧美激情| 亚洲人成电影观看| 国产人伦9x9x在线观看| 久久精品久久久久久久性| 一本一本久久a久久精品综合妖精| 香蕉国产在线看| 看免费av毛片| 蜜桃国产av成人99| 亚洲成人国产一区在线观看 | 极品人妻少妇av视频| 中文字幕高清在线视频| 日日啪夜夜爽| 欧美激情高清一区二区三区 | 天天影视国产精品| av网站免费在线观看视频| 校园人妻丝袜中文字幕| 涩涩av久久男人的天堂| 亚洲精品aⅴ在线观看| 亚洲精品国产一区二区精华液| 美女福利国产在线| 在线免费观看不下载黄p国产| 免费看不卡的av| 国产精品久久久久久人妻精品电影 | 一级爰片在线观看| 大话2 男鬼变身卡| 欧美中文综合在线视频| 777久久人妻少妇嫩草av网站| 人妻 亚洲 视频| 欧美日韩成人在线一区二区| 人人澡人人妻人| 国产一级毛片在线| 看免费成人av毛片| 最新在线观看一区二区三区 | 亚洲精品av麻豆狂野| 久久久久久人人人人人| 这个男人来自地球电影免费观看 | 在线观看免费日韩欧美大片| 国产黄色视频一区二区在线观看| 天堂中文最新版在线下载| 啦啦啦视频在线资源免费观看| 另类亚洲欧美激情| 国产精品欧美亚洲77777| 国产日韩欧美亚洲二区| 国产麻豆69| 亚洲av综合色区一区| 人成视频在线观看免费观看| 我要看黄色一级片免费的| www.av在线官网国产| 国产亚洲最大av| 天天影视国产精品| 日本vs欧美在线观看视频| 韩国高清视频一区二区三区| 亚洲少妇的诱惑av| 超色免费av| av女优亚洲男人天堂| 伦理电影大哥的女人| 亚洲精品国产一区二区精华液| av在线app专区| 亚洲一卡2卡3卡4卡5卡精品中文| 超碰97精品在线观看| 欧美精品一区二区大全| 日本欧美国产在线视频| 黄频高清免费视频| 我要看黄色一级片免费的| 久久久国产精品麻豆| 成人影院久久| 一区二区日韩欧美中文字幕| 中文精品一卡2卡3卡4更新| 女人高潮潮喷娇喘18禁视频| 亚洲伊人久久精品综合| 美女视频免费永久观看网站| 久久99精品国语久久久| 97精品久久久久久久久久精品| 国产成人精品久久久久久| 好男人视频免费观看在线| 国产一区二区三区综合在线观看| 欧美久久黑人一区二区| 久久精品aⅴ一区二区三区四区| 亚洲一级一片aⅴ在线观看| 熟妇人妻不卡中文字幕| 少妇人妻 视频| 欧美成人精品欧美一级黄| 只有这里有精品99| 少妇人妻精品综合一区二区| 午夜日本视频在线| 亚洲一区中文字幕在线| 国产有黄有色有爽视频| 日韩免费高清中文字幕av| 女性生殖器流出的白浆| 少妇猛男粗大的猛烈进出视频| 久久亚洲国产成人精品v| 欧美另类一区| 国产精品久久久人人做人人爽| av又黄又爽大尺度在线免费看| 亚洲国产精品国产精品| av片东京热男人的天堂| 日本91视频免费播放| 极品人妻少妇av视频| 18禁观看日本| 99香蕉大伊视频| 妹子高潮喷水视频| 免费观看a级毛片全部| 久久久精品区二区三区| 日本av免费视频播放| 亚洲五月色婷婷综合| 成人毛片60女人毛片免费| 丰满少妇做爰视频| 中文字幕人妻丝袜制服| 国产精品久久久久成人av| 国产精品av久久久久免费| 午夜日韩欧美国产| 国产男女超爽视频在线观看| 一区二区三区乱码不卡18| 久久久国产精品麻豆| 黄色毛片三级朝国网站| 人人妻人人澡人人爽人人夜夜| 久久 成人 亚洲| 欧美激情极品国产一区二区三区| 国产免费视频播放在线视频| 久久女婷五月综合色啪小说| 天天躁日日躁夜夜躁夜夜| 人人妻人人澡人人爽人人夜夜| 久久 成人 亚洲| 国产在线免费精品| 亚洲婷婷狠狠爱综合网| 国产有黄有色有爽视频| www日本在线高清视频| 亚洲欧美精品自产自拍| 曰老女人黄片| 天天躁夜夜躁狠狠久久av| 曰老女人黄片| 午夜日本视频在线| 一边摸一边做爽爽视频免费| 中文欧美无线码| 韩国高清视频一区二区三区| 一区在线观看完整版| 国产精品久久久久久久久免| 久久久久精品久久久久真实原创| 国产乱人偷精品视频| 日本vs欧美在线观看视频| 国产乱人偷精品视频| 日本vs欧美在线观看视频| 精品亚洲乱码少妇综合久久| 亚洲一级一片aⅴ在线观看| 国产野战对白在线观看| 日韩中文字幕视频在线看片| 成年动漫av网址| 国产成人啪精品午夜网站| 久久精品aⅴ一区二区三区四区| 欧美中文综合在线视频| 亚洲成人国产一区在线观看 | 亚洲精品中文字幕在线视频| 日本午夜av视频| 欧美 日韩 精品 国产| 欧美日韩一级在线毛片| 人人妻,人人澡人人爽秒播 | 丰满迷人的少妇在线观看| 亚洲精华国产精华液的使用体验| 欧美日韩视频高清一区二区三区二| 亚洲精品美女久久久久99蜜臀 | 国产爽快片一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲av成人不卡在线观看播放网 | 少妇被粗大猛烈的视频| 久久人人97超碰香蕉20202| 少妇人妻 视频| 国产乱人偷精品视频| 成人国语在线视频| 男女无遮挡免费网站观看| 国产成人啪精品午夜网站| 国产精品久久久人人做人人爽| 伦理电影免费视频| 不卡视频在线观看欧美| 好男人视频免费观看在线| 国产熟女欧美一区二区| 亚洲欧美中文字幕日韩二区| 欧美亚洲 丝袜 人妻 在线| 亚洲精品在线美女| 操出白浆在线播放| 免费黄频网站在线观看国产| 一级片'在线观看视频| 精品视频人人做人人爽| 国产日韩一区二区三区精品不卡| 黄色视频不卡| 99久久99久久久精品蜜桃| 一级毛片黄色毛片免费观看视频| 下体分泌物呈黄色| 免费人妻精品一区二区三区视频| 午夜精品国产一区二区电影| 99久久精品国产亚洲精品| 一个人免费看片子| 免费av中文字幕在线| 熟女av电影| 校园人妻丝袜中文字幕| 黑丝袜美女国产一区| 国语对白做爰xxxⅹ性视频网站| 一级毛片黄色毛片免费观看视频| 老司机在亚洲福利影院| 久久性视频一级片| www.自偷自拍.com| 各种免费的搞黄视频| avwww免费| 久久ye,这里只有精品| 亚洲,欧美,日韩| 国产1区2区3区精品| 99国产综合亚洲精品| 老司机影院成人| 亚洲精品视频女| 婷婷色av中文字幕| 亚洲国产av影院在线观看| 丁香六月天网| 一级毛片电影观看| 日本黄色日本黄色录像| 日韩免费高清中文字幕av| 交换朋友夫妻互换小说| 国产一区二区三区综合在线观看| 看免费av毛片| 香蕉丝袜av| 欧美日韩精品网址| 日本av免费视频播放| 嫩草影视91久久| 亚洲成国产人片在线观看| 久久婷婷青草| 精品少妇久久久久久888优播| 日韩熟女老妇一区二区性免费视频| 午夜福利影视在线免费观看| 日韩av不卡免费在线播放| 另类精品久久| e午夜精品久久久久久久| 日韩精品免费视频一区二区三区| 成人漫画全彩无遮挡| 亚洲精品美女久久av网站| 色播在线永久视频| 男的添女的下面高潮视频| 亚洲欧美激情在线| 制服诱惑二区| 精品人妻一区二区三区麻豆| 黑人猛操日本美女一级片| 99国产精品免费福利视频| 亚洲国产欧美网| av有码第一页| 2021少妇久久久久久久久久久| 久久久久久久久久久久大奶| 飞空精品影院首页| 狂野欧美激情性xxxx| 两个人免费观看高清视频| 日本av手机在线免费观看| 国产av精品麻豆| 午夜激情av网站| 国产 一区精品| 涩涩av久久男人的天堂| 欧美少妇被猛烈插入视频| 老汉色av国产亚洲站长工具| 欧美人与性动交α欧美软件| 只有这里有精品99| 十八禁高潮呻吟视频| 亚洲欧洲日产国产| 19禁男女啪啪无遮挡网站| 国语对白做爰xxxⅹ性视频网站| 亚洲国产欧美一区二区综合| 黄色 视频免费看| 国产精品亚洲av一区麻豆 | 精品少妇内射三级| 国产成人精品久久二区二区91 | 日韩一本色道免费dvd| 晚上一个人看的免费电影| 狂野欧美激情性xxxx| 日韩电影二区| 亚洲色图综合在线观看| 亚洲av日韩在线播放| 精品国产超薄肉色丝袜足j| 亚洲伊人久久精品综合| 国产又色又爽无遮挡免| 午夜福利视频精品| 一级片免费观看大全| 中文字幕亚洲精品专区| 大码成人一级视频| 深夜精品福利| 欧美人与性动交α欧美软件| 成人手机av| 精品亚洲成a人片在线观看| 男女高潮啪啪啪动态图| 亚洲成色77777| 97人妻天天添夜夜摸| 亚洲一卡2卡3卡4卡5卡精品中文| 国产色婷婷99| 一区二区三区四区激情视频| 男女下面插进去视频免费观看| 亚洲精品中文字幕在线视频| 精品久久蜜臀av无| 在线精品无人区一区二区三| 夜夜骑夜夜射夜夜干| 超色免费av| av在线app专区| 国产高清国产精品国产三级| 黑丝袜美女国产一区| 精品少妇一区二区三区视频日本电影 | 国产成人av激情在线播放| 国产有黄有色有爽视频| 精品国产超薄肉色丝袜足j| 一本—道久久a久久精品蜜桃钙片| 亚洲视频免费观看视频| av国产久精品久网站免费入址| 亚洲国产欧美一区二区综合| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产中文字幕在线视频| 女人久久www免费人成看片| 亚洲国产精品成人久久小说| 免费在线观看黄色视频的| 精品一品国产午夜福利视频| 国产成人免费无遮挡视频| 最近最新中文字幕大全免费视频 | 伦理电影大哥的女人| 黄色视频在线播放观看不卡| bbb黄色大片| 久热这里只有精品99| 嫩草影院入口| 十八禁人妻一区二区| 国产亚洲av片在线观看秒播厂| 两个人看的免费小视频| 一级毛片黄色毛片免费观看视频| 免费看av在线观看网站| 啦啦啦在线观看免费高清www| 欧美日本中文国产一区发布| 下体分泌物呈黄色| 国产麻豆69| 国产一卡二卡三卡精品 | 日韩熟女老妇一区二区性免费视频| 97人妻天天添夜夜摸| 亚洲国产av新网站| 久久天堂一区二区三区四区| 日韩免费高清中文字幕av| 91aial.com中文字幕在线观看| tube8黄色片| 久久青草综合色| 中文天堂在线官网| 国产在线一区二区三区精| 欧美精品av麻豆av| 欧美日韩亚洲综合一区二区三区_| 日韩一区二区三区影片| 免费不卡黄色视频| 国产伦理片在线播放av一区| 超碰成人久久| 操出白浆在线播放| 国产精品国产av在线观看| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| 精品国产一区二区三区四区第35| 卡戴珊不雅视频在线播放| 这个男人来自地球电影免费观看 | 搡老乐熟女国产| 国产有黄有色有爽视频| 如何舔出高潮| 久久久精品94久久精品| 人体艺术视频欧美日本| 黄色 视频免费看| 国产精品 欧美亚洲| 最近中文字幕2019免费版| 蜜桃国产av成人99| 国产亚洲午夜精品一区二区久久| 国产成人精品在线电影| 亚洲熟女毛片儿| 中文字幕色久视频| 女性生殖器流出的白浆| 黄网站色视频无遮挡免费观看| 在线观看三级黄色| 成人亚洲精品一区在线观看| 一边摸一边做爽爽视频免费| xxx大片免费视频| 国产免费现黄频在线看| 在线观看www视频免费| 美女扒开内裤让男人捅视频| 国产男女超爽视频在线观看| 蜜桃国产av成人99| 欧美少妇被猛烈插入视频| 中文字幕人妻丝袜制服| 中文字幕制服av| 在线观看国产h片| 亚洲欧美精品自产自拍| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区 视频在线| 国产免费福利视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 久久综合国产亚洲精品| 电影成人av| 91成人精品电影| 啦啦啦视频在线资源免费观看|