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

    地鐵列車通過隧道時的氣動性能研究*

    2016-12-15 03:06:19徐世南張繼業(yè)
    城市軌道交通研究 2016年9期
    關(guān)鍵詞:尾車頭車峰峰

    徐世南 張繼業(yè) 熊 駿 孟 添

    (西南交通大學(xué)牽引動力國家重點實驗室,610031,成都//第一作者,碩士研究生)

    ?

    地鐵列車通過隧道時的氣動性能研究*

    徐世南 張繼業(yè) 熊 駿 孟 添

    (西南交通大學(xué)牽引動力國家重點實驗室,610031,成都//第一作者,碩士研究生)

    列車通過隧道時引起的空氣動力效應(yīng)會對列車運行的安全性、乘客乘坐的舒適性等產(chǎn)生不良影響。基于列車空氣動力學(xué)理論,采用計算流體力學(xué)軟件FLUENT對某型號地鐵車輛通過最不利長度隧道時的空氣動力學(xué)性能進(jìn)行數(shù)值模擬,得到并分析了地鐵列車和隧道壁面監(jiān)測點的壓力時程曲線和分布特征。研究表明:車體表面壓力峰峰值、3 s內(nèi)車內(nèi)壓力波動最大值及隧道內(nèi)附屬物壓力峰峰值,與列車速度的平方近似成線性關(guān)系;隧道斷面凈空面積越小,車體承受的壓力越大;地鐵列車通過隧道時需限速,以達(dá)到人體舒適性評價標(biāo)準(zhǔn)。

    地鐵列車; 隧道; 空氣動力效應(yīng); 壓力時程曲線

    Author′s address Traction Power State Key Laboratory,Southwest Jiaotong University,610031,Chengdu,China

    列車通過隧道時,會引起隧道內(nèi)和車體內(nèi)外空氣壓力急劇變化,引起出口微氣壓波、列車風(fēng)等空氣動力效應(yīng),這將會對列車運行的安全性、乘客乘坐的舒適性等產(chǎn)生嚴(yán)重影響。因此,各國學(xué)者對列車隧道耦合空氣動力學(xué)問題進(jìn)行了比較深入的研究。

    文獻(xiàn)[1-2]研究了高速列車通過隧道時產(chǎn)生的微壓波和瞬態(tài)壓力;文獻(xiàn)[3]研究了高速列車通過大長隧道時如何減輕空氣阻力和水波特性;文獻(xiàn)[4-7]開展了許多實車或模型試驗,得到列車隧道耦合空氣動力學(xué)數(shù)據(jù);文獻(xiàn)[8]對列車通過隧道時隧道壁面、車體表面和車廂內(nèi)部的耦合關(guān)系做了研究;文獻(xiàn)[9-10]研究了高速動車組的隧道壓力波特性;文獻(xiàn)[11]研究了車速、隧道長度等對隧道洞口微氣壓波的影響;文獻(xiàn)[12]進(jìn)行了地鐵列車隧道內(nèi)運行的氣動力學(xué)試驗與仿真。

    我國城市軌道交通建設(shè)速度很快,但是,地鐵列車在隧道尤其是在最不利長度隧道內(nèi)行駛的相關(guān)研究文獻(xiàn)較少,很多研究主要局限于空氣動力效應(yīng)的某單一方面分析。此外,由于地鐵列車頭部為鈍型,且地鐵隧道界面遠(yuǎn)小于高鐵隧道;因此,雖然地鐵列車行駛速度不高,但空氣動力學(xué)效應(yīng)明顯,會產(chǎn)生很多不良影響。如地鐵列車前方空氣受到壓縮而劇烈波動形成壓力波,則當(dāng)車內(nèi)壓力波動超過一定值后會使乘客產(chǎn)生耳膜疼痛[13-14]。所以,對地鐵列車通過最不利長度隧道時的氣動性能進(jìn)行研究非常重要。本文采用數(shù)值模擬,對某型號國產(chǎn)地鐵車輛在不同車速下通過對應(yīng)的最不利長度的隧道進(jìn)行了三維模擬,分析了地鐵車輛車內(nèi)外壓力和壁面壓力,對地鐵車輛安全性進(jìn)行了分析并給出相關(guān)建議。

    1 計算模型和網(wǎng)格劃分

    1.1 數(shù)學(xué)模型

    地鐵列車快速通過隧道時引起周圍空氣的流動是粘性、可壓縮、非定常的三維湍流流動,它需要遵循質(zhì)量守恒定律、動量守恒定律和能量守恒定律這三大定律[15]。通過求解三維瞬態(tài)可壓縮的雷諾時均N-S和k-ε兩個方程湍流模型獲得列車通過隧道時底板的氣動壓力,其控制方程的運輸方程[16]形式為:

    式中:

    t——列車運行時間;

    φ——流場通量;

    u——流場速度矢量;

    ut——列車運動速度矢量;

    ρ——擴(kuò)散系數(shù);

    S——源項。

    1.2 幾何模型

    列車計算模型采用某國產(chǎn)地鐵車輛,列車為6節(jié)編組,即頭車、尾車和4節(jié)中間車。頭、尾車長度均為21.26 m,中間車長度為20.00 m,列車總長為122.52 m,車寬為2.80 m,車高為2.92 m,地鐵車輛最大橫斷面面積為8.51 m2;忽略車體外部復(fù)雜結(jié)構(gòu)細(xì)節(jié),如受電弓、轉(zhuǎn)向架等。列車幾何模型如圖1所示。

    圖1 列車幾何模型示意圖

    1.3 計算區(qū)域及網(wǎng)格劃分

    隧道計算模型采用某城市快速軌道交通線工程中的兩種矩形斷面的單線隧道:一種隧道高為5.876 m,寬為4.4 m,隧道斷面凈空面積為26.15 m2,標(biāo)記為TB隧道;另一種隧道高為4.676 m,寬為4.4 m,隧道斷面凈空面積為22.23 m2,標(biāo)記為TS隧道。

    列車隧道通過計算區(qū)域如圖2所示。整個計算區(qū)域長Lt+692 m(Lt為最不利隧道長度)、寬120 m、高60 m,列車初始位置距離隧道入口50 m。圖中Lv為地鐵車輛長度。由文獻(xiàn)[17]可知,當(dāng)列車速度分別為120、110、100、90、80 km/h時,對應(yīng)的最不利隧道長度分別為990、1 035、1 100、1 165和1 250 m。

    圖2 列車隧道通過計算區(qū)域

    流場邊界條件的具體設(shè)置如圖3所示。外流場邊界A、B、C、D為壓力出口邊界條件,交界處F、E為Interface邊界條件,地面、隧道壁面和隧道進(jìn)出口面G、H按光滑壁面處理,給定無滑移邊界條件。

    圖3 列車-隧道流場邊界圖

    采用網(wǎng)格劃分工具ICEM CFD劃分網(wǎng)格,固定區(qū)域采用結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分。由于高速列車轉(zhuǎn)向架區(qū)域結(jié)構(gòu)較為復(fù)雜,因此移動區(qū)域采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分。

    1.4 壓力測點分布

    為了研究列車在通過隧道過程中車體表面壓力波的變化,分別在頭車、中間車、尾車表面取如圖4所示的壓力測點。其中,測點head1、tail1位于頭車和尾車鼻尖處,測點head2、tail2位于頭車和尾車擋風(fēng)玻璃處,測點head3、tail3位于頭車和尾車側(cè)翼處,head4、mid1、mid4、tail4分別位于列車頭、中、尾車側(cè)面車窗處,head5、head6、mid2、mid3、mid5、mid6、tail5、tail6分別位于列車頭車、中車、尾車的車頂及車底。

    以990 m長隧道為例介紹隧道測點布置。測點在列車近壁面處位于高度3 m、離車體0.6 m位置:在隧道內(nèi)距離隧道入口和出口10、20、50 m處各布置3個測點,在隧道內(nèi)中間位置和中間位置兩側(cè)100 m處布置3個測點,測點示意圖和測點名稱如圖5所示。

    圖4 測點分布示意圖

    圖5 990 m長隧道測點分布示意圖

    2 計算結(jié)果分析

    分別計算地鐵列車以速度120、110、100、90、80 km/h通過TB隧和TS隧道時各測點的壓力,得到各測點壓力隨時間變化的曲線。

    2.1 計算結(jié)果驗證

    為驗證計算模型和計算方法的正確性,將數(shù)值計算結(jié)果與實測數(shù)據(jù)進(jìn)行對比。實測數(shù)據(jù)是列車以350 km/h速度通過線間距為5.0 m,長度為2 989 m的復(fù)線隧道時某型號高速列車某車窗測點的壓力時程曲線。實車試驗與數(shù)值計算壓力時程曲線如圖6所示。由圖6可以看出:壓力波動規(guī)律基本一致,波峰、波谷、峰峰值誤差都比較小,所以,數(shù)值計算基本滿足計算精度。

    列車表面壓力變化都是受到馬赫波傳播的影響,所以各個監(jiān)測點壓力的變化規(guī)律相同。以120 km/h速度的地鐵列車通過TB隧道時頭車側(cè)面車窗處的測點為例,分析其壓力變化與馬赫波的關(guān)系,如圖7所示。

    圖6 列車某車窗測點壓力時程曲線

    根據(jù)隧道內(nèi)瞬態(tài)壓力形成機(jī)理,列車頭部進(jìn)入隧道產(chǎn)生壓縮波,尾部進(jìn)入隧道產(chǎn)生膨脹波,這兩種波在隧道口會產(chǎn)生反射。當(dāng)列車前行時,列車頭車鼻尖處監(jiān)測點首先遇到①處的膨脹波,此時壓力開始下降;隨著列車前行,該處監(jiān)測點遇到②處的膨脹波,壓力再次減??;之后遇到③處的壓縮波,監(jiān)測點壓力開始增加;之后,隨著監(jiān)測點遇到不同位置處的壓縮波和膨脹波,監(jiān)測點處壓力相應(yīng)增大和減小。由此可見,該監(jiān)測點壓力變化趨勢滿足波的傳播理論,驗證了計算模型規(guī)律的正確性。

    2.2 車體表面壓力幅值分析

    地鐵列車以120 km/h速度通過TS和TB兩種隧道時,列車車頭側(cè)面車窗測點(head4)的壓力分布如圖8所示。從圖9可以看出,列車以相同速度通過兩種隧道截面時的壓力變化曲線的變化規(guī)律一致;通過TS隧道承受的壓力正峰值、負(fù)峰值以及峰峰值均比TB隧道的大。這是因為當(dāng)列車橫截面積一定時,隧道凈空面積越小,其阻塞比越大,列車承受的壓力波動越大。因此,隧道截面越小,車體壁面壓力峰值及幅值越大。

    圖7 地鐵列車頭車鼻尖監(jiān)測點壓力變化過程

    圖8 地鐵列車以120 km/h車速通過990 m長

    以地鐵列車以120 km/h速度通過990 m長TB隧道情況為例,頭車、中間車和尾車的各個測點壓力變化規(guī)律如圖9所示。由圖9可知,當(dāng)列車進(jìn)入隧道時,列車前端及周圍空氣受到隧道壁面的約束,使列車前端的空氣受到劇烈壓縮,空氣壓力驟然增加;隨著列車逐漸進(jìn)入隧道,壓力繼續(xù)增加,但幅度稍平緩,直到整列車完全進(jìn)入隧道。最大正壓位于頭車鼻尖處,達(dá)到1 600 Pa;負(fù)壓較大位置位于尾車側(cè)翼處,出現(xiàn)在中間時刻,達(dá)到1 227 Pa。此外,非流線型部位各測點的壓力峰峰值差別不大。列車測點的波動是由于列車在隧道中運行時會遇到壓縮波和膨脹波,前者導(dǎo)致壓力增加,后者使壓力減小;隨著壓縮波和膨脹波在隧道內(nèi)不斷反射及改變,會逐漸減弱,因此波動幅度逐漸減小。

    以頭車鼻尖處測點(head1)為例,分析其在不同車速下,分別通過TS和TB兩種隧道時,測點壓力峰峰值與速度情況,如圖10所示。由圖10可知,在相同速度下通過TS隧道的壓力幅值大于TB隧道;對于同一隧道而言,列車運行速度越大,壓力幅值越大,且壓力幅值與速度平方近似成線性關(guān)系。通過曲線擬合可得:當(dāng)?shù)罔F列車通過TS隧道時,y=0.08x2.13,其相關(guān)性系數(shù)為0.98;地鐵列車通過TB隧道時,y=0.06x2.12,其相關(guān)性系數(shù)為0.99。

    圖9 地鐵列車以120 km/h速度通過990 m長TB隧道壓力波動

    圖10 不同速度下車體表面壓力幅值

    2.3 車內(nèi)壓力幅值分析

    地鐵列車通過單線隧道時車內(nèi)壓力的變化率與車內(nèi)外壓差成正比,其表達(dá)式[8]為:

    式中:

    Pout——車外壓力;

    Pin——車內(nèi)壓力;

    τ——密封指數(shù),參考?xì)W洲標(biāo)準(zhǔn),對于地鐵車輛而言,一般選取τ=0.5 s。

    根據(jù)數(shù)值仿真的車外壓力和上式對車內(nèi)壓力進(jìn)行數(shù)值求解。當(dāng)求解車內(nèi)壓力時,選取車外壓力的參考點為頭車側(cè)面車窗附近,以120 km/h速度等級通過990 m長隧道時內(nèi)外壓力波動情況為例,如圖11所示。由圖11可知,由于車體的密封性能夠緩解和滯后壓力變化,因此車內(nèi)壓力的波動稍有滯后且小于車外的壓力波動;車內(nèi)壓力波動最大幅值和3 s內(nèi)車內(nèi)壓力波動最大值均為1 215 Pa。

    圖11 地鐵列車通過990 m長隧道時車內(nèi)外壓力波動情況

    根據(jù)人體舒適性評判標(biāo)準(zhǔn),列車通過單線隧道時車內(nèi)壓力波動應(yīng)該小于等于800 Pa/3 s[18],所以不同車速下對應(yīng)的車內(nèi)壓力波動最大值示意圖如圖12所示。根據(jù)評判規(guī)則,在密封指數(shù)為0.5 s時,該地鐵列車在通過TB與TS隧道時,車速要分別小于100 km/h和90 km/h,才符合相關(guān)標(biāo)準(zhǔn)。

    圖12 列車通過隧道時不同速度下3 s內(nèi)車內(nèi)壓力

    通過曲線擬合,可得列車通過隧道時3 s內(nèi)車內(nèi)壓力波動最大值與速度平方近似呈線性關(guān)系:當(dāng)列車過TS隧道時,y=0.77x2.08,其相關(guān)性系數(shù)為0.98;列車通過TB隧道時,y=0.68x2.04,其相關(guān)性系數(shù)為0.99。

    2.4 隧道內(nèi)附屬物壓力波動分析

    地鐵列車以120 km/h速度等級通過990 m長TB隧道時,隧道內(nèi)不同位置處的壓力波動情況如圖13所示。

    由圖13可知,當(dāng)列車頭部剛進(jìn)入隧道的時候,產(chǎn)生的初始壓縮波到達(dá)測點處,測點壓力值急劇增加,隨著列車逐漸進(jìn)入隧道,壓力也逐漸增加,然后當(dāng)壓縮波到達(dá)隧道出口后,部分以膨脹波的形式反射回來,且列車尾部進(jìn)入隧道時也會產(chǎn)生膨脹波,可使壓力減小;當(dāng)列車頭部經(jīng)過測點時,壓力會上升直到列車尾部經(jīng)過,之后測點的壓力就會隨著壓力波在隧道內(nèi)的傳播而不斷變化。最大正壓力為1 000 Pa,出現(xiàn)在入口處;最大負(fù)壓為700 Pa,出現(xiàn)在中間測點處。

    以隧道中間處壁面監(jiān)測點(in-496)的壓力峰峰值為例,做出速度與壓力峰峰值的關(guān)系圖,如圖14所示。由圖14可知,隧道壁面壓力峰峰值近似與地鐵列車速度平方呈正比。通過曲線擬合可得:當(dāng)列車通過TS隧道時,y=0.078x2.03,其相關(guān)性系數(shù)為0.99;列車通過TB隧道時;y=0.065x2.02,其相關(guān)性系數(shù)為0.99。

    圖14 列車以不同速度通過隧道時隧道壁面壓力峰峰值

    3 結(jié)論

    本文通過計算流體力學(xué),利用滑移網(wǎng)格技術(shù)對列車在隧道內(nèi)底板壓力的氣動性能進(jìn)行研究,得到如下結(jié)論:

    (1) 該型號地鐵車輛通過隧道時,最大正壓位于頭車鼻尖處,頭車和尾車側(cè)翼處負(fù)壓較大,非流線型部位各測點的壓力峰峰值差別不大。

    (2) 列車在隧道中行駛時,在速度相等的條件下,隧道斷面凈空面積越小,車體承受的壓力正峰值、負(fù)峰值以及峰峰值越大。

    (3) 為達(dá)到人體舒適性評價指標(biāo),該型號地鐵車輛在通過大斷面隧道和小斷面隧道時,建議車速應(yīng)分別滿足100 km/h和90 km/h的限速。

    (4) 地鐵列車通過隧道時,車體表面壓力峰峰值、3 s內(nèi)車內(nèi)壓力波動最大值、隧道內(nèi)附屬物壓力峰峰值與列車速度的平方近似成線性關(guān)系。

    [1] HOWE M S.Pressure transients generated when high-speed trains pass in a tunnel [J].IMA Journal of Applied Mathematics,2000,65(3):315.

    [2] HOWE M S.Influence of train Mach number on the compression wave generated in a tunnel-entrance hood [J].Journal of Engineering Mathematics,2003,46(2):147.

    [3] ATTURO B,MICHELE M,STEFANO S. The alleviation of the aerodynamic drag and wave effects of high-speed trains in very long tunnels[J].Journal of Wind Engineering and Industrial Aerodynamics,2001(89):365.

    [4] 軌道交通安全教育部重點實驗室.遂渝線200 km/h隧道空氣動力學(xué)實車試驗研究報告[R].長沙:中南大學(xué),2005.

    [5] 軌道交通安全教育部重點實驗室.遂渝線200 km/h無碴軌道隧道空氣動力學(xué)實車試驗研究報告[R].長沙:中南大學(xué),2006.

    [6] 王英學(xué),高波,駱建軍,等. 高速列車進(jìn)入隧道空氣動力學(xué)模型實驗分析[J].空氣動力學(xué),2004,22(3):346.

    [7] 余南陽,雷波,許志浩,等.高速列車隧道壓力波淺水槽模擬試驗研究[J].鐵道學(xué)報,2004,26(2):52.

    [8] 田紅旗.列車空氣動力學(xué)[M].北京:中國鐵道出版社,2004.

    [9] 王瀟芹,梅元貴,周朝軍. CRH3高速列車隧道壓力波特性淺析[J].蘭州交通大學(xué)學(xué)報,2011,30(1):117.

    [10] 郭曉穎.CRH5動車組隧道運行壓力波數(shù)值研究[D].大連:大連交通大學(xué),2013.

    [11] 馬偉斌,張千里,程愛君,等. 高速鐵路隧道洞口微氣壓波影響因素與變化規(guī)律研究[J].鐵道學(xué)報,2013,35(5):97.

    [12] 劉鳳華,余以正.地鐵列車隧道氣動力學(xué)試驗與仿真[J].大連交通大學(xué)學(xué)報,2013,34(4):7.

    [13] 林世生.關(guān)于廣州地鐵隧道空氣動力學(xué)效應(yīng)緩解措施的研究[J].隧道建設(shè),2013,33(8):650.

    [14] 吳煒,彭金龍.快速地鐵隧道空氣動力學(xué)效應(yīng)[J].城市軌道交通研究,2011,(12):37.

    [15] GAYLARD A P.The application of computational fluid dynamics to railway acrodynamics [J].Journal of Rail and Rapid Transit,Proc Instn Mech Engrs,1993,207(7):133.

    [16] LI T,ZHANG J Y,ZHANG W H. An improved algorithm for fluid-structure interaction of high-speed trains under crosswind[J].Journal of Modern Transportation,2011,19(2):75.

    [17] 王建宇,萬曉燕,吳劍.隧道長度對瞬變壓力的影響[J].現(xiàn)代隧道技術(shù),2008,45(6):1.

    [18] 王建宇,萬曉燕,吳劍.高速鐵路隧道內(nèi)瞬變氣壓和乘車舒適度準(zhǔn)則[J].現(xiàn)代隧道技術(shù),2008,45(2):1.

    Aerodynamic Performance of Metro Vehicle Passing Through Tunnels

    XU Shinan, ZHANG Jiye, XIONG Jun, MENG Tian

    The aerodynamic effect induced by a train when entering a tunnel has negative influence on the train running safety and passengers' comfort. Based on aerodynamics, the pressure temporal curves and pressure distribution characteristics of metro vehicle and long tunnels are analyzed according to the numerical simulation with software FLUENT. The results show that the amplitude of pressure on the carbody, the maximum of internal pressure fluctuation during three seconds, the peak-to-peak pressure value of ancillary facilities in tunnels are approximately proportional to the square of car speed. In addition, the smaller the tunnel cross-section area is, the lager pressure that the carbody has to bear. In order to satisfy the standard requirements of human body comfort, the train speed should be constrainted when the metro vehicle is running in tunnels.

    metro train; tunnel; aerodynamic effect; pressure temporal curve

    *高速鐵路基礎(chǔ)研究聯(lián)合基金資助項目(U1234208);國家自然科學(xué)基金面上項目(51475394);牽引動力國家重點實驗室自主課題(2016TPL_T02)

    V 211

    10.16037/j.1007-869x.2016.09.021

    2014-11-07)

    猜你喜歡
    尾車頭車峰峰
    混勻堆料機(jī)尾車調(diào)偏研究
    單堆單取斗輪機(jī)尾車防脫鉤保護(hù)裝置及連鎖
    丁輝:閱兵坦克“頭車”駕駛員
    晚晴(2020年8期)2020-12-03 13:49:57
    頭車半自動鉤緩裝置傾斜問題研究
    斗輪堆取料機(jī)尾車改造研究
    愛遲到
    年輕真好
    Numerical simulation of Gurney flap on SFYT15thick airfoil
    幼年綠
    白字小先生
    成人三级黄色视频| 国产精品精品国产色婷婷| 亚洲精品影视一区二区三区av| 成年免费大片在线观看| 99热精品在线国产| 中文亚洲av片在线观看爽| 最新中文字幕久久久久| 国产一区二区三区在线臀色熟女| 久久九九热精品免费| 国产精品伦人一区二区| 亚洲国产日韩欧美精品在线观看| 亚洲电影在线观看av| 黄色女人牲交| 91久久精品国产一区二区成人| 在线免费观看不下载黄p国产 | 老熟妇仑乱视频hdxx| 99久久中文字幕三级久久日本| 校园人妻丝袜中文字幕| 色视频www国产| 精品久久久久久久人妻蜜臀av| 日本-黄色视频高清免费观看| 搡女人真爽免费视频火全软件 | 中亚洲国语对白在线视频| 成人二区视频| av在线亚洲专区| ponron亚洲| 日本精品一区二区三区蜜桃| 91在线精品国自产拍蜜月| xxxwww97欧美| 看片在线看免费视频| 中文字幕av在线有码专区| 别揉我奶头 嗯啊视频| 两人在一起打扑克的视频| 午夜爱爱视频在线播放| 干丝袜人妻中文字幕| 精品无人区乱码1区二区| 亚洲四区av| 亚洲图色成人| 国语自产精品视频在线第100页| 成人三级黄色视频| 男女之事视频高清在线观看| 啦啦啦啦在线视频资源| 成人永久免费在线观看视频| 亚洲乱码一区二区免费版| 色综合站精品国产| 人妻少妇偷人精品九色| 欧美潮喷喷水| 久久久久久九九精品二区国产| 国产高清视频在线观看网站| 给我免费播放毛片高清在线观看| 日韩欧美国产在线观看| 国产熟女欧美一区二区| 久久精品国产清高在天天线| 日韩高清综合在线| 老师上课跳d突然被开到最大视频| 一个人看的www免费观看视频| 日韩中字成人| 成人午夜高清在线视频| 婷婷丁香在线五月| 欧美激情国产日韩精品一区| 亚洲国产精品sss在线观看| 欧美性猛交╳xxx乱大交人| 非洲黑人性xxxx精品又粗又长| 成人欧美大片| 91午夜精品亚洲一区二区三区 | 伦精品一区二区三区| 日本五十路高清| 嫩草影院精品99| 国产中年淑女户外野战色| 国产激情偷乱视频一区二区| 大型黄色视频在线免费观看| 亚洲无线在线观看| 麻豆精品久久久久久蜜桃| 在线国产一区二区在线| av中文乱码字幕在线| 亚洲最大成人av| aaaaa片日本免费| 可以在线观看毛片的网站| 国产免费一级a男人的天堂| 国产精品一区二区三区四区久久| 国产精品一区www在线观看 | 欧美成人性av电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩东京热| 深夜a级毛片| 听说在线观看完整版免费高清| 中文字幕精品亚洲无线码一区| 免费一级毛片在线播放高清视频| 亚洲精品亚洲一区二区| 我要看日韩黄色一级片| 夜夜夜夜夜久久久久| 中文资源天堂在线| 亚洲国产精品sss在线观看| 亚洲七黄色美女视频| 人妻久久中文字幕网| 亚洲av五月六月丁香网| 日韩精品有码人妻一区| 精品久久久久久久久av| 久久久久国内视频| 国内精品久久久久精免费| 亚洲电影在线观看av| 亚洲av成人av| 狂野欧美白嫩少妇大欣赏| 可以在线观看毛片的网站| 国产亚洲精品久久久com| 在线播放国产精品三级| 日本-黄色视频高清免费观看| 亚洲av成人av| 一级黄片播放器| av在线亚洲专区| av国产免费在线观看| 色播亚洲综合网| 日本在线视频免费播放| 成人特级黄色片久久久久久久| 日本精品一区二区三区蜜桃| 18禁黄网站禁片午夜丰满| 久久久精品大字幕| 国产老妇女一区| 日本 av在线| 色综合站精品国产| 美女高潮的动态| 伊人久久精品亚洲午夜| 久久精品国产亚洲av香蕉五月| 久久天躁狠狠躁夜夜2o2o| 国产综合懂色| 成年女人看的毛片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 色视频www国产| 黄色一级大片看看| 丰满人妻一区二区三区视频av| 99久久精品国产国产毛片| 99精品在免费线老司机午夜| 亚洲真实伦在线观看| а√天堂www在线а√下载| 亚洲精品乱码久久久v下载方式| 99视频精品全部免费 在线| 亚洲精品色激情综合| 午夜爱爱视频在线播放| 亚洲,欧美,日韩| 国产不卡一卡二| 成人无遮挡网站| 精品久久久噜噜| 国产极品精品免费视频能看的| 久久中文看片网| 18禁裸乳无遮挡免费网站照片| 国产aⅴ精品一区二区三区波| 身体一侧抽搐| 俄罗斯特黄特色一大片| 亚洲不卡免费看| 成年女人看的毛片在线观看| 亚洲最大成人av| 亚洲av不卡在线观看| 国产精品一区二区三区四区免费观看 | 国产爱豆传媒在线观看| 天美传媒精品一区二区| 亚洲精品乱码久久久v下载方式| 亚洲国产高清在线一区二区三| 国产黄a三级三级三级人| 99热这里只有是精品在线观看| 久久久久国产精品人妻aⅴ院| 久久久久久伊人网av| 麻豆国产av国片精品| 在线播放无遮挡| 久久精品国产自在天天线| 男人的好看免费观看在线视频| 国产精品野战在线观看| 国产女主播在线喷水免费视频网站 | 国产久久久一区二区三区| 亚洲不卡免费看| 国产精品不卡视频一区二区| 在线观看66精品国产| 小说图片视频综合网站| 嫩草影视91久久| 欧美一区二区亚洲| 欧美色欧美亚洲另类二区| 亚洲精品影视一区二区三区av| 简卡轻食公司| 精品99又大又爽又粗少妇毛片 | 欧美日韩国产亚洲二区| 亚洲av二区三区四区| 女同久久另类99精品国产91| 亚洲成人中文字幕在线播放| 国产一区二区在线av高清观看| 91在线精品国自产拍蜜月| 日本欧美国产在线视频| 天天躁日日操中文字幕| 日韩亚洲欧美综合| 亚洲天堂国产精品一区在线| 欧美不卡视频在线免费观看| 日本与韩国留学比较| 国产黄a三级三级三级人| 99久久精品国产国产毛片| 又爽又黄a免费视频| 一个人看的www免费观看视频| 国产高清视频在线播放一区| 国产精品一区二区三区四区久久| 日韩欧美精品免费久久| 内地一区二区视频在线| 亚洲经典国产精华液单| 成人三级黄色视频| 男人和女人高潮做爰伦理| 国内精品久久久久精免费| 一进一出抽搐gif免费好疼| 亚洲av.av天堂| www.www免费av| 亚洲第一电影网av| 日本成人三级电影网站| www.色视频.com| 日韩亚洲欧美综合| 亚洲成av人片在线播放无| 成人三级黄色视频| 一区二区三区四区激情视频 | 国产三级中文精品| 九九久久精品国产亚洲av麻豆| 九色成人免费人妻av| 久9热在线精品视频| 99精品在免费线老司机午夜| 国内精品美女久久久久久| 深夜精品福利| 国产精品久久视频播放| 欧洲精品卡2卡3卡4卡5卡区| 少妇高潮的动态图| 国产蜜桃级精品一区二区三区| 夜夜爽天天搞| 色尼玛亚洲综合影院| 99久久九九国产精品国产免费| 无遮挡黄片免费观看| av在线天堂中文字幕| 婷婷丁香在线五月| 国内精品美女久久久久久| 男女之事视频高清在线观看| 久久久久久久亚洲中文字幕| 日韩欧美精品v在线| 免费观看的影片在线观看| 91狼人影院| 色吧在线观看| 精品日产1卡2卡| 日本黄大片高清| 国产黄色小视频在线观看| 成人综合一区亚洲| 久久这里只有精品中国| 韩国av在线不卡| 九九在线视频观看精品| 免费观看人在逋| 国内精品久久久久精免费| 精华霜和精华液先用哪个| 亚洲国产色片| 内射极品少妇av片p| 精品午夜福利在线看| 国产精品久久久久久久久免| 欧美+亚洲+日韩+国产| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| 精品久久久久久,| av视频在线观看入口| 十八禁国产超污无遮挡网站| 国内少妇人妻偷人精品xxx网站| 亚洲狠狠婷婷综合久久图片| 午夜福利成人在线免费观看| 在线天堂最新版资源| 长腿黑丝高跟| 如何舔出高潮| 日韩中文字幕欧美一区二区| 国产国拍精品亚洲av在线观看| 久久精品综合一区二区三区| 成人av在线播放网站| 少妇的逼水好多| 91午夜精品亚洲一区二区三区 | 99热这里只有是精品50| 两个人的视频大全免费| 精品久久久久久久久av| 国内精品久久久久久久电影| 欧美人与善性xxx| 精品久久久久久久久亚洲 | 免费在线观看日本一区| 国产乱人视频| 一个人免费在线观看电影| 国产久久久一区二区三区| 日本黄色视频三级网站网址| 麻豆国产97在线/欧美| 久久精品国产亚洲网站| 免费看光身美女| 久久久成人免费电影| 简卡轻食公司| 在线观看免费视频日本深夜| 中文字幕熟女人妻在线| 直男gayav资源| 欧美日本亚洲视频在线播放| 欧美xxxx性猛交bbbb| 在线天堂最新版资源| 久久6这里有精品| 波多野结衣巨乳人妻| 久久精品国产清高在天天线| 一级黄片播放器| 老司机深夜福利视频在线观看| 色哟哟·www| 观看免费一级毛片| 色综合亚洲欧美另类图片| 国产v大片淫在线免费观看| 久久热精品热| 别揉我奶头~嗯~啊~动态视频| 特级一级黄色大片| 欧美性猛交╳xxx乱大交人| 99精品在免费线老司机午夜| 深夜精品福利| 久久人妻av系列| 热99在线观看视频| 成熟少妇高潮喷水视频| 一级a爱片免费观看的视频| 精品无人区乱码1区二区| 精品一区二区三区人妻视频| 性色avwww在线观看| 欧美激情久久久久久爽电影| 成年女人永久免费观看视频| 亚洲av成人精品一区久久| 国产不卡一卡二| 成人鲁丝片一二三区免费| 亚洲专区中文字幕在线| 午夜爱爱视频在线播放| 日本-黄色视频高清免费观看| 欧美成人a在线观看| 成人精品一区二区免费| 麻豆国产97在线/欧美| 国产精品一区二区免费欧美| 亚洲综合色惰| 99久久成人亚洲精品观看| 五月伊人婷婷丁香| 成人二区视频| 国产三级中文精品| 少妇人妻精品综合一区二区 | 精品99又大又爽又粗少妇毛片 | 日本一本二区三区精品| 精品乱码久久久久久99久播| 国产大屁股一区二区在线视频| 国产一区二区三区av在线 | 色哟哟·www| 天天躁日日操中文字幕| 国产精品伦人一区二区| 亚洲中文字幕日韩| 久久人人爽人人爽人人片va| 非洲黑人性xxxx精品又粗又长| 婷婷亚洲欧美| 久久人人爽人人爽人人片va| 亚洲一区高清亚洲精品| 一个人看视频在线观看www免费| 久久精品人妻少妇| а√天堂www在线а√下载| 99久国产av精品| 国产精品久久久久久久电影| 国产精品,欧美在线| 国产日本99.免费观看| 天天一区二区日本电影三级| 国产精品国产三级国产av玫瑰| 国产精品一区www在线观看 | 我要搜黄色片| 十八禁国产超污无遮挡网站| 最好的美女福利视频网| 高清毛片免费观看视频网站| 在线观看舔阴道视频| 国产精品久久久久久精品电影| 给我免费播放毛片高清在线观看| 91狼人影院| 久久中文看片网| 免费看av在线观看网站| 国产一级毛片七仙女欲春2| 性色avwww在线观看| 一级黄色大片毛片| 美女cb高潮喷水在线观看| 久久香蕉精品热| 99riav亚洲国产免费| av福利片在线观看| 欧美成人一区二区免费高清观看| 中文资源天堂在线| 一个人免费在线观看电影| 亚洲,欧美,日韩| 我要搜黄色片| 成人国产综合亚洲| 欧美xxxx黑人xx丫x性爽| 在线观看66精品国产| 精品久久久久久久末码| .国产精品久久| 亚洲专区国产一区二区| 中国美女看黄片| 女同久久另类99精品国产91| 亚洲四区av| 国产一区二区三区视频了| 亚洲最大成人手机在线| av在线蜜桃| 久9热在线精品视频| 国产老妇女一区| 嫁个100分男人电影在线观看| 亚洲成人免费电影在线观看| 精品一区二区三区视频在线| 日韩精品青青久久久久久| 99热6这里只有精品| 嫩草影视91久久| aaaaa片日本免费| 久久人人精品亚洲av| 国产视频内射| 男女啪啪激烈高潮av片| 桃色一区二区三区在线观看| 3wmmmm亚洲av在线观看| 99久久成人亚洲精品观看| 久久中文看片网| 午夜精品一区二区三区免费看| 天堂√8在线中文| 美女被艹到高潮喷水动态| 丰满乱子伦码专区| 国产成人一区二区在线| 免费人成视频x8x8入口观看| 全区人妻精品视频| 麻豆国产av国片精品| 久久热精品热| 麻豆精品久久久久久蜜桃| 精品久久久久久,| 成熟少妇高潮喷水视频| 国产私拍福利视频在线观看| 免费人成在线观看视频色| 国产一区二区激情短视频| 欧美高清成人免费视频www| 国产精品爽爽va在线观看网站| 中国美白少妇内射xxxbb| 超碰av人人做人人爽久久| 日韩,欧美,国产一区二区三区 | 一区二区三区激情视频| 精品人妻熟女av久视频| 午夜影院日韩av| 老司机福利观看| 久久99热这里只有精品18| 国产大屁股一区二区在线视频| 国产高清视频在线观看网站| 午夜免费男女啪啪视频观看 | 久久久国产成人精品二区| av国产免费在线观看| 免费黄网站久久成人精品| 免费av毛片视频| 国产精品一区二区免费欧美| a级毛片a级免费在线| 亚洲性久久影院| 日本爱情动作片www.在线观看 | 精品国产三级普通话版| 国产精品一区二区免费欧美| 淫秽高清视频在线观看| 少妇熟女aⅴ在线视频| 精品久久久久久久末码| 悠悠久久av| 韩国av在线不卡| 亚洲欧美日韩无卡精品| 又爽又黄a免费视频| 欧美精品国产亚洲| 国产色爽女视频免费观看| 少妇丰满av| 啦啦啦啦在线视频资源| 小蜜桃在线观看免费完整版高清| 欧美日韩中文字幕国产精品一区二区三区| 老司机深夜福利视频在线观看| 亚洲精品影视一区二区三区av| 成年女人看的毛片在线观看| 精品无人区乱码1区二区| 97人妻精品一区二区三区麻豆| 亚洲国产精品久久男人天堂| 国产高清视频在线观看网站| 国产亚洲精品av在线| 成人午夜高清在线视频| 国产午夜福利久久久久久| 午夜激情欧美在线| 欧美bdsm另类| 99热这里只有是精品50| 听说在线观看完整版免费高清| 国产毛片a区久久久久| 亚洲欧美清纯卡通| 91午夜精品亚洲一区二区三区 | 亚洲一级一片aⅴ在线观看| 联通29元200g的流量卡| 日韩欧美免费精品| 在线观看66精品国产| 国产淫片久久久久久久久| 2021天堂中文幕一二区在线观| 欧美极品一区二区三区四区| 欧美黑人欧美精品刺激| 中文资源天堂在线| 欧美成人一区二区免费高清观看| 日本免费一区二区三区高清不卡| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久精品电影| 国产精品电影一区二区三区| 亚洲七黄色美女视频| 日韩中文字幕欧美一区二区| 欧美色欧美亚洲另类二区| 热99在线观看视频| 日韩一区二区视频免费看| 免费看美女性在线毛片视频| 欧美区成人在线视频| 能在线免费观看的黄片| 黄片wwwwww| 精品一区二区三区av网在线观看| av在线天堂中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 久9热在线精品视频| 无遮挡黄片免费观看| 一本精品99久久精品77| 国产一区二区三区在线臀色熟女| 黄色日韩在线| 可以在线观看的亚洲视频| 一级a爱片免费观看的视频| 久久人妻av系列| 最近视频中文字幕2019在线8| 看十八女毛片水多多多| 国产欧美日韩一区二区精品| 自拍偷自拍亚洲精品老妇| 高清在线国产一区| 淫秽高清视频在线观看| 少妇熟女aⅴ在线视频| 亚洲最大成人手机在线| 亚洲自拍偷在线| 国产aⅴ精品一区二区三区波| 亚洲熟妇熟女久久| 国产色婷婷99| АⅤ资源中文在线天堂| 免费人成在线观看视频色| 久久久久久久午夜电影| 黄色欧美视频在线观看| 波多野结衣高清作品| 国产高清三级在线| 亚州av有码| 久久精品国产鲁丝片午夜精品 | 99久国产av精品| 欧美日韩精品成人综合77777| 婷婷丁香在线五月| 人人妻人人澡欧美一区二区| 日本黄色片子视频| 少妇人妻一区二区三区视频| 免费看av在线观看网站| 一级a爱片免费观看的视频| 特大巨黑吊av在线直播| 亚洲成人中文字幕在线播放| 成人亚洲精品av一区二区| 国产主播在线观看一区二区| 亚洲不卡免费看| 高清日韩中文字幕在线| 精品久久久久久久久久免费视频| 五月玫瑰六月丁香| 成人一区二区视频在线观看| 亚洲无线在线观看| 搡老熟女国产l中国老女人| 校园人妻丝袜中文字幕| 国产精品国产高清国产av| videossex国产| 在线天堂最新版资源| 夜夜夜夜夜久久久久| 久久中文看片网| av专区在线播放| 少妇猛男粗大的猛烈进出视频 | 一级黄色大片毛片| 麻豆精品久久久久久蜜桃| 色吧在线观看| 熟妇人妻久久中文字幕3abv| 1024手机看黄色片| www日本黄色视频网| 窝窝影院91人妻| 99久久无色码亚洲精品果冻| 日韩精品青青久久久久久| 在线播放国产精品三级| 极品教师在线免费播放| 日本-黄色视频高清免费观看| 成年女人毛片免费观看观看9| 国产欧美日韩精品亚洲av| 婷婷亚洲欧美| 99热6这里只有精品| 免费电影在线观看免费观看| 老司机深夜福利视频在线观看| 欧美成人a在线观看| 日日啪夜夜撸| 久久人妻av系列| 亚洲成人免费电影在线观看| 欧美3d第一页| 哪里可以看免费的av片| 国产精品国产高清国产av| 三级毛片av免费| 99久久中文字幕三级久久日本| 亚洲,欧美,日韩| 精品一区二区三区视频在线| 久久这里只有精品中国| 搡老熟女国产l中国老女人| 亚洲av免费高清在线观看| 少妇被粗大猛烈的视频| 国产精品野战在线观看| 中国美女看黄片| 亚洲自拍偷在线| 欧美最黄视频在线播放免费| 天天躁日日操中文字幕| 自拍偷自拍亚洲精品老妇| 欧美一级a爱片免费观看看| 99国产精品一区二区蜜桃av| 波多野结衣高清无吗| 男女那种视频在线观看| 亚洲av第一区精品v没综合| 国产爱豆传媒在线观看| 精品久久久久久久久久久久久| 国内毛片毛片毛片毛片毛片| 悠悠久久av| 久久久久精品国产欧美久久久| 亚洲av一区综合| 国产爱豆传媒在线观看| 亚洲天堂国产精品一区在线| а√天堂www在线а√下载| 久久精品人妻少妇| 十八禁国产超污无遮挡网站| 欧美zozozo另类| 日本五十路高清| 免费观看精品视频网站| 欧美成人一区二区免费高清观看| 成人三级黄色视频| 午夜视频国产福利| 国产亚洲av嫩草精品影院|