• <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
    幼年綠
    白字小先生
    国产精品国产三级专区第一集| 成年美女黄网站色视频大全免费 | 久久 成人 亚洲| 久久久久久久亚洲中文字幕| 国产亚洲欧美精品永久| 国语对白做爰xxxⅹ性视频网站| 草草在线视频免费看| 亚洲国产欧美日韩在线播放| 国产精品女同一区二区软件| 亚洲色图 男人天堂 中文字幕 | 久久韩国三级中文字幕| 爱豆传媒免费全集在线观看| 亚洲精品日韩在线中文字幕| 国产在线免费精品| 黄色怎么调成土黄色| 久久av网站| 久久人人爽人人片av| 午夜91福利影院| 少妇的逼好多水| 国产黄频视频在线观看| 人人妻人人澡人人爽人人夜夜| 一区二区日韩欧美中文字幕 | 久久久国产一区二区| 免费人妻精品一区二区三区视频| 免费少妇av软件| 丰满饥渴人妻一区二区三| 91久久精品国产一区二区成人| 夫妻午夜视频| 国产 一区精品| 国产男女超爽视频在线观看| 日韩精品免费视频一区二区三区 | 欧美精品一区二区大全| 日韩成人伦理影院| 18禁动态无遮挡网站| 欧美日韩av久久| 午夜免费鲁丝| 国产一区二区在线观看av| 肉色欧美久久久久久久蜜桃| 国产精品不卡视频一区二区| 99热这里只有精品一区| 亚洲四区av| 人人澡人人妻人| 大片电影免费在线观看免费| 91aial.com中文字幕在线观看| 波野结衣二区三区在线| 国产成人午夜福利电影在线观看| 亚洲欧美色中文字幕在线| 国产精品久久久久久精品古装| 最近手机中文字幕大全| √禁漫天堂资源中文www| 国产成人免费观看mmmm| 永久免费av网站大全| 国产在线免费精品| 久久久欧美国产精品| 最后的刺客免费高清国语| 亚洲欧美一区二区三区黑人 | 久久久久久久久大av| 国产极品天堂在线| 少妇高潮的动态图| 男人添女人高潮全过程视频| 天堂俺去俺来也www色官网| 午夜激情av网站| 久久韩国三级中文字幕| 搡老乐熟女国产| 亚洲美女视频黄频| av免费观看日本| 精品熟女少妇av免费看| 日韩在线高清观看一区二区三区| 热re99久久精品国产66热6| 亚洲欧美清纯卡通| 免费大片18禁| 在线观看三级黄色| 18在线观看网站| 亚洲欧美精品自产自拍| 国产伦理片在线播放av一区| 国产高清三级在线| 另类精品久久| 2021少妇久久久久久久久久久| 成人影院久久| 国产伦精品一区二区三区视频9| 少妇人妻久久综合中文| 我的老师免费观看完整版| 97超碰精品成人国产| 97超碰精品成人国产| 黄色欧美视频在线观看| 国产精品蜜桃在线观看| 熟女人妻精品中文字幕| 国模一区二区三区四区视频| 久久99热这里只频精品6学生| 精品久久久久久久久av| 亚洲天堂av无毛| 久久久久久久久久久久大奶| 精品国产一区二区三区久久久樱花| 久久鲁丝午夜福利片| 国产伦理片在线播放av一区| 精品视频人人做人人爽| 国产乱人偷精品视频| 国产精品一二三区在线看| 国产精品嫩草影院av在线观看| 国产在视频线精品| 中文字幕av电影在线播放| 日本av手机在线免费观看| 中文字幕最新亚洲高清| 十分钟在线观看高清视频www| 国产极品粉嫩免费观看在线 | 免费av中文字幕在线| 国产精品国产三级国产av玫瑰| 99热国产这里只有精品6| 欧美日韩在线观看h| 两个人的视频大全免费| 黄色配什么色好看| 91精品一卡2卡3卡4卡| 女的被弄到高潮叫床怎么办| 最近的中文字幕免费完整| 曰老女人黄片| 国产精品熟女久久久久浪| 春色校园在线视频观看| 国产成人a∨麻豆精品| 天堂8中文在线网| 亚洲av电影在线观看一区二区三区| 亚洲内射少妇av| 一级片'在线观看视频| 精品久久国产蜜桃| 国产国拍精品亚洲av在线观看| 亚洲精品色激情综合| 丰满迷人的少妇在线观看| 久久久国产一区二区| 各种免费的搞黄视频| 亚洲综合色网址| 成人国产麻豆网| 免费高清在线观看视频在线观看| 伊人亚洲综合成人网| 国产精品一区二区在线观看99| 不卡视频在线观看欧美| 国产又色又爽无遮挡免| 日韩精品免费视频一区二区三区 | 人人妻人人添人人爽欧美一区卜| 一级a做视频免费观看| 精品国产乱码久久久久久小说| 亚洲国产精品一区二区三区在线| 国产在线视频一区二区| 男女无遮挡免费网站观看| 久久久久久久国产电影| 美女视频免费永久观看网站| 少妇精品久久久久久久| 国产精品久久久久久精品古装| 欧美 亚洲 国产 日韩一| av卡一久久| 色视频在线一区二区三区| 亚洲av成人精品一二三区| 一级,二级,三级黄色视频| 中文字幕人妻丝袜制服| www.av在线官网国产| 亚洲丝袜综合中文字幕| 色婷婷久久久亚洲欧美| 午夜福利视频精品| 免费久久久久久久精品成人欧美视频 | 国产男女内射视频| 亚洲精品乱久久久久久| 国产伦精品一区二区三区视频9| 2022亚洲国产成人精品| 两个人的视频大全免费| 99精国产麻豆久久婷婷| 在线播放无遮挡| 精品亚洲成a人片在线观看| 亚洲精品国产av蜜桃| 一级黄片播放器| 九色亚洲精品在线播放| 熟女人妻精品中文字幕| 久久久久精品性色| 啦啦啦视频在线资源免费观看| 国产片内射在线| 18禁裸乳无遮挡动漫免费视频| 国产精品国产三级专区第一集| 亚洲精品视频女| 大话2 男鬼变身卡| 日韩人妻高清精品专区| 纵有疾风起免费观看全集完整版| 日本免费在线观看一区| 国产成人精品在线电影| 中文字幕人妻丝袜制服| 欧美三级亚洲精品| 欧美激情 高清一区二区三区| 一区二区日韩欧美中文字幕 | 日本av手机在线免费观看| 精品少妇久久久久久888优播| 亚洲色图综合在线观看| 国产成人a∨麻豆精品| 日本欧美国产在线视频| 在线 av 中文字幕| 成人毛片a级毛片在线播放| 亚洲精品第二区| 亚洲人与动物交配视频| 国产精品99久久99久久久不卡 | 日韩 亚洲 欧美在线| 少妇熟女欧美另类| 国产精品久久久久久精品古装| 丝袜在线中文字幕| 纯流量卡能插随身wifi吗| 又黄又爽又刺激的免费视频.| 在线观看一区二区三区激情| 日韩精品免费视频一区二区三区 | 老司机影院成人| 国产 精品1| 色婷婷久久久亚洲欧美| 丝袜脚勾引网站| 一级黄片播放器| 国产欧美日韩一区二区三区在线 | 欧美日韩在线观看h| 精品国产一区二区久久| 五月玫瑰六月丁香| 国产精品一区二区在线不卡| 日本av手机在线免费观看| 99九九在线精品视频| 最近的中文字幕免费完整| 精品人妻熟女av久视频| 国产 精品1| 99热全是精品| 亚洲av不卡在线观看| 99热这里只有是精品在线观看| 欧美日韩亚洲高清精品| 亚洲美女搞黄在线观看| 国产色婷婷99| 久久毛片免费看一区二区三区| 亚洲,欧美,日韩| 亚洲第一区二区三区不卡| 视频在线观看一区二区三区| 街头女战士在线观看网站| 99热网站在线观看| 五月玫瑰六月丁香| 大片电影免费在线观看免费| 久久久精品94久久精品| 中国三级夫妇交换| 亚洲国产日韩一区二区| 菩萨蛮人人尽说江南好唐韦庄| 性色avwww在线观看| 欧美成人精品欧美一级黄| 一区二区日韩欧美中文字幕 | 中文字幕久久专区| 精品少妇黑人巨大在线播放| 22中文网久久字幕| www.色视频.com| 欧美国产精品一级二级三级| 亚洲高清免费不卡视频| 麻豆精品久久久久久蜜桃| 久久国产精品男人的天堂亚洲 | 亚洲,欧美,日韩| 天天躁夜夜躁狠狠久久av| 亚洲欧美成人综合另类久久久| 99九九线精品视频在线观看视频| 午夜福利影视在线免费观看| 日韩 亚洲 欧美在线| av线在线观看网站| 国产一区二区三区综合在线观看 | 久久久午夜欧美精品| 亚洲欧美成人综合另类久久久| 久久综合国产亚洲精品| 一本一本综合久久| 亚洲欧美日韩另类电影网站| 国产日韩欧美视频二区| 免费黄网站久久成人精品| 99久久精品一区二区三区| 夫妻性生交免费视频一级片| 91午夜精品亚洲一区二区三区| 精品久久久久久电影网| 人人澡人人妻人| 欧美 亚洲 国产 日韩一| 国产精品久久久久久精品古装| 久久狼人影院| 色94色欧美一区二区| 亚洲天堂av无毛| 丰满乱子伦码专区| 久久久久网色| 丝瓜视频免费看黄片| 亚洲欧美成人精品一区二区| 熟女电影av网| 日本wwww免费看| 国产高清国产精品国产三级| 国语对白做爰xxxⅹ性视频网站| 久久久久精品性色| 亚洲人成网站在线播| 成年女人在线观看亚洲视频| 大陆偷拍与自拍| 3wmmmm亚洲av在线观看| 国产成人91sexporn| 午夜久久久在线观看| 下体分泌物呈黄色| 少妇精品久久久久久久| 免费看光身美女| 亚洲精品亚洲一区二区| 国产成人91sexporn| 卡戴珊不雅视频在线播放| 麻豆乱淫一区二区| 欧美一级a爱片免费观看看| 国产伦理片在线播放av一区| 中文字幕人妻丝袜制服| 美女主播在线视频| 国产片内射在线| 亚洲国产色片| 下体分泌物呈黄色| 成年av动漫网址| av专区在线播放| 国产午夜精品一二区理论片| 亚洲成人av在线免费| 亚洲熟女精品中文字幕| 青春草亚洲视频在线观看| 女性被躁到高潮视频| 日韩精品免费视频一区二区三区 | 一本久久精品| 蜜桃国产av成人99| 在线观看一区二区三区激情| 亚洲成色77777| 国产熟女午夜一区二区三区 | 免费人成在线观看视频色| 亚洲国产精品一区二区三区在线| 久久久久精品久久久久真实原创| 亚洲欧美日韩卡通动漫| 人妻系列 视频| 日韩成人伦理影院| 日韩av不卡免费在线播放| 精品久久国产蜜桃| 精品一区二区三区视频在线| 亚洲少妇的诱惑av| 日韩,欧美,国产一区二区三区| 欧美日韩在线观看h| 青春草亚洲视频在线观看| av卡一久久| 欧美日韩亚洲高清精品| 亚洲国产精品999| 只有这里有精品99| 成人18禁高潮啪啪吃奶动态图 | 国产日韩欧美在线精品| 午夜福利在线观看免费完整高清在| 视频区图区小说| 国产av国产精品国产| 妹子高潮喷水视频| 日韩在线高清观看一区二区三区| 欧美日韩成人在线一区二区| 亚洲欧洲精品一区二区精品久久久 | www.av在线官网国产| 日本av免费视频播放| xxx大片免费视频| 一区二区三区免费毛片| 亚洲欧美精品自产自拍| 国产精品一国产av| 亚洲情色 制服丝袜| 久久久久久人妻| 看免费成人av毛片| 男人操女人黄网站| 亚洲婷婷狠狠爱综合网| 久久婷婷青草| 成年美女黄网站色视频大全免费 | 老司机影院成人| 亚洲精品一区蜜桃| 国产成人一区二区在线| 免费观看av网站的网址| 女人精品久久久久毛片| 日日爽夜夜爽网站| 国产无遮挡羞羞视频在线观看| 免费人成在线观看视频色| 欧美日韩综合久久久久久| 亚洲色图 男人天堂 中文字幕 | 国产精品.久久久| 亚洲精品日韩av片在线观看| 老司机亚洲免费影院| 亚洲精华国产精华液的使用体验| 天堂中文最新版在线下载| 七月丁香在线播放| 啦啦啦中文免费视频观看日本| 国产乱来视频区| 亚洲人成网站在线播| 久久久久久久久久人人人人人人| 午夜激情久久久久久久| 亚洲欧美清纯卡通| a级片在线免费高清观看视频| 国产国语露脸激情在线看| 99热6这里只有精品| 欧美最新免费一区二区三区| xxx大片免费视频| 欧美成人精品欧美一级黄| 成人毛片60女人毛片免费| 日日摸夜夜添夜夜爱| 成人二区视频| 女性生殖器流出的白浆| 免费av中文字幕在线| 亚洲少妇的诱惑av| 肉色欧美久久久久久久蜜桃| 一级二级三级毛片免费看| 高清毛片免费看| 亚洲精品日韩在线中文字幕| 性色avwww在线观看| 精品人妻熟女毛片av久久网站| 久久ye,这里只有精品| 在线观看人妻少妇| 日本wwww免费看| 18禁观看日本| 国产 一区精品| 国产日韩欧美视频二区| 边亲边吃奶的免费视频| 亚洲国产最新在线播放| 夜夜爽夜夜爽视频| 亚洲美女视频黄频| 国产精品久久久久久久久免| 日本与韩国留学比较| 丝袜脚勾引网站| 亚洲成色77777| 久久国产亚洲av麻豆专区| 免费看不卡的av| 久久人人爽人人爽人人片va| 18禁观看日本| 久久久久久久久久久久大奶| 久久国产亚洲av麻豆专区| 亚洲精品久久成人aⅴ小说 | 亚洲欧洲国产日韩| 国产黄片视频在线免费观看| 涩涩av久久男人的天堂| 精品国产乱码久久久久久小说| 97精品久久久久久久久久精品| 你懂的网址亚洲精品在线观看| 青青草视频在线视频观看| 日韩精品免费视频一区二区三区 | av在线观看视频网站免费| 99久久精品国产国产毛片| 亚洲国产日韩一区二区| 又粗又硬又长又爽又黄的视频| 欧美日韩亚洲高清精品| 视频在线观看一区二区三区| 99九九在线精品视频| 国产伦理片在线播放av一区| 色婷婷av一区二区三区视频| 亚洲国产精品成人久久小说| 尾随美女入室| 日本wwww免费看| 男人爽女人下面视频在线观看| 欧美97在线视频| 亚洲一区二区三区欧美精品| 91精品国产国语对白视频| 国产亚洲av片在线观看秒播厂| 纵有疾风起免费观看全集完整版| 9色porny在线观看| 涩涩av久久男人的天堂| 亚洲av中文av极速乱| 久久97久久精品| 免费人成在线观看视频色| 日本色播在线视频| 久久久久精品性色| 久久午夜综合久久蜜桃| 亚洲国产精品国产精品| 五月伊人婷婷丁香| 热re99久久精品国产66热6| 午夜日本视频在线| 曰老女人黄片| 满18在线观看网站| 国产成人av激情在线播放 | 国产免费现黄频在线看| 免费大片18禁| 午夜免费观看性视频| 国产一区亚洲一区在线观看| 美女视频免费永久观看网站| 国产高清不卡午夜福利| 国产视频首页在线观看| 老熟女久久久| 国产精品国产三级专区第一集| 我要看黄色一级片免费的| 人人妻人人爽人人添夜夜欢视频| 七月丁香在线播放| 欧美 日韩 精品 国产| 日本与韩国留学比较| 啦啦啦啦在线视频资源| 欧美精品一区二区免费开放| 日产精品乱码卡一卡2卡三| 观看美女的网站| 亚洲精品一区蜜桃| 国产成人精品无人区| 一本一本综合久久| 精品国产露脸久久av麻豆| 日日摸夜夜添夜夜爱| 精品一品国产午夜福利视频| 亚洲精品乱久久久久久| 中文字幕免费在线视频6| 欧美精品一区二区免费开放| 欧美日韩av久久| 免费av不卡在线播放| 国产片内射在线| videos熟女内射| 国产欧美日韩综合在线一区二区| 老司机亚洲免费影院| 女性生殖器流出的白浆| 成人18禁高潮啪啪吃奶动态图 | 在线观看免费视频网站a站| 亚洲av.av天堂| 国产成人aa在线观看| 永久免费av网站大全| 国产免费一区二区三区四区乱码| 日日摸夜夜添夜夜添av毛片| 一本色道久久久久久精品综合| 成人二区视频| 欧美精品人与动牲交sv欧美| 性色avwww在线观看| 国产成人精品福利久久| 伦理电影免费视频| 在线观看一区二区三区激情| 99热国产这里只有精品6| 国产视频内射| 22中文网久久字幕| 亚洲高清免费不卡视频| 97超碰精品成人国产| 国产精品免费大片| 国产色婷婷99| 欧美人与性动交α欧美精品济南到 | 大片电影免费在线观看免费| 中文字幕免费在线视频6| 日本黄大片高清| 日本av手机在线免费观看| 国产日韩欧美在线精品| 日本av免费视频播放| 人妻系列 视频| 日韩一本色道免费dvd| 热99国产精品久久久久久7| 免费观看在线日韩| av不卡在线播放| 欧美日韩综合久久久久久| av天堂久久9| 三上悠亚av全集在线观看| 日韩人妻高清精品专区| 热re99久久国产66热| 热99久久久久精品小说推荐| 国产精品人妻久久久影院| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品视频女| 亚洲av免费高清在线观看| 人体艺术视频欧美日本| 毛片一级片免费看久久久久| 欧美激情国产日韩精品一区| 久久99蜜桃精品久久| 国产黄色视频一区二区在线观看| 在线天堂最新版资源| videossex国产| videosex国产| 亚洲,一卡二卡三卡| 成年人免费黄色播放视频| 日韩成人av中文字幕在线观看| 国产又色又爽无遮挡免| 黑人高潮一二区| 成年av动漫网址| 日韩中字成人| 少妇熟女欧美另类| 18在线观看网站| 久久精品国产亚洲网站| av播播在线观看一区| 在线观看www视频免费| 一区二区三区精品91| 亚洲精品第二区| 亚洲精品自拍成人| 免费黄色在线免费观看| freevideosex欧美| 久久99热6这里只有精品| 欧美日韩精品成人综合77777| 自线自在国产av| 一本—道久久a久久精品蜜桃钙片| 国产日韩欧美在线精品| 男女国产视频网站| 日韩成人伦理影院| 啦啦啦中文免费视频观看日本| 色94色欧美一区二区| 我要看黄色一级片免费的| 伊人亚洲综合成人网| 亚洲欧美一区二区三区国产| 亚洲精品乱码久久久久久按摩| 2021少妇久久久久久久久久久| 亚洲欧美精品自产自拍| 99热6这里只有精品| 国产一区二区三区av在线| 午夜精品国产一区二区电影| 国产国语露脸激情在线看| 精品亚洲成国产av| 在线观看三级黄色| 亚洲色图综合在线观看| 国产视频内射| 国产午夜精品久久久久久一区二区三区| 青春草国产在线视频| av播播在线观看一区| 亚洲成人一二三区av| 国产熟女午夜一区二区三区 | 免费av不卡在线播放| 精品久久久精品久久久| 免费观看在线日韩| 亚洲精品日本国产第一区| 中文天堂在线官网| 国产视频内射| 成人黄色视频免费在线看| 亚洲人成网站在线播| 日本wwww免费看| 国产亚洲精品第一综合不卡 | 亚洲精品乱久久久久久| 免费少妇av软件| √禁漫天堂资源中文www| 亚洲国产日韩一区二区| 亚洲精品国产av成人精品| 纵有疾风起免费观看全集完整版| 国产乱来视频区| 国产国语露脸激情在线看| 精品亚洲乱码少妇综合久久| 我的女老师完整版在线观看| 男人添女人高潮全过程视频| 精品一区在线观看国产| 精品人妻熟女毛片av久久网站| 国产成人免费无遮挡视频| 美女脱内裤让男人舔精品视频| 中文乱码字字幕精品一区二区三区| 午夜福利网站1000一区二区三区| 在线 av 中文字幕| 只有这里有精品99| 中文字幕免费在线视频6| 亚洲伊人久久精品综合|