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

    不同運行方式對高速地鐵氣動效應的影響

    2019-04-26 05:19:16冉騰飛梁習鋒熊小慧
    鐵道科學與工程學報 2019年4期
    關鍵詞:明線站間風井

    冉騰飛,梁習鋒,熊小慧

    ?

    不同運行方式對高速地鐵氣動效應的影響

    冉騰飛1, 2, 3,梁習鋒1, 2, 3,熊小慧1, 2, 3

    (1. 中南大學 交通運輸工程學院,軌道交通安全教育部重點實驗室,湖南 長沙 410075;2. 中南大學 軌道交通安全關鍵技術國際合作聯(lián)合實驗室,湖南 長沙 410075; 3. 中南大學 軌道交通列車安全保障技術國家地方聯(lián)合工程研究中心,湖南 長沙 410075)

    為了研究時速140 km/h高速地鐵列車以不同運行方式在隧道中運行時的氣動效應,采用三維、可壓、非定常N-S方程的數(shù)值計算方法,對地鐵列車由明線駛入隧道及站間運行時產(chǎn)生的氣動效應進行數(shù)值模擬,分析不同運行方式對高速地鐵隧道氣動效應的影響。研究結果表明:列車站間運行時,車體表面測點壓力峰峰值沿車長方向基本不變;而列車由明線駛入隧道時,車體表面測點壓力峰峰值從頭車向尾車逐漸降低。2種運行方式下的隧道壁面測點壓力峰峰值均在中間風井處達到最小值。并且列車由明線駛入隧道時的最大車體表面和隧道壁面壓力峰峰值分別為列車站間運行時的1.37倍與1.49倍。不同列車密封指數(shù)下,列車由明線駛入隧道時的車內(nèi)壓力變化均大于列車站間運行時的車內(nèi)壓力變化。因此,地鐵列車由明線駛入隧道時的空氣動力學效應比站間運行時更加不利。

    高速地鐵;運行方式;氣動效應;站間運行;中間風井;壓力波;數(shù)值模擬;密封指數(shù)

    地鐵隧道的運行環(huán)境十分復雜,在運行過程中,地鐵列車需要不斷地在車站間啟動加速與制動停車,并多次駛入與駛出隧道,以便在地下線路與高架線路間轉換。列車運行速度的改變會引起隧道內(nèi)壓力變化,隧道入口及中間風井處隧道截面的突變會引起隧道內(nèi)壓力劇烈波動。由于隧道內(nèi)壓力波動幅值和列車運行速度的平方近似成正比[1],因此列車運行速度的提高也會引起隧道內(nèi)壓力劇烈變化。當劇烈的壓力波動傳入車內(nèi)時,會引起乘客耳膜壓痛,甚至嘔吐,產(chǎn)生乘坐舒適性等問題[2]。因此有必要研究高速地鐵列車以不同方式在隧道中運行時的氣動效應。本文依托于下一代高速地鐵項目的研究,分別對時速140 km/h的下一代高速地鐵列車由明線駛入隧道與站間運行時引起的空氣動力學效應進行研究。數(shù)值計算和試驗是國內(nèi)外學者研究地鐵隧道空氣動力學問題的主要手段。試驗主要分為實車試驗和模型試驗2類,模型試驗分為風洞試驗、水槽試驗和動模型試驗等。但是,由于試驗具有經(jīng)濟成本高、操作不便、容易受客觀環(huán)境影響等特點,大部分研究人員采用數(shù)值計算的方法研究地鐵隧道空氣動力學。王秀珍[3]采用三維數(shù)值模擬方法,研究了車速、隧道截面形狀和隧道長度等因素對地鐵列車氣動效應的影響,沒有考慮列車站間運行方式對隧道內(nèi)壓力的影響,也沒有考慮地鐵隧道結構的復雜性和特殊性。祝嵐等[4]等利用三維、可壓縮和非定常的湍流模型,對S型地鐵列車勻速進、出隧道以及在隧道中加減速行駛進行數(shù)值計算。車輪飛[5]利用動網(wǎng)格和三維數(shù)值模擬方法計算地鐵列車通過中間風井時的隧道內(nèi)壓力波動,但是列車最大運行速度為100 km/h,且未考慮列車由明線駛入隧道及加減速的情況。劉伊江[6]計算地鐵列車在車站間勻速運行時的隧道內(nèi)風壓,未充分考慮中間風井所引起的氣動效應,也沒有考慮列車由明線駛入隧道的情況。NIU等[7]計算地鐵列車在隧道中的站間運行,卻忽略了中間風井對隧道內(nèi)氣動效應的影響。本文借鑒國內(nèi)外學者對高速鐵路隧道空氣動力學的研究方法及理論,結合地鐵隧道結構的復雜性和特殊性,對高速地鐵隧道空氣動力學進行研究。利用ICEM CFD對地鐵列車與隧道進行結構網(wǎng)格劃分,采用可壓、非定常N-S方程及?湍流方程的數(shù)值模擬方法及滑移網(wǎng)格技術,對時速140 km/h高速地鐵列車以不同方式在隧道中運行時的空氣動力學效應進行數(shù)值計算,得到隧道內(nèi)和車內(nèi)的壓力變化情況,研究不同運行方式對高速地鐵氣動效應的影響,并判斷哪種運行方式引起的空氣動力學效應更不利。

    1 運行區(qū)間的確定

    地鐵列車站間運行的速度圖如圖1(a)所示,運行過程分為起動加速、惰性運行和制動停車3個階段[8]。為確保數(shù)值模擬計算簡單、有效,可以將列車的站間運行過程看成勻加速運行、勻速運行和勻減速運行[7],如圖1(b)所示。其中啟動加速、勻速運行和制動減速的加速度分別為1=1 m/s2,2=0 m/s2,3=?1 m/s2。

    (a) 站間運行速度圖;(b) 站間運行簡化速度圖

    為了對比分析地鐵列車站間運行與勻速駛入隧道過程中的氣動效應,地鐵列車經(jīng)過中間風井時需保持勻速運行(即惰性運行)。因此,可以將地鐵列車的運行區(qū)間設置為最短運行區(qū)間的基礎上適當增加列車惰性運行區(qū)間,具體分析如下。

    當?shù)罔F列車的惰性運行時間為0時,可以計算得到列車站間運行的最短運行區(qū)間長度為:

    其中:1=1 m/s2,2=?1 m/s2,v=140 km/h=38.89 m/s,1=2=38.89 s,3=77.78 s。

    在最短運行區(qū)間的基礎上,適當增加12.54 s的惰性運行時間,可以計算得到地鐵隧道的運行區(qū)間長度為:

    其中:1=1 m/s2,2=?1 m/s2,v=140 km/h=38.89 m/s,1=38.89 s,2=51.43 s,3=90.32 s。

    2 計算條件

    采用滑移網(wǎng)格法模擬地鐵列車在隧道中的運行,計算區(qū)域如圖2所示。中間風井是在隧道中部預留的通道,主要用于通風、防水,危急時刻也可用于排煙、逃生與消防。為滿足地鐵隧道環(huán)境的通風換氣、環(huán)控和消防方面的要求,在地鐵隧道的普通長區(qū)間(大于等于1.8 km)中均需設置中間風井。本文中,列車在隧道中運行的區(qū)間長度為2 km,在普通長區(qū)間范圍內(nèi)。因此需要在隧道運行區(qū)間的中間位置設置矩形截面的中間風井[4],其中,風井壁面的尺寸為4 m×4 m×15 m,風井出口域的尺寸為50 m×50 m×15 m。

    為了真實模擬地鐵列車駛入隧道洞口的情況,保證列車周圍流場充分發(fā)展[9],設置車頭距離隧道入口50 m,車尾距離流域入口處150 m。

    單位:m

    地鐵隧道長度很大,隧道網(wǎng)絡復雜[6],因此,對地鐵隧道環(huán)境進行了簡化。為了模擬地鐵長大的隧道環(huán)境,避免出口邊界受尾流的影響[9],在地鐵隧道2 km運行區(qū)間的基礎上添加1 km長的隧道,如圖2所示;同時,為了讓隧道出口(列車站間運行時為隧道入口和隧道出口)免受反射波帶來的影響,將其設為壓力遠場P-far[7]。

    給定列車運動邊界條件:方向速度分量等于列車運行速度=140 km/h,和方向速度分量等于0。流域側面、頂面給定為對稱面symmetry,流域的底面、風井壁面、隧道壁面及地面均給定無滑移邊界條件wall。同時,為了使流場得到充分發(fā)展,將進口邊界條件設為in=0,將風井出口邊界條件設為out=0。

    數(shù)值計算采用的6車編組A型地鐵列車如圖3所示。圖中,對地鐵列車的外形進行了簡化,忽略轉向架、受電弓的影響[11?12]。其中,列車高=3.8 m,寬=3 m,整車長tr=139.5 m,列車橫截面面積v=9.785 m2。地鐵列車表面共布置10個測點,測點布置情況如圖4所示,其中頭尾車測點對稱布置。

    單位:m

    數(shù)值計算采用的隧道凈空斷面如圖4所示。圖中,數(shù)值計算采用隧道的盾構直徑為6 m,對應的凈空斷面面積A為26 m2,對應的阻塞比(即列車橫截面積與隧道凈空斷面面積之比)為0.376。

    隧道壁面壓力測點布置在距離地面3.94 m的隧道壁面上,并在中間風井兩側對稱布置,共有22個測點,其分布情況如圖4所示。

    在CFD計算中,計算網(wǎng)格質量的好壞會直接影響數(shù)值計算結果的效率與精度[13],因此,對隧道中的流體區(qū)域進行結構網(wǎng)格劃分:采用外C網(wǎng)格與O網(wǎng)格將列車壁面、隧道壁面和相應的BLOCK的surface進行映射,建立無厚度壁面;中間風井與隧道連接處的結構類似于T型管,可以劃分為三通網(wǎng)格,部分網(wǎng)格劃分如圖5所示。

    圖5 計算網(wǎng)格分布示意圖

    3 數(shù)值算法驗證

    為驗證數(shù)值計算結果的正確性與可信度,利用中南大學軌道交通安全教育部重點實驗室自主研制1:20縮比的動模型試驗裝置,對時速100 km/h的地鐵列車勻速通過盾構直徑為6 m,長度為1 km隧道時的場景進行動模型試驗,并進行相同工況的數(shù)值計算驗證。其中,動模型試驗裝置的列車和隧道模型如圖6所示,數(shù)值計算方法為可壓、非定常N-S方程與?湍流方程。

    圖6 動模型試驗模型

    動模型試驗中的測點布置示意圖如圖7所示。其中,隧道壁面上共布置6個測點,車體表面上共布置6個測點。

    單位:m

    數(shù)值計算與動模型試驗中的車體表面2號測點和隧道壁面3號測點的壓力對比波形如圖8所示;相應的測點壓力變化的結果對比如表1所示。

    從表1中數(shù)據(jù)可以看出,數(shù)值計算所得測點壓力的最大值、最小值和峰峰值與動模型試驗結果中最大值、最小值和峰峰值的最大差異為5.9%,在數(shù)值計算誤差容許范圍內(nèi),從而驗證了本文數(shù)值計算方法的正確性與可信度。

    (a) 車體表面2號測點;(b) 隧道壁面3號測點

    表1 壓力變化結果對比

    4 列車在隧道中運行時氣動效應研究

    分別對地鐵列車站間運行和由明線駛入隧道時的隧道內(nèi)和車內(nèi)氣動效應進行對比分析,研究不同運行方式對高速地鐵氣動效應的影響。

    圖9和圖10分別為列車站間運行和由明線駛入隧道時,不同時刻隧道內(nèi)的壓力分布云圖,其中列車的運行區(qū)間均為2 km。

    (a) t =2.1 s;(b) t =44.1 s;(c) t =48.3 s;(d) t =58.8 s

    從圖9可以看出,=2.1 s時,列車啟動加速時隧道內(nèi)壓力變化范圍很小,中間風井處壓力基本無變化。=44.1 s時,列車勻速運行至中間風井前側,列車速度的提高使隧道內(nèi)壓力變化范圍明顯增大,其中尾車的車頭與車身過渡處負壓值達到?1.7 kPa;中間風井處斷面突變,導致附近的壓力均在0.4 kPa以下。=48.3 s時,列車勻速通過中間風井,頭車經(jīng)過中間風井時會產(chǎn)生壓縮波,使隧道內(nèi)壓力增大,其中頭車的鼻尖以及相應的隧道壁面正壓值均達到2.5 kPa;中間風井斷面突變,導致附近的負壓至迅速減小至0.4 kPa以下;而尾車通過中間風井后會形成膨脹波,使隧道內(nèi)壓力下降。=58.8 s時,列車減速制動,列車速度的降低導致隧道內(nèi)壓力變化范圍變小,中間風井處受列車運動的影響很小,壓力基本無變化。

    (a) t =4.4 s;(b) t =26.4 s;(c) t =35.2 s;(d) t =58.8 s

    從圖10可以看出,=4.4 s時,列車頭部駛入隧道,頭車會形成初始壓縮波,使隧道內(nèi)壓力驟增。=26.4 s時,列車在到達中間風井前,由于中間風井斷面突變的影響,車體表面和隧道壁面的壓力變小,均在0.55 kPa以下。在=35.2 s,列車已經(jīng)通過中間風井,在頭車通過中間風井時會產(chǎn)生壓縮波,其以聲速在隧道內(nèi)傳播;在尾車通過中間風井后,列車尾部產(chǎn)生膨脹波,壓力下降。=58.8 s時,列車繼續(xù)在隧道中傳播,車體表面壓力變化范圍減小,并且頭車及隧道壁面的正壓值急劇降低。

    4.1 車體表面壓力對比

    當列車由明線駛入隧道時,由于壓力遠場的設置,隧道出口沒有反射波返回,壓力波在隧道入口和中間風井之間不斷反射。而列車站間運行時,隧道入口和隧道出口均沒有反射波返回,壓力波只能在中間風井處反射。

    本節(jié)選取車體表面2號測點,分析以不同運行方式對車體表面測點壓力變化的影響,結果如圖11所示。

    從圖11可以看出,不同運行方式對車體表面測點壓力變化的影響很大:在列車勻速駛入隧道的過程中,測點壓力的第1次上升①是由于列車頭部由明線駛入隧道時產(chǎn)生的初始壓縮波傳播到車體表面測點引起的。測點壓力的第1次下降②是由于列車尾部駛入隧道時,產(chǎn)生的初始膨脹波傳播到車體表面測點引起的。測點壓力波的第2次下降③是因為初始壓縮波在傳播時遇到中間風井時,一部分壓縮波以膨脹波的形式向隧道入口返回時傳播到車體表面測點引起的。這之后,隧道內(nèi)的壓縮波與膨脹波在隧道入口與中間風井之間不斷地反射,這些反射回來的壓力波傳播到車體表面測點,導致測點壓力不斷變化。列車經(jīng)過中間風井④時,車體表面測點壓力劇烈變化:頭車經(jīng)過中間風井時產(chǎn)生壓縮波,導致測點壓力迅速增大;列車經(jīng)過中間風井過程中,中間風井處氣流紊亂,導致車體表面壓力劇烈變化;尾車經(jīng)過中間風井時產(chǎn)生膨脹波,導致車體表面測點壓力迅速下降。這之后車體表面測點壓力變化逐漸趨于穩(wěn)定。

    圖11 不同運行方式下車體表面測點壓力變化曲線

    列車站間運行時,在加速運行過程⑤中,車體表面測點主要在負壓區(qū)波動,并且負壓值緩慢降低,這是由于列車啟動加速過程中,所引起的壓縮波不明顯,從而削弱了測點在正壓區(qū)的波動;列車在時刻⑥開始勻速運行,測點壓力基本不變;當列車經(jīng)過中間風井⑦時,車體表面壓力會迅速增大,并劇烈變化;列車通過中間風井后,測點壓力迅速減小至0附近,并保持穩(wěn)定。從時刻⑧開始,列車在隧道內(nèi)減速運行⑨,由于地鐵隧道較長且隧道出口無反射波返回,車體表面測點壓力經(jīng)歷較小的正壓波動后逐漸趨于平緩。

    不同運行方式下沿車長方向的壓力峰峰值變化規(guī)律如圖12所示。其中,距鼻尖距離為車體表面各測點距離鼻尖的水平距離與列車高度的比值,列車高度=3.8 m。從圖12可以看出,列車站間運行時,車體表面測點壓力峰峰值沿車長方向基本不變;而列車由明線駛入隧道時,車體表面測點壓力峰峰值沿車長方向從頭車向尾車呈下降趨勢;并且列車在駛入隧道時車體表面測點最大壓力峰峰值比列車站間運行時的車體表面測點壓力峰峰值大37%。因此,地鐵列車由明線駛入隧道時車體表面壓力變化超過列車站間運行時車體表面壓力變化。

    圖12 不同運行方式下車體表面測點壓力峰峰值變化規(guī)律

    4.2 隧道壁面壓力對比

    選取距離列車啟動位置500 m處的隧道壁面測點,分析以不同運行方式對隧道壁面測點壓力變化的影響,結果如圖13所示。

    從圖13可以看出,不同運行方式對隧道壁面測點壓力變化的影響很大:在列車勻速駛入隧道的過程中,測點壓力的第一次上升①,是由于列車頭部由明線駛入隧道時產(chǎn)生的初始壓縮波傳播到隧道壁面測點引起的。測點壓力第1次下降②,是引起的列車尾部駛入隧道時產(chǎn)生的膨脹波傳播到隧道壁面測點引起的。這之后,隧道內(nèi)的壓縮波與膨脹波在隧道入口與中間風井之間不斷地反射,這些反射回來的壓力波傳播到隧道壁面測點,導致測點壓力不斷變化。列車經(jīng)過隧道壁面測點③時,隧道壁面測點不斷受到列車帶來的壓縮波與膨脹波,測點壓力發(fā)生輕微抖動。列車經(jīng)過風井④時,隧道壁面測點壓力劇烈變化,這之后壓力變化逐漸趨于 平緩。

    圖13 不同運行方式下隧道壁面測點壓力變化曲線

    列車站間運行時,在加速運行過程⑤中,測點在正壓區(qū)產(chǎn)生微小波動,這是由于列車在車站啟動加速時,產(chǎn)生的壓縮波傳播到隧道壁面測點引起的;在列車經(jīng)過隧道壁面測點⑥時,隧道壁面測點壓力會發(fā)生輕微抖動;在列車在⑦處開始勻速運行,測點壓力基本不變化;列車經(jīng)過中間風井⑧時,隧道壁面壓力迅速增大,并發(fā)生劇烈變化,在列車通過中間風井后,測點壓力又迅速減小至0值附近,并保持穩(wěn)定。從時刻⑨開始,列車在隧道內(nèi)減速運行⑩,由于地鐵隧道較長且出口處無反射波,列車減速過程中無膨脹波反射回來,隧道測點壓力經(jīng)歷較小的負壓波動后逐漸趨于平緩。

    從圖14可以看出,列車以不同運行方式在隧道中運行時,受中間風井處斷面突變的影響,隧道壁面測點壓力峰峰值在中間風井附近達到最小值,最小值為1.47 kPa左右。列車由明線駛入隧道時,壓力峰峰值在距離隧道入口500 m和1 500 m附近達到最大值,為4.83 kPa。列車站間運行時,在距離列車啟動位置840 m和1 160 m附近達到最大值,為3.23 kPa。并且列車駛入隧道時的最大壓力峰峰值比站間運行時的最大壓力峰峰值大49%。因此,地鐵列車由明線駛入隧道時隧道壁面壓力變化超過列車站間運行時隧道壁面壓力變化。

    圖14 不同運行方式下隧道壁面測點壓力峰峰值變化規(guī)律

    4.3 更不利地鐵列車運行方式的選擇

    車內(nèi)壓力變化是衡量乘客乘坐舒適性的重要指標之一,也是判斷空氣動力學效應的重要標準。車內(nèi)壓力變化可以由車外壓力變化通過相應的迭代公式計算得到。

    假定車內(nèi)壓力變化率與車內(nèi)外壓差成正比:

    其中:p為車內(nèi)部壓力;p為車外部壓力;為列車密封指數(shù)。

    式(1)可以等效成

    其中:Δ為數(shù)值計算過程中的時間步長。

    在初始時刻=0 s時,p0=0,式(2)可以等效成

    式(3)即為車內(nèi)壓力變化的迭代公式。

    當?shù)罔F列車的密封指數(shù)分別為=0.5 s,=0.7 s,=1.5 s時,將列車以不同運行方式運行時測得的車體表面2號測點對應的車內(nèi)測點壓力變化率和3 s壓力變化幅值(即任意3 s內(nèi)壓力變化峰峰值)進行對比分析,結果如圖15所示。

    從圖15可以看出,地鐵列車以不同方式在隧道中運行時,車內(nèi)壓力變化差距很大。不同密封指數(shù)的地鐵列車以不同運行方式在隧道中運行時,車內(nèi)最大壓力變化率和3 s壓力變化幅值均沿車長方向從頭車向尾車逐漸減小。當密封指數(shù)分別為=0.5 s,=0.7 s,=1.5 s時,地鐵列車由明線駛入隧道時的車內(nèi)最大壓力變化率分別是列車站間運行時車內(nèi)壓力變化率的1.14倍,1.10倍和1.15倍;3 s壓力變化幅值分別是列車站間運行時的1.37倍,1.29倍和1.11倍。因此,地鐵列車由明線駛入隧道時車內(nèi)壓力變化超過列車站間運行時車內(nèi)壓力變化。

    綜上所述,地鐵列車由明線駛入隧道時的空氣動力學效應比站間運行時更不利。

    (a1), (a2) τ=0.5 s;(b1), (b2) τ=0.7 s;(c1), (c2) τ=1.5 s

    5 結論

    1) 列車站間運行時,車體表面壓力峰峰值沿車長方向基本不變;列車由明線駛入隧道時,車體表面壓力峰峰值從頭車向尾車逐漸降低。

    2) 不同運行方式下,隧道壁面測點壓力峰峰值均在中間風井處達到最小值。

    3) 列車由明線駛入隧道時的隧道內(nèi)壓力峰峰值和車內(nèi)壓力變化均超過列車站間運行時的隧道內(nèi)壓力峰峰值和車內(nèi)壓力變化。

    [1] 李新華. 高速列車過隧道、會車等典型狀態(tài)下非定??諝鈩恿W效應研究[D]. 杭州: 浙江大學, 2012. LI Xinhua. Unsteady simulation of a train passing a tunnel and two trains passing by each other[D]. Hangzhou: Zhejiang University, 2012.

    [2] 鈴木浩明, 高魁源. 車內(nèi)壓力波動引起耳鳴的研究[J]. 國外鐵道車輛, 1999(5): 15?18. Suzuki Hiroa, GAO Kuiyuan. A review of research trends on passengers’ aural discomfort caused by rail tunnel pressure change[J]. Foreign Rolling Stock, 1999(5): 15?18.

    [3] 王秀珍. 地鐵列車氣動效應分析[J]. 中國科技信息, 2011(21): 85?86.WANG Xiuzhen. The aerodynamic effect analysis of metro train[J]. China Science and Technology Information, 2011(21): 85?86.

    [4] 祝嵐, 張東, 孫振旭, 等. 基于乘客舒適性的快速地鐵隧道壓力波分析[J]. 都市快軌交通, 2015, 28(1): 87?91.ZHU Lan, ZHANG Dong, SUN Zhenxu, et al. Analysis on the pressure wave of rapid subway tunnel based on passenger comfort[J]. Urban Rapid Rail Transit, 2015, 28(1): 87?91.

    [5] 車輪飛. 城市地鐵隧道中間風井處車箱內(nèi)瞬變壓力模擬分析[J]. 建筑熱能通風空調(diào), 2011, 30(3): 89?91.CHE Lunfei. Simulation and analysis of the transient pressure variation of the carriage at the middle ventilation shaft of the urban subway tunnel[J]. Building Energy & Environment, 2011, 30(3): 89?91.

    [6] 劉伊江. 高速地鐵隧道壓力波研究及隧道斷面的擬定[J]. 鐵道標準設計, 2010(增2): 119?124.LIU Yijiang. Study on the pressure wave of high-speed subway tunnel and the design of tunnel section[J]. Railway Standard Design, 2010(Suppl 2): 119?124.

    [7] NIU J, ZHOU D, LIANG X, et al. Numerical study on the aerodynamic pressure of a metro train running between two adjacent platforms[J]. Tunnelling & Underground Space Technology, 2017, 65: 187?199.

    [8] 毛勵良. 地鐵車站合理站間距探討[J]. 都市快軌交通, 1996(3): 33?37. MAO Liliang.Reasonable station spacing of subway stations[J]. Urban Rapid Rail Transit, 1996(3): 33?37.

    [9] Ogawa T, Fujii K. Numerical investigation of three-dimensional compressible flows induced by a train moving into a tunnel[J]. Computers & Fluids, 1997, 26(6): 565?585.

    [10] González M L, Vega M G, Oro J M F, et al. Numerical modeling of the piston effect in longitudinal ventilation systems for subway tunnels[J]. Tunnelling & Underground Space Technology, 2014, 40(2): 22?37.

    [11] YANG X, JIN J, SHI G. Preliminary study on streamlined design of longitudinal profile of high-speed train head shape[J]. Procedia-Social and Behavioral Sciences, 2013, 96: 1469?1476.

    [12] Samane F. Aerodynamic of the trains in tunnels [D]. Birmingham: University of Birmingham, 2014.

    [13] 孫幫成, 李明高. ANSYS FLUENT 14.0仿真分析與優(yōu)化設計[M]. 北京: 機械工業(yè)出版社, 2014. SUN Bangcheng, LI Minggao. ANSYS FLUENT 14.0 simulation analysis and optimization design[M]. Beijing: China Mechanical Press, 2014.

    [14] CHU C R, CHIEN S Y, WANG C Y, et al. Numerical simulation of two trains intersecting in a tunnel[J]. Tunnelling & Underground Space Technology Incorporating Trenchless Technology Research, 2014, 42(5): 161?174.

    [15] XIANG X T, XUE L P. Tunnel hood effects on high speed train-tunnel compression wave[J]. Journal of Hydrodynamics Ser B, 2010, 22(5): 940?947.

    [16] WANG Yiwei, YANG Guowei, HUANG Chengguang, et al. Influence of tunnel length on the pressure wave generated by high-speed trains passing each other[J]. Science China Technological Sciences, 2012, 55(1): 255? 263.

    [17] 齊江浩, 趙蕾, 王君, 等. 地鐵隧道活塞風實測及特征分析[J]. 鐵道科學與工程學報, 2016, 13(4): 740?747.QI Jianghao, ZHAO Lei, WANG Jun, et al. Subway piston wind measurement and analysis of Xi’an metro Line 2[J]. Journal of Railway Science and Engineering, 2016, 13(4): 740?747.

    [18] 劉希賢, 郭安寧, 梅元貴, 等. 高速鐵路復線隧道內(nèi)壓力波特性的數(shù)值模擬研究[J]. 鐵道科學與工程學報, 2015, 12(1): 20?27. LIU Xixian, GUO Anning, MEI Yuangui, et al. Numerical simulation of pressure wave characteristics generated by a high-speed train passing through a double- track tunnel[J]. Journal of Railway Science and Engineering, 2015, 12(1): 20?27.

    Influence of different operation modes on the aerodynamic effect of high-speed subway

    RAN Tengfei1, 2, 3, LIANG Xifeng1, 2, 3, XIONG Xiaohui1, 2, 3

    (1.Key Laboratory of Traffic Safety on Track (Central South University), Ministry of Education, School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China; 2. Joint International Research Laboratory of Key Technology for Rail Traffic Safety, Central South University, Changsha 410075, China; 3. National & Local Joint Engineering Research Center of Safety Technology for Rail Vehicle, Central South University, Changsha 410075, China)

    To research the aerodynamic effect of high-speed trains running with different operation modes in tunnel whose speed is 140 km/h, it simulated aerodynamic problems caused by trains when they enter tunnel from open air and run between stations based on the three-dimensional, compressible and unsteady N-S equation method. The influence of aerodynamic effect with different trains’ operation modes in high-speed subway tunnel was analyzed. The results show: The peak-to-peak value of measurement points on train’s surface is constant along the direction of train length by trains entering tunnel from open air. But the peak-to-peak value of measurement points on train’s surface decreases from head carriage to rear carriage by trains running between stations. The peak-to-peak value of measurement points on tunnel wall is the minimum at the place of middle shaft by trains running with different operating modes. And the maximal peak-to-peak value on train’s surface and tunnel wall caused by trains entering tunnel from open air is respectively 1.37 times and 1.49 times compared with trains running between stations. The pressure changes inside carriages caused by trains entering tunnel from open air is larger than that caused by trains running between stations with different tightness index. As a result, the aerodynamic effect caused by trains entering tunnel from open air is worse than the aerodynamic effect caused by trains running between stations.

    high-speed subway; operation modes; aerodynamic effect; running between stations; middle shaft; pressure wave; numerical simulation; tightness index

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

    U270

    A

    1672 ? 7029(2019)04 ? 0860 ? 11

    2018?04?26

    國家科技支撐計劃項目(2015BAG12B01-24)

    梁習鋒(1963?),男,湖南長沙人,教授,從事列車空氣動力學研究;E?mail:gszxlxf@163.com

    (編輯 陽麗霞)

    猜你喜歡
    明線站間風井
    課堂教學的“明線”與“暗線”
    遼寧教育(2022年19期)2022-11-18 07:20:18
    地庫風井與景觀的和諧共生
    工程與建設(2019年3期)2019-10-10 01:40:20
    站間未設通過信號機的區(qū)間紅燈轉移問題探討
    行車對數(shù)和位置對活塞風井通風特性的影響
    一元微積分教學之思考
    商情(2017年6期)2017-04-18 00:51:17
    《拯救溺水魚》中的敘事內(nèi)容的明線與暗線
    求知導刊(2016年34期)2017-03-02 23:12:15
    《清水洗塵》的網(wǎng)狀結構分析
    文學教育(2017年1期)2017-02-10 23:00:23
    單線自動站間聯(lián)系電路的改進
    ZPW-2000A站間聯(lián)系電路的改進
    立井壁后注漿技術在白象山鐵礦風井的實踐
    日韩高清综合在线| 欧美一区二区国产精品久久精品| 九草在线视频观看| 高清在线视频一区二区三区 | 亚洲欧美日韩东京热| 国产精品不卡视频一区二区| 久久精品国产亚洲av涩爱 | 久久人人爽人人爽人人片va| АⅤ资源中文在线天堂| 亚洲av第一区精品v没综合| 欧美zozozo另类| 日日干狠狠操夜夜爽| 18禁在线无遮挡免费观看视频| 欧美xxxx黑人xx丫x性爽| 最近的中文字幕免费完整| 国产精品蜜桃在线观看 | 亚洲精品影视一区二区三区av| 日韩欧美精品免费久久| 久久人人爽人人片av| 精品久久久噜噜| 两性午夜刺激爽爽歪歪视频在线观看| 悠悠久久av| 国产午夜精品久久久久久一区二区三区| 草草在线视频免费看| 国产成人精品婷婷| 免费观看人在逋| 欧美成人一区二区免费高清观看| 国产色爽女视频免费观看| 国产成人精品一,二区 | 青春草亚洲视频在线观看| 午夜亚洲福利在线播放| 午夜爱爱视频在线播放| 国产黄a三级三级三级人| 99国产精品一区二区蜜桃av| 亚洲精品456在线播放app| 男女那种视频在线观看| 久久鲁丝午夜福利片| av.在线天堂| av在线观看视频网站免费| 欧美bdsm另类| 成年女人永久免费观看视频| 久久人人爽人人爽人人片va| 成人特级av手机在线观看| 国产一区二区三区在线臀色熟女| 日韩欧美精品免费久久| 亚洲成a人片在线一区二区| 国内久久婷婷六月综合欲色啪| 免费观看精品视频网站| 在线a可以看的网站| 又粗又爽又猛毛片免费看| 欧美日韩在线观看h| eeuss影院久久| 男的添女的下面高潮视频| 高清在线视频一区二区三区 | 欧美日韩在线观看h| 丰满的人妻完整版| 一级毛片久久久久久久久女| 毛片女人毛片| 又爽又黄无遮挡网站| 国产在线精品亚洲第一网站| 亚洲欧美日韩东京热| 1000部很黄的大片| www日本黄色视频网| 99热只有精品国产| 日韩成人av中文字幕在线观看| 能在线免费看毛片的网站| 久久久久性生活片| 大又大粗又爽又黄少妇毛片口| 国产毛片a区久久久久| 久久精品91蜜桃| 人人妻人人澡欧美一区二区| 亚洲国产精品sss在线观看| 成人特级av手机在线观看| 男插女下体视频免费在线播放| 色视频www国产| 久久久久久久亚洲中文字幕| 黄色配什么色好看| 亚洲真实伦在线观看| 国产精品免费一区二区三区在线| 日日摸夜夜添夜夜爱| 三级毛片av免费| 国产乱人偷精品视频| 日韩成人伦理影院| 天堂中文最新版在线下载 | 日本色播在线视频| 久久热精品热| 亚洲色图av天堂| 亚洲人成网站在线播放欧美日韩| 看非洲黑人一级黄片| 亚洲国产欧美人成| 亚洲高清免费不卡视频| 最近的中文字幕免费完整| 别揉我奶头 嗯啊视频| 国产成人一区二区在线| 国内精品美女久久久久久| 一区福利在线观看| 国产黄色小视频在线观看| 麻豆国产av国片精品| 亚洲久久久久久中文字幕| 中文字幕av在线有码专区| 国产精品日韩av在线免费观看| 好男人视频免费观看在线| 热99re8久久精品国产| 国产在视频线在精品| 成人永久免费在线观看视频| 成人毛片a级毛片在线播放| 热99在线观看视频| 国产精品一区二区三区四区久久| 69av精品久久久久久| 狠狠狠狠99中文字幕| 一本一本综合久久| 丝袜喷水一区| 黄色配什么色好看| 日韩一本色道免费dvd| 亚洲精品日韩在线中文字幕 | 国产单亲对白刺激| 成人欧美大片| 插逼视频在线观看| 日本-黄色视频高清免费观看| 熟女电影av网| 国语自产精品视频在线第100页| 熟妇人妻久久中文字幕3abv| 国产精品一二三区在线看| 欧美色视频一区免费| 一个人免费在线观看电影| 久久综合国产亚洲精品| 综合色丁香网| 久久综合国产亚洲精品| 偷拍熟女少妇极品色| 超碰av人人做人人爽久久| 18禁在线播放成人免费| 成人特级av手机在线观看| 亚洲无线在线观看| 69av精品久久久久久| 美女高潮的动态| 中文字幕久久专区| 黄色日韩在线| 国产三级在线视频| 九色成人免费人妻av| 国产三级在线视频| 小说图片视频综合网站| 一区二区三区免费毛片| 欧美精品国产亚洲| 我的老师免费观看完整版| 高清午夜精品一区二区三区 | 色尼玛亚洲综合影院| 国产在视频线在精品| 午夜精品一区二区三区免费看| 国产亚洲av片在线观看秒播厂 | 国产高潮美女av| 人妻夜夜爽99麻豆av| 久久久久久九九精品二区国产| 尾随美女入室| 观看美女的网站| 亚洲国产欧美在线一区| 久久久精品欧美日韩精品| 久久久久久久午夜电影| 搡女人真爽免费视频火全软件| 免费av不卡在线播放| 国产高清不卡午夜福利| 国产爱豆传媒在线观看| 日韩精品青青久久久久久| av免费观看日本| 3wmmmm亚洲av在线观看| 久久精品综合一区二区三区| 国产亚洲av嫩草精品影院| 成人美女网站在线观看视频| 精品少妇黑人巨大在线播放 | av天堂中文字幕网| 亚洲美女视频黄频| 国产蜜桃级精品一区二区三区| 草草在线视频免费看| 美女国产视频在线观看| 精品久久久久久久久av| 免费大片18禁| 18禁在线无遮挡免费观看视频| 老司机影院成人| 一边摸一边抽搐一进一小说| 午夜福利视频1000在线观看| 久久久国产成人免费| 一进一出抽搐动态| 欧美3d第一页| 人妻久久中文字幕网| 高清毛片免费看| 中出人妻视频一区二区| 99久久人妻综合| 人妻制服诱惑在线中文字幕| 99久久久亚洲精品蜜臀av| 亚洲av.av天堂| 国产黄片视频在线免费观看| 国产片特级美女逼逼视频| 男女那种视频在线观看| 日本成人三级电影网站| 亚洲av免费在线观看| 亚洲精品影视一区二区三区av| 亚洲一区高清亚洲精品| 日韩成人伦理影院| 亚洲激情五月婷婷啪啪| 亚洲最大成人av| 亚洲内射少妇av| 国产成人影院久久av| 国产大屁股一区二区在线视频| 亚洲色图av天堂| 国产不卡一卡二| 亚洲七黄色美女视频| 三级毛片av免费| 成人一区二区视频在线观看| 成人毛片60女人毛片免费| 久久久久久久久大av| 亚洲国产精品成人综合色| 亚洲成a人片在线一区二区| 亚洲av中文字字幕乱码综合| 99久久久亚洲精品蜜臀av| 婷婷色综合大香蕉| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 国产又黄又爽又无遮挡在线| 国语自产精品视频在线第100页| 26uuu在线亚洲综合色| 国产高清有码在线观看视频| 日韩亚洲欧美综合| 免费一级毛片在线播放高清视频| av.在线天堂| 爱豆传媒免费全集在线观看| 国模一区二区三区四区视频| 天堂影院成人在线观看| 国产女主播在线喷水免费视频网站 | 国产伦在线观看视频一区| 99久久九九国产精品国产免费| 日日啪夜夜撸| 免费观看人在逋| 一本久久中文字幕| a级毛色黄片| 欧美潮喷喷水| 两个人的视频大全免费| 亚洲中文字幕日韩| 国产蜜桃级精品一区二区三区| 一个人免费在线观看电影| 日本撒尿小便嘘嘘汇集6| 国产成人精品一,二区 | 久久久国产成人精品二区| 我要搜黄色片| 99久国产av精品国产电影| 国产精品久久久久久精品电影小说 | 久久久久久久久中文| 国产老妇伦熟女老妇高清| 欧美色视频一区免费| 精品免费久久久久久久清纯| 一级黄色大片毛片| 国产伦精品一区二区三区视频9| 久久人人爽人人爽人人片va| 国产精品精品国产色婷婷| 好男人视频免费观看在线| 性插视频无遮挡在线免费观看| 麻豆国产av国片精品| 日韩欧美国产在线观看| 亚洲综合色惰| 91久久精品国产一区二区三区| 啦啦啦啦在线视频资源| 婷婷亚洲欧美| 99久久中文字幕三级久久日本| 久久午夜亚洲精品久久| 欧美潮喷喷水| 欧美日韩在线观看h| 一夜夜www| 亚洲精品乱码久久久v下载方式| 夜夜看夜夜爽夜夜摸| 国产成人福利小说| 国产成人aa在线观看| 中文资源天堂在线| 美女大奶头视频| 欧美在线一区亚洲| 高清毛片免费看| 亚洲人成网站在线播放欧美日韩| 有码 亚洲区| 夜夜夜夜夜久久久久| 黄色配什么色好看| 日韩一本色道免费dvd| 国产精品一区二区性色av| 亚洲最大成人av| 欧美另类亚洲清纯唯美| 国产精品1区2区在线观看.| 国产精品国产三级国产av玫瑰| 一本一本综合久久| 久久久欧美国产精品| av在线观看视频网站免费| 99九九线精品视频在线观看视频| 人妻系列 视频| 伦理电影大哥的女人| 免费av毛片视频| 国产淫片久久久久久久久| 亚洲成av人片在线播放无| 久久99蜜桃精品久久| 久久久色成人| 99视频精品全部免费 在线| 一夜夜www| 亚洲成av人片在线播放无| 99在线人妻在线中文字幕| 亚洲天堂国产精品一区在线| 99久久精品热视频| 亚洲精品乱码久久久v下载方式| 国产真实伦视频高清在线观看| 久久国内精品自在自线图片| 老熟妇乱子伦视频在线观看| 欧美激情在线99| 国产伦一二天堂av在线观看| 中国国产av一级| 亚洲三级黄色毛片| 亚洲精品色激情综合| 久久久久久久久久久丰满| 深夜a级毛片| 青春草国产在线视频 | 国产女主播在线喷水免费视频网站 | 国产伦理片在线播放av一区 | 亚洲天堂国产精品一区在线| 久久精品人妻少妇| 欧美在线一区亚洲| 久久久精品大字幕| 成人三级黄色视频| 在线观看午夜福利视频| 国产成人午夜福利电影在线观看| 婷婷六月久久综合丁香| 成人亚洲精品av一区二区| av在线观看视频网站免费| 国产成人aa在线观看| 麻豆久久精品国产亚洲av| 午夜老司机福利剧场| 免费看光身美女| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线观看播放| 国产私拍福利视频在线观看| 夜夜夜夜夜久久久久| 免费无遮挡裸体视频| 久久婷婷人人爽人人干人人爱| 1024手机看黄色片| 91精品一卡2卡3卡4卡| 久久精品夜夜夜夜夜久久蜜豆| 亚洲婷婷狠狠爱综合网| 99久久精品一区二区三区| 成人三级黄色视频| 亚洲高清免费不卡视频| 免费观看的影片在线观看| 国产色婷婷99| 国产高潮美女av| 欧美日韩乱码在线| 99热这里只有是精品50| 自拍偷自拍亚洲精品老妇| av卡一久久| 小蜜桃在线观看免费完整版高清| 黑人高潮一二区| 精品久久久久久久久久久久久| 亚洲欧洲国产日韩| 高清在线视频一区二区三区 | 天堂√8在线中文| 十八禁国产超污无遮挡网站| 波多野结衣高清作品| 亚洲在线自拍视频| 五月玫瑰六月丁香| 中文字幕久久专区| 男女边吃奶边做爰视频| 18+在线观看网站| 国产视频首页在线观看| 91狼人影院| 18禁在线播放成人免费| 国产成人freesex在线| 欧美性猛交╳xxx乱大交人| 少妇的逼水好多| 少妇的逼好多水| 成人一区二区视频在线观看| 九九热线精品视视频播放| 国产欧美日韩精品一区二区| 国产私拍福利视频在线观看| 国产日韩欧美在线精品| 国产国拍精品亚洲av在线观看| 精品久久久久久久久久免费视频| 变态另类丝袜制服| 国产激情偷乱视频一区二区| av国产免费在线观看| 成人性生交大片免费视频hd| 天天躁夜夜躁狠狠久久av| 亚洲国产色片| 国产乱人视频| 成年免费大片在线观看| 久久国产乱子免费精品| 内地一区二区视频在线| 国产精品,欧美在线| 国产在视频线在精品| av免费在线看不卡| 国产探花在线观看一区二区| 国产精品爽爽va在线观看网站| 小蜜桃在线观看免费完整版高清| 国产精品女同一区二区软件| 岛国毛片在线播放| 能在线免费看毛片的网站| 全区人妻精品视频| av女优亚洲男人天堂| 在线观看av片永久免费下载| 午夜免费激情av| 麻豆成人av视频| 一进一出抽搐gif免费好疼| 黄片wwwwww| 熟女人妻精品中文字幕| 三级男女做爰猛烈吃奶摸视频| 欧美日韩乱码在线| 久久这里有精品视频免费| www.色视频.com| 桃色一区二区三区在线观看| 国产极品精品免费视频能看的| 成人午夜高清在线视频| 欧美激情国产日韩精品一区| 久久久久久国产a免费观看| 高清在线视频一区二区三区 | 久久精品国产亚洲av天美| av.在线天堂| 国产av不卡久久| 国产v大片淫在线免费观看| 亚洲成人久久性| 免费搜索国产男女视频| 91精品一卡2卡3卡4卡| 有码 亚洲区| 亚洲在线观看片| 亚洲国产欧美在线一区| 久久久久久久亚洲中文字幕| 日韩制服骚丝袜av| 国产久久久一区二区三区| 欧美日韩在线观看h| 久久99热6这里只有精品| 成人毛片60女人毛片免费| 久久午夜亚洲精品久久| 女同久久另类99精品国产91| 日韩成人av中文字幕在线观看| 久久久久久久久中文| 床上黄色一级片| 天堂网av新在线| 99久久中文字幕三级久久日本| 亚洲一区高清亚洲精品| 91久久精品电影网| 午夜福利在线观看吧| 国产精品99久久久久久久久| 久久6这里有精品| 网址你懂的国产日韩在线| 97热精品久久久久久| 亚洲电影在线观看av| 国产成人aa在线观看| 久久99蜜桃精品久久| 亚洲人成网站在线播| 免费av毛片视频| 联通29元200g的流量卡| 最近视频中文字幕2019在线8| 成人永久免费在线观看视频| 精品久久久久久久末码| 欧美精品一区二区大全| 色视频www国产| 久久精品夜夜夜夜夜久久蜜豆| 国产综合懂色| 综合色av麻豆| 少妇熟女aⅴ在线视频| 天天躁日日操中文字幕| 网址你懂的国产日韩在线| 黑人高潮一二区| 天堂影院成人在线观看| 黄片wwwwww| 日韩一区二区三区影片| 九九在线视频观看精品| 成人毛片a级毛片在线播放| 国产亚洲精品久久久com| 日韩一本色道免费dvd| 男人狂女人下面高潮的视频| 精品国产三级普通话版| 精品久久久久久久久久免费视频| 国产精品爽爽va在线观看网站| 国内精品一区二区在线观看| 能在线免费看毛片的网站| 一区二区三区高清视频在线| 干丝袜人妻中文字幕| avwww免费| 国产大屁股一区二区在线视频| 亚洲国产欧美人成| 亚洲最大成人中文| 一边摸一边抽搐一进一小说| 国产爱豆传媒在线观看| 蜜桃亚洲精品一区二区三区| 久久精品国产清高在天天线| 成人午夜精彩视频在线观看| а√天堂www在线а√下载| 亚洲人成网站在线播放欧美日韩| 热99re8久久精品国产| eeuss影院久久| 99热这里只有是精品50| 欧美性猛交╳xxx乱大交人| 国语自产精品视频在线第100页| 亚洲欧美日韩高清专用| 欧美日韩在线观看h| 99热这里只有是精品在线观看| 亚洲真实伦在线观看| 国产蜜桃级精品一区二区三区| 免费观看的影片在线观看| 久久久久久久久久久免费av| 一夜夜www| 欧美人与善性xxx| 成人漫画全彩无遮挡| 在线播放无遮挡| 欧美性猛交╳xxx乱大交人| 91精品一卡2卡3卡4卡| 亚洲最大成人av| 在线播放无遮挡| 深爱激情五月婷婷| 亚洲性久久影院| 美女xxoo啪啪120秒动态图| 精品久久久久久久人妻蜜臀av| 国产视频首页在线观看| 国语自产精品视频在线第100页| av女优亚洲男人天堂| 免费观看a级毛片全部| 亚洲不卡免费看| 成年女人永久免费观看视频| 老熟妇乱子伦视频在线观看| 成人av在线播放网站| 亚洲欧美精品专区久久| 99久久成人亚洲精品观看| 九九热线精品视视频播放| 亚洲经典国产精华液单| 舔av片在线| 又黄又爽又刺激的免费视频.| 婷婷精品国产亚洲av| 亚洲四区av| 国产一级毛片七仙女欲春2| 91久久精品国产一区二区成人| 在线a可以看的网站| 99久国产av精品国产电影| 国产高清不卡午夜福利| 高清毛片免费观看视频网站| 国产三级在线视频| 别揉我奶头 嗯啊视频| 变态另类丝袜制服| 国产精品一区二区三区四区免费观看| 男人和女人高潮做爰伦理| 村上凉子中文字幕在线| 国产成人a区在线观看| 欧美日韩乱码在线| 国内精品美女久久久久久| 嘟嘟电影网在线观看| 亚洲国产欧洲综合997久久,| 久久久久国产网址| 午夜爱爱视频在线播放| 搡女人真爽免费视频火全软件| 啦啦啦韩国在线观看视频| 99国产精品一区二区蜜桃av| 22中文网久久字幕| 色尼玛亚洲综合影院| 国产一区二区三区在线臀色熟女| 国产精品1区2区在线观看.| 亚洲va在线va天堂va国产| 亚洲中文字幕一区二区三区有码在线看| 69av精品久久久久久| 中文字幕精品亚洲无线码一区| 国产乱人视频| 人妻系列 视频| 久久人人精品亚洲av| 日本黄色片子视频| 村上凉子中文字幕在线| 婷婷色av中文字幕| 男女做爰动态图高潮gif福利片| 国内久久婷婷六月综合欲色啪| 日日撸夜夜添| 天堂网av新在线| 色综合站精品国产| 久久久久久久久久成人| 免费人成视频x8x8入口观看| 九九在线视频观看精品| 国内精品宾馆在线| 久久久久久九九精品二区国产| 啦啦啦啦在线视频资源| 欧美bdsm另类| 乱码一卡2卡4卡精品| 国产成人精品婷婷| 成人午夜精彩视频在线观看| 国产成人freesex在线| 日本免费a在线| 在线a可以看的网站| 最近的中文字幕免费完整| 最新中文字幕久久久久| 啦啦啦观看免费观看视频高清| 青青草视频在线视频观看| 国产毛片a区久久久久| 老司机影院成人| 男女啪啪激烈高潮av片| 性欧美人与动物交配| 99热网站在线观看| 国产欧美日韩精品一区二区| 综合色丁香网| 国产精品一区二区三区四区免费观看| 日本一二三区视频观看| 亚洲精品亚洲一区二区| 午夜福利成人在线免费观看| 97超视频在线观看视频| 精品国内亚洲2022精品成人| 国产精品爽爽va在线观看网站| 久久鲁丝午夜福利片| 一区二区三区四区激情视频 | 少妇丰满av| 一本精品99久久精品77| 色哟哟·www| 久久久久久久久久成人| 午夜老司机福利剧场| 1000部很黄的大片| 日本撒尿小便嘘嘘汇集6| 日韩精品有码人妻一区| 国产成人a区在线观看| 黄色视频,在线免费观看| 岛国毛片在线播放| 春色校园在线视频观看| 日韩av不卡免费在线播放| 亚洲欧美精品综合久久99| 97超视频在线观看视频|