• <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免费观看日本| 精品国产一区二区三区久久久樱花 | 成年人午夜在线观看视频| 国产精品久久久久久精品电影| 国产精品蜜桃在线观看| 成人免费观看视频高清| 国产精品国产三级国产av玫瑰| 久久久久国产精品人妻一区二区| 欧美激情久久久久久爽电影| 一级毛片我不卡| 精品久久久久久久久亚洲| 一级毛片电影观看| 日本黄大片高清| 中文字幕人妻熟人妻熟丝袜美| 亚洲aⅴ乱码一区二区在线播放| 国产极品天堂在线| 人人妻人人看人人澡| 国产熟女欧美一区二区| 高清毛片免费看| 亚洲国产精品成人综合色| 97人妻精品一区二区三区麻豆| 久久久久久久久久久免费av| 蜜桃亚洲精品一区二区三区| 日本午夜av视频| 网址你懂的国产日韩在线| 女的被弄到高潮叫床怎么办| 天天躁日日操中文字幕| 日韩av免费高清视频| 永久网站在线| 97精品久久久久久久久久精品| 寂寞人妻少妇视频99o| 欧美高清性xxxxhd video| 国产一区亚洲一区在线观看| 一级毛片aaaaaa免费看小| 麻豆成人午夜福利视频| 国产男女内射视频| 国产精品一区二区在线观看99| 亚洲国产精品国产精品| 亚洲精品久久午夜乱码| 深夜a级毛片| 成人欧美大片| 欧美高清性xxxxhd video| 三级国产精品欧美在线观看| 久久精品国产鲁丝片午夜精品| 午夜亚洲福利在线播放| 黄色视频在线播放观看不卡| 国产精品伦人一区二区| 亚洲高清免费不卡视频| 国产综合懂色| 综合色丁香网| 日韩精品有码人妻一区| 亚洲欧美一区二区三区黑人 | 精品少妇久久久久久888优播| 国产黄频视频在线观看| 久久精品国产亚洲av涩爱| 亚洲人成网站在线播| 在线观看免费高清a一片| 麻豆国产97在线/欧美| 亚洲自拍偷在线| 精品国产一区二区三区久久久樱花 | 日韩一区二区三区影片| 视频中文字幕在线观看| 一区二区av电影网| 亚洲激情五月婷婷啪啪| 成人综合一区亚洲| 久久99热6这里只有精品| 精品久久久久久久末码| 少妇人妻 视频| 男的添女的下面高潮视频| 国产老妇伦熟女老妇高清| 久久久久久伊人网av| 人人妻人人澡人人爽人人夜夜| 久久久久九九精品影院| 欧美亚洲 丝袜 人妻 在线| 91狼人影院| 亚洲av在线观看美女高潮| 亚洲国产最新在线播放| 草草在线视频免费看| 寂寞人妻少妇视频99o| 人妻 亚洲 视频| 91精品国产九色| 欧美zozozo另类| 99久久中文字幕三级久久日本| 边亲边吃奶的免费视频| 在线精品无人区一区二区三 | 菩萨蛮人人尽说江南好唐韦庄| 九草在线视频观看| 欧美97在线视频| 伊人久久国产一区二区| 中文乱码字字幕精品一区二区三区| 听说在线观看完整版免费高清| 一区二区av电影网| 免费高清在线观看视频在线观看| 亚洲天堂av无毛| 自拍偷自拍亚洲精品老妇| 成人免费观看视频高清| 夜夜爽夜夜爽视频| 亚洲成人一二三区av| 老司机影院成人| 毛片一级片免费看久久久久| 欧美国产精品一级二级三级 | 夫妻午夜视频| 国产 一区精品| 日日啪夜夜撸| 一级毛片黄色毛片免费观看视频| 日日摸夜夜添夜夜添av毛片| 高清在线视频一区二区三区| 少妇 在线观看| 午夜福利视频精品| 国国产精品蜜臀av免费| 精品国产三级普通话版| 亚洲欧美精品自产自拍| 香蕉精品网在线| 午夜免费鲁丝| 国产欧美另类精品又又久久亚洲欧美| 国产真实伦视频高清在线观看| 国产精品.久久久| 国产成人a区在线观看| 极品教师在线视频| 久久久久网色| 亚洲欧美日韩卡通动漫| 18禁裸乳无遮挡动漫免费视频 | 欧美性感艳星| 国产 一区精品| 亚洲欧美一区二区三区黑人 | 色视频www国产| 有码 亚洲区| 精品一区二区免费观看| 亚洲欧美精品专区久久| 男的添女的下面高潮视频| 精品国产一区二区三区久久久樱花 | 亚洲,欧美,日韩| 久久久久久久久久久丰满| 午夜激情福利司机影院| 中文乱码字字幕精品一区二区三区| 午夜福利视频精品| 久久久久久久精品精品| 久久人人爽人人爽人人片va| 日韩人妻高清精品专区| 99久久九九国产精品国产免费| 成人欧美大片| 国产在线一区二区三区精| av女优亚洲男人天堂| 啦啦啦中文免费视频观看日本| 久久精品综合一区二区三区| 99热网站在线观看| 欧美日韩视频精品一区| 亚洲精品乱码久久久久久按摩| 九九在线视频观看精品| 看十八女毛片水多多多| 国产又色又爽无遮挡免| 一级毛片aaaaaa免费看小| 青春草视频在线免费观看| 午夜精品一区二区三区免费看| 久久精品国产亚洲av天美| 啦啦啦啦在线视频资源| 97精品久久久久久久久久精品| 大片免费播放器 马上看| 久久99热这里只频精品6学生| 男女边吃奶边做爰视频| 欧美精品国产亚洲| 高清日韩中文字幕在线| 美女主播在线视频| 久久久久久久亚洲中文字幕| 在线免费观看不下载黄p国产| 一级毛片我不卡| 毛片一级片免费看久久久久| 久久韩国三级中文字幕| 国产黄片视频在线免费观看| av一本久久久久| 亚洲内射少妇av| 免费不卡的大黄色大毛片视频在线观看| 精华霜和精华液先用哪个| 国产精品偷伦视频观看了| 免费看a级黄色片| 在线观看一区二区三区| 久久久成人免费电影| 看非洲黑人一级黄片| 3wmmmm亚洲av在线观看| 精华霜和精华液先用哪个| av播播在线观看一区| 深夜a级毛片| 国产一区二区三区av在线| 国产亚洲5aaaaa淫片| 国产一区二区三区av在线| 成人黄色视频免费在线看| 国产爱豆传媒在线观看| 边亲边吃奶的免费视频| 久久热精品热| 午夜日本视频在线| 亚洲精品,欧美精品| 国产午夜精品久久久久久一区二区三区| 中文字幕av成人在线电影| 亚洲av欧美aⅴ国产| 久久久亚洲精品成人影院| 午夜福利在线在线| 少妇人妻一区二区三区视频| 精品人妻视频免费看| av国产免费在线观看| 少妇人妻 视频| 国产成人91sexporn| 日韩av免费高清视频| 久久久精品94久久精品| 国产高清国产精品国产三级 | 九九久久精品国产亚洲av麻豆| 精品久久久久久电影网| 少妇人妻久久综合中文| 丰满人妻一区二区三区视频av| 少妇人妻久久综合中文| 又粗又硬又长又爽又黄的视频| av一本久久久久| 嫩草影院精品99| 超碰av人人做人人爽久久| 啦啦啦中文免费视频观看日本| 少妇裸体淫交视频免费看高清| 69人妻影院| av在线天堂中文字幕| 国产综合精华液| 亚洲色图综合在线观看| 毛片一级片免费看久久久久| 一本久久精品| 2021少妇久久久久久久久久久| 春色校园在线视频观看| 男插女下体视频免费在线播放| 日本欧美国产在线视频| 久久精品夜色国产| 久久久久性生活片| 亚洲久久久久久中文字幕| 国产 一区 欧美 日韩| 免费人成在线观看视频色| 欧美老熟妇乱子伦牲交| h日本视频在线播放| 国产免费福利视频在线观看| 国产精品偷伦视频观看了| 丰满少妇做爰视频| 国产精品熟女久久久久浪| 亚洲欧美日韩另类电影网站 | 一区二区三区免费毛片| 国产一区二区在线观看日韩| 色网站视频免费| 91aial.com中文字幕在线观看| 网址你懂的国产日韩在线| 国产熟女欧美一区二区| 国产成人freesex在线| 波野结衣二区三区在线| 熟女av电影| 国产毛片在线视频| 久久亚洲国产成人精品v| 高清视频免费观看一区二区| 少妇高潮的动态图| 日本免费在线观看一区| av.在线天堂| 蜜臀久久99精品久久宅男| 日韩一区二区视频免费看| 最近2019中文字幕mv第一页| 成人二区视频| 成人亚洲精品av一区二区| 99热这里只有精品一区| 国产日韩欧美在线精品| 人妻 亚洲 视频| 久久久久久久久久成人| 99久久九九国产精品国产免费| 美女脱内裤让男人舔精品视频| 亚洲激情五月婷婷啪啪| 国产一区二区在线观看日韩| 爱豆传媒免费全集在线观看| 久久久久久国产a免费观看| 亚洲国产欧美人成| 美女cb高潮喷水在线观看| 亚洲久久久久久中文字幕| 亚洲在线观看片| 三级经典国产精品| 成人毛片60女人毛片免费| 国产亚洲5aaaaa淫片| 最近最新中文字幕免费大全7| 亚洲精品国产成人久久av| 亚洲aⅴ乱码一区二区在线播放| 在现免费观看毛片| 女的被弄到高潮叫床怎么办| 青春草视频在线免费观看| 人人妻人人爽人人添夜夜欢视频 | 国产精品久久久久久久电影| 男女边摸边吃奶| 亚洲欧美精品专区久久| 免费av毛片视频| 18禁在线播放成人免费| 丰满少妇做爰视频| 国产69精品久久久久777片| 丝袜喷水一区| 伦精品一区二区三区| 身体一侧抽搐| 亚洲欧洲国产日韩| 大码成人一级视频| 久久人人爽av亚洲精品天堂 | 日本猛色少妇xxxxx猛交久久| 丝袜美腿在线中文| 熟女人妻精品中文字幕| 日日啪夜夜爽| 久久久精品免费免费高清| 三级男女做爰猛烈吃奶摸视频| 午夜福利在线在线| 一级毛片 在线播放| 久久综合国产亚洲精品| 亚洲精品日本国产第一区| 熟妇人妻不卡中文字幕| 亚洲精华国产精华液的使用体验| 99久久九九国产精品国产免费| 日韩欧美 国产精品| 男女下面进入的视频免费午夜| 久久久精品94久久精品| 亚洲精品国产av成人精品| 亚洲精品456在线播放app| 秋霞伦理黄片| 十八禁网站网址无遮挡 | 午夜激情久久久久久久| 性插视频无遮挡在线免费观看| 麻豆成人午夜福利视频| 国产爽快片一区二区三区| 国产伦理片在线播放av一区| 午夜福利网站1000一区二区三区| 亚洲怡红院男人天堂| 日本爱情动作片www.在线观看| 一本一本综合久久| 欧美成人午夜免费资源| 看黄色毛片网站| 亚洲精品日韩在线中文字幕| 国产精品av视频在线免费观看| 国产成人福利小说| 精品久久国产蜜桃| 性色avwww在线观看| 乱系列少妇在线播放| 高清毛片免费看| 亚洲成人一二三区av| 久久精品熟女亚洲av麻豆精品| 乱码一卡2卡4卡精品| 国产毛片a区久久久久| 2022亚洲国产成人精品| 麻豆国产97在线/欧美| 春色校园在线视频观看| 久久精品国产亚洲av天美| 欧美老熟妇乱子伦牲交| 国内少妇人妻偷人精品xxx网站| 中国三级夫妇交换| 嫩草影院入口| 又爽又黄无遮挡网站| 欧美老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 九色成人免费人妻av| 性色avwww在线观看| 日本三级黄在线观看| 精品酒店卫生间| 国产精品偷伦视频观看了| 精品久久国产蜜桃| 亚洲最大成人手机在线| 亚洲精品一区蜜桃| 久久久久久久久久久免费av| 免费观看a级毛片全部| 听说在线观看完整版免费高清| 秋霞在线观看毛片| 成人一区二区视频在线观看| 日本黄大片高清| 亚洲国产欧美在线一区| 大片免费播放器 马上看| 国产一级毛片在线| 国产爽快片一区二区三区| 久久久欧美国产精品| 22中文网久久字幕| 亚洲欧洲日产国产| 熟女人妻精品中文字幕| 涩涩av久久男人的天堂| 免费av观看视频| 麻豆久久精品国产亚洲av| 国产成人免费观看mmmm| 最后的刺客免费高清国语| 女的被弄到高潮叫床怎么办| 内射极品少妇av片p| 久久热精品热| 精品久久久久久久久av| 欧美日韩国产mv在线观看视频 | 97热精品久久久久久| 可以在线观看毛片的网站| 国产永久视频网站| 爱豆传媒免费全集在线观看| 美女内射精品一级片tv| 香蕉精品网在线| 在线播放无遮挡| 91久久精品电影网| 美女国产视频在线观看| 一级毛片久久久久久久久女| 久久久久性生活片| 欧美变态另类bdsm刘玥| 天堂中文最新版在线下载 | 色吧在线观看| 校园人妻丝袜中文字幕| 韩国高清视频一区二区三区| 18+在线观看网站| 交换朋友夫妻互换小说| 国国产精品蜜臀av免费| 色婷婷久久久亚洲欧美| 成人综合一区亚洲| 中文天堂在线官网| 国产一区二区亚洲精品在线观看| 国产片特级美女逼逼视频| 欧美潮喷喷水| 亚洲av中文av极速乱| 久久综合国产亚洲精品| 免费看日本二区| 黄色一级大片看看| 国产成人91sexporn| 日本一二三区视频观看| 欧美97在线视频| 午夜福利高清视频| a级毛片免费高清观看在线播放| 99热全是精品| 久久鲁丝午夜福利片| 男女下面进入的视频免费午夜| 精品久久久久久久末码| av在线app专区| 天堂中文最新版在线下载 | videossex国产| 真实男女啪啪啪动态图| av在线蜜桃| 欧美精品国产亚洲| 在线免费观看不下载黄p国产| 尾随美女入室| 少妇人妻一区二区三区视频| 国产黄片视频在线免费观看| 九九爱精品视频在线观看| 天堂中文最新版在线下载 | 观看美女的网站| 免费黄网站久久成人精品| 人人妻人人澡人人爽人人夜夜| 大码成人一级视频| 欧美丝袜亚洲另类| 午夜免费观看性视频| 国产女主播在线喷水免费视频网站| 日本猛色少妇xxxxx猛交久久| 午夜激情福利司机影院| 日韩成人伦理影院| 亚洲精品成人久久久久久| 天堂俺去俺来也www色官网| 久久6这里有精品| 欧美高清性xxxxhd video| 在线天堂最新版资源| 大片电影免费在线观看免费| 狂野欧美白嫩少妇大欣赏| 内射极品少妇av片p| 国产精品秋霞免费鲁丝片| 18禁在线无遮挡免费观看视频| 久久99热6这里只有精品| av在线亚洲专区| 一区二区av电影网| 国产精品熟女久久久久浪| a级一级毛片免费在线观看| 久久热精品热| 综合色丁香网| av线在线观看网站| 亚洲欧美一区二区三区黑人 | 搡老乐熟女国产| 一级av片app| 99视频精品全部免费 在线| 亚洲自偷自拍三级| a级毛色黄片| 亚洲av二区三区四区| 亚洲国产欧美人成| 99热这里只有是精品在线观看| 国产精品99久久久久久久久| 亚洲怡红院男人天堂| 大香蕉97超碰在线| 欧美人与善性xxx| 日韩精品有码人妻一区| 性色avwww在线观看| 波多野结衣巨乳人妻| 亚洲欧洲日产国产| 午夜激情福利司机影院| av.在线天堂| 国产精品久久久久久av不卡| 久久鲁丝午夜福利片| 国产男人的电影天堂91| 熟妇人妻不卡中文字幕| 久久精品人妻少妇| 国产亚洲av嫩草精品影院| 99re6热这里在线精品视频| 成人国产av品久久久| 三级国产精品欧美在线观看| 久久ye,这里只有精品| 真实男女啪啪啪动态图| av免费观看日本| 简卡轻食公司| 免费观看av网站的网址| 久久6这里有精品| 不卡视频在线观看欧美| 欧美成人一区二区免费高清观看| 22中文网久久字幕| 又爽又黄无遮挡网站| 日本wwww免费看| 精品视频人人做人人爽| 深爱激情五月婷婷| 成人免费观看视频高清| 国产精品熟女久久久久浪| 亚洲av中文av极速乱| 极品教师在线视频| 国产精品人妻久久久影院| 久热这里只有精品99| 欧美性猛交╳xxx乱大交人| 国产精品.久久久| 色婷婷久久久亚洲欧美| 日本与韩国留学比较| 免费黄网站久久成人精品| 久久人人爽人人爽人人片va| 国产伦在线观看视频一区| 97热精品久久久久久| 欧美精品一区二区大全| 亚洲精品456在线播放app| 丰满人妻一区二区三区视频av| 欧美潮喷喷水| 69人妻影院| 永久网站在线| 久久精品国产亚洲网站| 男女下面进入的视频免费午夜| 黄色欧美视频在线观看| 欧美一区二区亚洲| 欧美日韩综合久久久久久| 最近的中文字幕免费完整| 美女主播在线视频| 大又大粗又爽又黄少妇毛片口| 噜噜噜噜噜久久久久久91| 成人免费观看视频高清| 水蜜桃什么品种好| 欧美日韩精品成人综合77777| 精品人妻一区二区三区麻豆| 激情 狠狠 欧美| 一个人看的www免费观看视频| 精品一区二区三区视频在线| 精品久久久久久久久亚洲| 国产精品麻豆人妻色哟哟久久| 青春草亚洲视频在线观看| 国产精品一区二区在线观看99| 麻豆成人午夜福利视频| 97人妻精品一区二区三区麻豆| 国产大屁股一区二区在线视频| 综合色av麻豆| 三级男女做爰猛烈吃奶摸视频| 精品人妻熟女av久视频| 男女下面进入的视频免费午夜| 欧美97在线视频| 校园人妻丝袜中文字幕| 三级国产精品欧美在线观看| 国产精品伦人一区二区| 国产一区亚洲一区在线观看| 亚洲性久久影院| 18禁裸乳无遮挡免费网站照片| 2022亚洲国产成人精品| 精品亚洲乱码少妇综合久久| 欧美另类一区| 三级经典国产精品| 国产黄色视频一区二区在线观看| 亚洲精品成人久久久久久| 精品视频人人做人人爽| 久久97久久精品| 午夜福利网站1000一区二区三区| 国产精品蜜桃在线观看| 国产av不卡久久| 国产欧美亚洲国产| 麻豆成人午夜福利视频| 亚洲欧美成人综合另类久久久| 可以在线观看毛片的网站| 校园人妻丝袜中文字幕| 国产女主播在线喷水免费视频网站| 亚洲精品456在线播放app| av又黄又爽大尺度在线免费看| 欧美人与善性xxx| 大香蕉97超碰在线| 日韩中字成人| av专区在线播放| 卡戴珊不雅视频在线播放| 亚洲精品亚洲一区二区| 久久精品国产a三级三级三级| 亚洲最大成人av| 99精国产麻豆久久婷婷| a级毛色黄片| 国产欧美另类精品又又久久亚洲欧美| 在线精品无人区一区二区三 | 久久久久久久亚洲中文字幕| 偷拍熟女少妇极品色| 国产av码专区亚洲av| 性色av一级| 国产一区亚洲一区在线观看| 不卡视频在线观看欧美| 水蜜桃什么品种好| 国产精品av视频在线免费观看| 人人妻人人看人人澡| 男插女下体视频免费在线播放| 91aial.com中文字幕在线观看| 国产亚洲av片在线观看秒播厂| 嫩草影院新地址| 国产免费福利视频在线观看| 成年av动漫网址| 精品久久国产蜜桃| 久久女婷五月综合色啪小说 | 人妻夜夜爽99麻豆av| 中文字幕制服av| 久久久精品94久久精品| 欧美最新免费一区二区三区| 青春草国产在线视频| 亚洲自偷自拍三级| 九草在线视频观看| 午夜激情福利司机影院| 色吧在线观看| 国产伦理片在线播放av一区| 久久久色成人| 五月伊人婷婷丁香| 免费看不卡的av| 最近手机中文字幕大全| 日本色播在线视频|