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

    高鐵擋風(fēng)墻對強(qiáng)風(fēng)沙響應(yīng)規(guī)律的數(shù)值模擬研究?

    2020-03-06 02:21:58岳海燕金阿芳高衛(wèi)強(qiáng)
    關(guān)鍵詞:擋風(fēng)墻沙粒坡腳

    岳海燕,金阿芳,高衛(wèi)強(qiáng)

    (新疆大學(xué)機(jī)械工程學(xué)院,新疆烏魯木齊830047)

    0 引言

    我國面積遼闊,鐵路運(yùn)輸具有不可替代的作用,20世紀(jì)90年代以來,中國對高速鐵路的建設(shè)開展了大量的科學(xué)研究和技術(shù)攻關(guān),我國鐵路形成了具有中國特色的高鐵技術(shù)體系[1].蘭新高速鐵路于2014年通車,是中國《中長期鐵路網(wǎng)規(guī)劃》的重點(diǎn)項目,也是歐亞大陸橋鐵路走廊不可或缺的部分.

    新疆植被稀少,沙源豐富,在大風(fēng)經(jīng)過的地帶就會形成沙塵天氣[2].大風(fēng)區(qū)內(nèi)自然條件惡劣,風(fēng)速高、風(fēng)期長、起風(fēng)速度快,尤其是百里風(fēng)區(qū)和三十里風(fēng)區(qū),屬于內(nèi)陸風(fēng)力最為強(qiáng)勁的地區(qū).蘭新線百里風(fēng)區(qū)每年都會出現(xiàn)速度超過40 m/s 的大風(fēng),實(shí)測最大風(fēng)速(瞬時)可達(dá)64 m/s,給鐵路運(yùn)輸造成了巨大的安全隱患.目前高速列車的行駛速度不斷提高,而列車由于大風(fēng)影響出現(xiàn)事故的情況頻發(fā)[3].強(qiáng)風(fēng)容易造成諸如火車傾覆,車窗玻璃破碎和鐵路交通設(shè)施損壞等事故,而高速列車普遍運(yùn)行速度快,對環(huán)境要求高,強(qiáng)橫風(fēng)將嚴(yán)重影響列車的安全運(yùn)營[4].

    根據(jù)以往的經(jīng)驗(yàn)總結(jié)和大風(fēng)對既有線路的危害分析,鐵路運(yùn)營需要解決列車傾覆、軌道積沙、沙石擊碎玻璃等諸多問題,以免影響鐵路安全運(yùn)營.在橫風(fēng)等惡劣天氣下,列車能否盡可能抵抗其影響從而安全運(yùn)營是衡量列車運(yùn)行的重要特征[5].為保障列車在大風(fēng)天氣下能正常運(yùn)行,國內(nèi)外多采取建造擋風(fēng)墻來減小橫風(fēng)對列車運(yùn)行影響的措施[6,7],擋風(fēng)墻是行之有效的防風(fēng)手段.

    研究人員對此進(jìn)行了眾多探索,MARIJO 等[8]運(yùn)用試驗(yàn)與數(shù)值計算方法研究了三維擋風(fēng)墻后紊流的特性.張軍平等[9]通過對路基周圍風(fēng)沙流場數(shù)值分析,發(fā)現(xiàn)在路基周圍會出現(xiàn)氣流運(yùn)動的速度分區(qū),并針對不同速度分區(qū)進(jìn)行研究.許志峰[10]研究了擋風(fēng)墻的疏透度,并提出使列車平穩(wěn)運(yùn)行的擋風(fēng)墻設(shè)置.張潔等[11]提出兩種優(yōu)化擋風(fēng)墻坡腳的方案,研究其對列車氣動力系數(shù)的影響.鄭曉靜等[12]對鐵路沿線已有擋風(fēng)墻在受強(qiáng)側(cè)風(fēng)作用時避免行進(jìn)列車傾覆方面的有效性以及擋風(fēng)墻沿線所存在的積沙問題進(jìn)行了定量分析.牛波等[13]通過數(shù)值模擬方法研究高速列車在不同車速和不同沙塵暴等級下運(yùn)行時,沙塵對擋風(fēng)墻高度和距離等參數(shù)的影響.

    上述文獻(xiàn)主要從設(shè)置擋風(fēng)墻后,防止列車傾覆的角度去研究,忽略了擋風(fēng)墻的阻沙效果及流場特性.蘭新高鐵線路所處的特殊地理環(huán)境和氣象條件,不僅要求沿線的擋風(fēng)墻具有擋風(fēng)作用外,還要具備防沙阻沙效果.本文利用CFD方法,借助湍流模型與DPM模型,模擬了不同初始風(fēng)速的強(qiáng)風(fēng)沙途經(jīng)不同參數(shù)的擋風(fēng)墻時路堤周圍的流場情況,繼而分析不同工況的擋風(fēng)墻對強(qiáng)風(fēng)沙的響應(yīng)規(guī)律,為今后西北多風(fēng)沙地區(qū)鐵道工程設(shè)置擋風(fēng)墻提供參考.

    1 數(shù)值模型的建立與分析

    本文主要計算風(fēng)沙條件下,擋風(fēng)墻周圍的流場特性,采用DPM離散相模型進(jìn)行數(shù)值模擬.控制方程包括質(zhì)量方程、動量方程與湍流模型方程,本文選擇k-ε 湍流模型.

    連續(xù)方程

    式中,ρ為密度,t為時間,u為速度,單位都為國際單位.

    動量方程

    式中,p為流體微單元上的壓強(qiáng),τxx,τyy,τzz等是因分子粘性作用而產(chǎn)生的作用在微單元體表面上的粘性應(yīng)力τ的分量,fx,fy,fz為三個方向的單位質(zhì)量力.

    K ?ε方程

    式中,k為湍動能,μt為湍動強(qiáng)度,ε為湍動耗散率,u為速度分量,G由平均速度梯度引起的湍動能產(chǎn)生,σt,σε為響應(yīng)的普朗特系數(shù),C為經(jīng)驗(yàn)常數(shù).

    2 計算模型及方法

    在大風(fēng)地區(qū)需要設(shè)置擋風(fēng)墻以減弱強(qiáng)風(fēng)對高速列車的影響,文中擋風(fēng)墻高度設(shè)定為H m.經(jīng)查閱文獻(xiàn),風(fēng)區(qū)多為季風(fēng)且風(fēng)向穩(wěn)定,因此,設(shè)計模型為單側(cè)擋風(fēng)墻.借助CFD方法對路堤及擋風(fēng)墻流場進(jìn)行數(shù)值計算,攜帶不同體積含量沙粒的風(fēng)沙流以一定的速度沿水平方向進(jìn)入流場左側(cè),因空氣質(zhì)量較輕,可忽略其重力對速度的影響.

    2.1 計算區(qū)域與邊界條件

    采用直立式擋風(fēng)墻,以蘭新高鐵沿線中某一典型擋風(fēng)墻幾何尺寸為設(shè)計原型建立數(shù)值模型,如圖1所示.設(shè)定路堤高度為3 m,路堤邊坡坡比為1:1,擋風(fēng)墻位置設(shè)置在距離路堤迎風(fēng)側(cè)坡腳3 m處,墻高為H m,路堤兩側(cè)坡腳都距離模型邊界20 m.

    為了讓氣流的繞流和流場充分發(fā)展,避免擋風(fēng)墻背風(fēng)側(cè)渦旋流對出口邊界條件的影響,數(shù)值計算模型長度為55 m,高度為15 m,寬度為25 m.

    模型左側(cè)為速度進(jìn)口邊界條件(VELOCITY-INLET),設(shè)置為含有沙粒的均勻風(fēng)場,風(fēng)速均勻,沙粒直徑一般在0.1 mm~0.25 mm;模型右側(cè)為出口,出口采用自由出口邊界條件(OUTLET-VENT),其余設(shè)置為無滑移壁面邊界條件(WALL).

    圖1 計算區(qū)域與邊界條件Fig 1 Calculating area and boundary conditions

    2.2 網(wǎng)格劃分及數(shù)值模擬工況

    本文為保證數(shù)值模擬的準(zhǔn)確性和計算效率,不過多占用計算資源,使用ICEM CFD軟件進(jìn)行網(wǎng)格劃分,如圖2所示.由風(fēng)沙流場的運(yùn)動特性可知,靠近擋風(fēng)墻,路堤的運(yùn)動更為復(fù)雜,因此在進(jìn)口、擋風(fēng)墻與路堤周圍處的網(wǎng)格劃分較為稠密,在出口、地面、壁面的網(wǎng)格較為稀疏,這樣既可以節(jié)約計算時間,加快數(shù)值模擬的收斂速度,又充分考慮到重點(diǎn)觀測部位以保證模擬的有效性和準(zhǔn)確性,提高了計算精度.計算工況如表1所示.

    圖2 計算域網(wǎng)格Fig 2 The grids of calculation field

    表1 計算工況Tab 1 Calculation conditions

    3 模擬結(jié)果及分析討論

    3.1 路堤周圍流場特性分析

    國內(nèi)關(guān)于風(fēng)沙流運(yùn)動特性研究的歷史較短,大多數(shù)研究機(jī)構(gòu)通過風(fēng)洞試驗(yàn)或者實(shí)際野外觀測的手段,本文使用數(shù)值模擬方法,研究未布設(shè)擋風(fēng)墻的情況下,路堤周圍的流場特性.氣流和沙粒從流場左側(cè)以20 m/s的速度平行于地表進(jìn)入,氣流的速度分布云圖和速度矢量圖如圖3所示.

    圖3 路堤周圍流場分布圖Fig 3 Flow field distribution around the embankment

    如圖3所示,氣流遇到路堤阻擋后,運(yùn)動方向發(fā)生改變,在路堤迎風(fēng)側(cè)路肩處開始形成速度分區(qū),在路堤上方速度變大,在迎風(fēng)側(cè)坡腳和背風(fēng)側(cè)處形成低速區(qū),且背風(fēng)側(cè)的風(fēng)速減小的范圍要遠(yuǎn)大于迎風(fēng)側(cè)坡腳處,由于受到迎風(fēng)側(cè)坡面的影響,風(fēng)速沿著坡面逐漸疊加,在迎風(fēng)側(cè)路肩上方位置風(fēng)速急劇增大,對路肩的侵蝕也最為嚴(yán)重,且沙粒也極易堆積在迎風(fēng)側(cè)與路堤處,影響列車的安全通行.

    圖4為路堤迎風(fēng)側(cè)坡腳和路肩處斷面上的風(fēng)速廓線圖,由圖4可知,在迎風(fēng)側(cè)坡腳處,自0 m~3.4 m范圍內(nèi),與路肩高度平齊處速度增長速度較快,從0 m/s增長至20 m/s左右,進(jìn)入高速區(qū),自此大約維持在22 m/s左右,風(fēng)速相對穩(wěn)定;在迎風(fēng)側(cè)路肩處,速度高達(dá)32 m/s,而后又沿著高度逐漸減小,直到距離路肩上方3 m處,速度減小緩慢,風(fēng)速達(dá)到相對穩(wěn)定的狀態(tài),此現(xiàn)象也與圖3描繪的現(xiàn)象相符.

    圖4 路堤不同位置處的風(fēng)速廓線圖Fig 4 Wind speed profile at different locations of the embankment

    3.2 擋風(fēng)墻背風(fēng)側(cè)流場特征分析

    為了研究擋風(fēng)墻背風(fēng)側(cè)流場特性,設(shè)定入口風(fēng)速為20 m/s,擋風(fēng)墻高度為3 m,對流場進(jìn)行模擬,結(jié)果如圖5所示.

    由圖5(a)和圖5(b)可以看出,攜帶沙粒的氣流在受到擋風(fēng)墻和路堤的阻礙時,氣流運(yùn)動方向發(fā)生明顯變化,過流斷面減小,氣流被壓縮發(fā)生繞流,且攜沙氣流的速度也形成較為明顯的分區(qū),在擋風(fēng)墻頂端上方區(qū)域,形成高速區(qū),在擋風(fēng)墻與路堤迎風(fēng)側(cè)之間、路堤上方、路堤背風(fēng)側(cè)形成低速區(qū),在擋風(fēng)墻與路堤迎風(fēng)側(cè)之間形成渦流,此處沙粒運(yùn)動也較為復(fù)雜,速度減小,一部分沙??赡軙逊e,同時也起到了凈化氣流的作用,且背風(fēng)側(cè)坡腳氣流速度減小的范圍要大于迎風(fēng)側(cè)坡腳的范圍,當(dāng)氣流繞過擋風(fēng)墻后,過流斷面增大,氣流會攜帶一部分沙粒越過路堤上方,一些吹向遠(yuǎn)處,一些落在擋風(fēng)墻背風(fēng)側(cè),形成堆積.經(jīng)過一段時間后,速度與方向又逐漸恢復(fù)到初始狀態(tài),繼續(xù)流動.

    圖5 設(shè)置3 m高擋風(fēng)墻周圍流場圖Fig 5 Flow field diagram around a 3 m high windshield wall

    3.3 不同高度的擋風(fēng)墻背風(fēng)側(cè)流場特征

    在高速鐵路軌道旁設(shè)置擋風(fēng)墻,在保護(hù)列車盡可能少的受到風(fēng)沙侵蝕的同時,又可減少鐵軌路面上的積沙面積,不讓列車運(yùn)行受到積沙的影響.擋風(fēng)墻的高度也直接影響著工程造價和保護(hù)列車及軌道的效果.擋風(fēng)墻越高,工程造價越高,需要的抗傾覆力矩也越大,擋風(fēng)墻太低起不到應(yīng)有的效果,所以,一個合理的擋風(fēng)墻高度起到了至關(guān)重要的作用.

    風(fēng)沙流是指含有沙粒的氣流,其形成要有較為密集的沙粒和一定速度的風(fēng)力.風(fēng)攜帶沙粒吹向擋風(fēng)墻周圍,研究不同高度條件下?lián)躏L(fēng)墻周圍的流場主要是觀察氣流通過擋風(fēng)墻之后的速度變化情況.為了研究不同高度條件下?lián)躏L(fēng)墻背風(fēng)側(cè)的風(fēng)速特征,分別取擋風(fēng)墻高度為2.5 m,3 m,3.5 m,4 m,攜帶沙粒的風(fēng)速分別以20 m/s,30 m/s,40 m/s的工況進(jìn)行數(shù)值模擬.截取同一個面處,對距離擋風(fēng)墻背風(fēng)側(cè)1.5 m處進(jìn)行數(shù)據(jù)提取.

    圖6 擋風(fēng)墻高度2.5 m時不同初始風(fēng)速的風(fēng)速廓線圖Fig 6 Wind speed profile of different initial windspeeds when the heigh of the windshield is 2.5 m

    圖7 擋風(fēng)墻高度3 m時不同初始風(fēng)速的風(fēng)速廓線圖Fig 7 Wind speed profile of different initial windspeeds when the height of the windshield is 3 m

    研究擋風(fēng)墻背風(fēng)側(cè)風(fēng)速變化情況有著重要的作用,當(dāng)擋風(fēng)墻高度為2.5 m時,如圖6所示,三種初始風(fēng)速的變化曲線規(guī)律大致相似,在0~2.5 m的高度范圍內(nèi),當(dāng)初始風(fēng)速為20 m/s時,最大風(fēng)速為11 m/s左右;初始風(fēng)速為30 m/s時,最大速度為16 m/s左右;初始風(fēng)速為40 m/s時,最大速度為21 m/s左右.在2.5 m~15 m的高度范圍內(nèi),風(fēng)速廓線都經(jīng)歷了先減小后增大的變化過程,初始風(fēng)速分別為20 m/s,30 m/s,40 m/s時,期間最小速度分別為6 m/s,8 m/s,12 m/s,在高速區(qū)時,速度又分別最大增加至31 m/s,47 m/s,63 m/s,而后隨著風(fēng)沙的流動,慢慢減速至初始速度.

    當(dāng)擋風(fēng)墻高度為3 m時,如圖7所示,在0~5 m的高度內(nèi),三種初始風(fēng)速的曲線區(qū)別很小,風(fēng)速也普遍小,初始風(fēng)速為20 m/s時,最大速度為12 m/s;初始風(fēng)速為30 m/s時,最大風(fēng)速為17 m/s;初始風(fēng)速為40 m/s時,最大風(fēng)速為15 m/s.在5 m~15 m的范圍內(nèi),速度先呈指數(shù)規(guī)律增大,分別增大至32 m/s,48 m/s,64 m/s,在擋風(fēng)墻高度2 m以上為加速區(qū),隨著高度的增加風(fēng)速呈指數(shù)分布規(guī)律上升,然后氣流速度逐漸變的平緩,漸漸趨于初始風(fēng)速.

    擋風(fēng)墻高度為3.5 m時,如圖8所示,與墻高為3 m的情況相比較,大致規(guī)律相同,只是在0 ~6 m的高度范圍內(nèi),三種初始風(fēng)速的曲線略有波動,速度都較小,速度保持在10 m/s以內(nèi);在6 m~15 m的高度范圍內(nèi),在高度為7 m的時候,都達(dá)到了流場內(nèi)的最大風(fēng)速.初始風(fēng)速為20 m/s時,風(fēng)速最大為34 m/s;初始風(fēng)速為30 m/s時,風(fēng)速最大為50 m/s;初始風(fēng)速為40 m/s時,風(fēng)速最大為68 m/s,然后隨著高度的增加逐漸遠(yuǎn)離路堤,速度又漸漸較為平緩的降低至初始風(fēng)速.

    當(dāng)擋風(fēng)墻高度為4 m時,如圖9所示,在距地面0 ~7 m的高度范圍內(nèi),速度波動平緩且風(fēng)速較小,最大風(fēng)速基本在11 m/s之內(nèi);在高度7 m~15 m范圍內(nèi),風(fēng)速先是呈指數(shù)規(guī)律上升.初始風(fēng)速為20 m/s時,最大風(fēng)速為35 m/s;初始風(fēng)速為30 m/s時,最大風(fēng)速為54 m/s;初始風(fēng)速為40 m/s時,最大風(fēng)速為72 m/s,然后隨著高度的增加逐漸遠(yuǎn)離路堤,速度又漸漸較為平緩的降低至初始風(fēng)速左右.

    圖8 擋風(fēng)墻高度3.5 m時不同初始風(fēng)速的風(fēng)速廓線圖Fig 8 Wind speed profile of different initial wind speeds when the height of the windshield is 3.5m

    圖9 擋風(fēng)墻高度4 m時不同初始風(fēng)速的風(fēng)速廓線Fig 9 Wind speed profile of different initial wind speeds when the height of the windshield is 4m

    3.4 擋風(fēng)墻合理高度研究

    擋風(fēng)墻的高度對傾覆力矩、建造成本、阻沙效果都有至關(guān)重要的影響.在前文中,主要討論了不同風(fēng)速下,各高度擋風(fēng)墻的背風(fēng)側(cè)流場特性,為了研究擋風(fēng)墻的合理高度,將四種不同高度的擋風(fēng)墻同時做對比,進(jìn)行合理高度的比較.

    如圖10所示,當(dāng)初始風(fēng)速為20 m/s時,可以看出,擋風(fēng)墻高度為2.5 m的情況與其余三種情況區(qū)別比較明顯,流場的速度波動更大,在墻的高度范圍內(nèi)速度的波動更為明顯,穩(wěn)定性較低,且離墻體近的地方風(fēng)速較大,不利于列車通行.擋風(fēng)墻高度為3 m,3.5 m,4 m時,通過比較,可以看出當(dāng)擋風(fēng)墻高度為4 m時,在相對離墻較遠(yuǎn)且遠(yuǎn)離高鐵的地方速度開始增長,安全性更高.當(dāng)初始風(fēng)速為30 m/s,40 m/s時,變化規(guī)律相似,從而說明擋風(fēng)墻高度為4 m時,對列車的保護(hù)性能更好,流場更為穩(wěn)定.

    圖10 不同初始風(fēng)速下?lián)躏L(fēng)墻風(fēng)速廓線Fig 10 Wind speed profile of windshield wall under different initial wind speeds

    4 結(jié)論

    在以上的研究分析中,未布設(shè)擋風(fēng)墻與布設(shè)擋風(fēng)墻之后的流場特性有明顯不同之處,擋風(fēng)墻改變了路堤周圍的流場分布,阻擋沙粒的沉積.

    (1)在未布設(shè)擋風(fēng)墻時,攜沙氣流遇到路堤后在路堤的路肩處形成分區(qū),在迎風(fēng)側(cè)路肩上方位置風(fēng)速急劇增大,對路肩的侵蝕最為嚴(yán)重,沙粒也極易堆積在迎風(fēng)側(cè)與路堤處,影響列車安全運(yùn)行.

    (2)在布設(shè)擋風(fēng)墻后,攜沙氣流遇到擋風(fēng)墻形成速度分區(qū),高速區(qū)由路堤的路肩前移到擋風(fēng)墻上方區(qū)域,且背風(fēng)側(cè)坡腳氣流速度減小的范圍要大于迎風(fēng)側(cè)坡腳的范圍.同時在擋風(fēng)墻與路堤迎風(fēng)側(cè)之間形成渦流,一部分沙??赡軙虼硕逊e,同時也起到了凈化氣流的作用,也減少了沙粒沉積在軌道的數(shù)量.

    (3)不同風(fēng)速的攜沙氣流通過擋風(fēng)墻時,風(fēng)速垂直廓線的規(guī)律大致相同,呈現(xiàn)先緩慢波動增長后指數(shù)規(guī)律增長的趨勢,最后逐漸靠近初始風(fēng)速.初始風(fēng)速越高,達(dá)到的最大風(fēng)速也越大.且出現(xiàn)的最大風(fēng)速約是初始風(fēng)速的1.7倍.

    (4)初始風(fēng)速一定,攜沙氣流通過不同高度擋風(fēng)墻時,風(fēng)速垂直廓線的規(guī)律大致相同.擋風(fēng)墻高度越高,在自身高度范圍內(nèi),速度波動越小.

    猜你喜歡
    擋風(fēng)墻沙粒坡腳
    軟土路基施工對鄰近管線的影響及保護(hù)措施
    軟弱結(jié)構(gòu)面位置對巖質(zhì)順傾邊坡穩(wěn)定性的影響
    擋風(fēng)墻在篦冷機(jī)中的應(yīng)用
    既有南疆線大風(fēng)區(qū)預(yù)制拼裝化擋風(fēng)墻設(shè)計研究
    沙粒和水珠
    想看山的小沙粒
    想看山的小沙粒
    順層巖質(zhì)邊坡坡腳開挖穩(wěn)定性數(shù)值模擬分析
    直接空冷體系擋風(fēng)墻體型系數(shù)研究
    坡腳圓法驗(yàn)算路基邊坡的穩(wěn)定性
    午夜久久久在线观看| 久久99热这里只有精品18| 长腿黑丝高跟| 亚洲人成电影免费在线| 看免费av毛片| 国内少妇人妻偷人精品xxx网站 | 亚洲三区欧美一区| 久热这里只有精品99| 黄色视频,在线免费观看| 午夜精品久久久久久毛片777| 国产久久久一区二区三区| 国产亚洲精品一区二区www| 国产99白浆流出| 一a级毛片在线观看| 丝袜美腿诱惑在线| 在线av久久热| 99国产精品一区二区三区| 日本五十路高清| 女性生殖器流出的白浆| x7x7x7水蜜桃| 欧美日本视频| 在线视频色国产色| 黄频高清免费视频| 精品午夜福利视频在线观看一区| 免费在线观看视频国产中文字幕亚洲| 久久久久国内视频| 在线av久久热| 欧美日韩精品网址| 亚洲激情在线av| 亚洲黑人精品在线| 亚洲国产精品久久男人天堂| 欧美激情久久久久久爽电影| avwww免费| 日本a在线网址| 狂野欧美激情性xxxx| 韩国av一区二区三区四区| 波多野结衣高清无吗| 欧美激情极品国产一区二区三区| 亚洲成人久久爱视频| 久久青草综合色| 久久精品91蜜桃| 丰满的人妻完整版| 成年人黄色毛片网站| 最近最新中文字幕大全电影3 | 19禁男女啪啪无遮挡网站| 久久九九热精品免费| 级片在线观看| 在线观看免费视频日本深夜| 国产精品,欧美在线| 午夜福利在线观看吧| 日本在线视频免费播放| 可以免费在线观看a视频的电影网站| 亚洲第一电影网av| 久热爱精品视频在线9| 久久狼人影院| 99riav亚洲国产免费| 国产人伦9x9x在线观看| 色哟哟哟哟哟哟| 日本精品一区二区三区蜜桃| 女警被强在线播放| 欧美日本视频| 国内揄拍国产精品人妻在线 | 国产亚洲精品av在线| 日韩欧美国产一区二区入口| 国产亚洲精品久久久久5区| 首页视频小说图片口味搜索| 亚洲色图av天堂| 日韩成人在线观看一区二区三区| 啪啪无遮挡十八禁网站| 成人18禁高潮啪啪吃奶动态图| 黑人巨大精品欧美一区二区mp4| 国产精品久久久久久精品电影 | 午夜福利高清视频| 久久久精品欧美日韩精品| 亚洲欧美激情综合另类| 手机成人av网站| av在线播放免费不卡| 人成视频在线观看免费观看| 男女那种视频在线观看| 久久久久亚洲av毛片大全| 久久久精品欧美日韩精品| 国产亚洲欧美98| 熟妇人妻久久中文字幕3abv| 日韩 欧美 亚洲 中文字幕| 久久精品国产亚洲av高清一级| 欧美日韩亚洲国产一区二区在线观看| 亚洲熟女毛片儿| 久久久久久久久久黄片| 99热6这里只有精品| 9191精品国产免费久久| 国产一区在线观看成人免费| 一个人观看的视频www高清免费观看 | 每晚都被弄得嗷嗷叫到高潮| 久久天堂一区二区三区四区| 欧美黑人巨大hd| 久久草成人影院| 久久国产精品人妻蜜桃| 日韩精品青青久久久久久| 午夜福利免费观看在线| 91麻豆av在线| 动漫黄色视频在线观看| 免费高清视频大片| 十分钟在线观看高清视频www| 熟妇人妻久久中文字幕3abv| 黄片小视频在线播放| 午夜精品在线福利| 精品一区二区三区视频在线观看免费| 国产三级黄色录像| 国产v大片淫在线免费观看| 手机成人av网站| 18禁国产床啪视频网站| 青草久久国产| 国产熟女xx| 巨乳人妻的诱惑在线观看| 国产精品久久久久久人妻精品电影| 亚洲国产中文字幕在线视频| 日韩欧美 国产精品| 久久久久久大精品| √禁漫天堂资源中文www| 国产精品亚洲一级av第二区| 琪琪午夜伦伦电影理论片6080| 免费在线观看成人毛片| 国产精品99久久99久久久不卡| 亚洲第一青青草原| 99久久国产精品久久久| 亚洲色图av天堂| 岛国视频午夜一区免费看| 欧美日韩一级在线毛片| 国产亚洲精品一区二区www| 精品国产一区二区三区四区第35| 久久久久九九精品影院| 国产熟女xx| 在线观看一区二区三区| 午夜福利视频1000在线观看| 国产精品野战在线观看| 国产一区在线观看成人免费| 国产真实乱freesex| 香蕉久久夜色| 免费在线观看完整版高清| 一级作爱视频免费观看| 欧美激情久久久久久爽电影| 国产一区二区在线av高清观看| av片东京热男人的天堂| 亚洲专区国产一区二区| 亚洲av成人av| 免费在线观看亚洲国产| 国产精品亚洲一级av第二区| 国内毛片毛片毛片毛片毛片| www.熟女人妻精品国产| av中文乱码字幕在线| 久久精品国产亚洲av高清一级| 精品久久久久久久久久久久久 | 午夜免费激情av| 中文字幕av电影在线播放| 亚洲自拍偷在线| 最近最新免费中文字幕在线| 欧美一级a爱片免费观看看 | 老鸭窝网址在线观看| 男女那种视频在线观看| 巨乳人妻的诱惑在线观看| 99久久精品国产亚洲精品| 99在线视频只有这里精品首页| 99久久99久久久精品蜜桃| 国产又色又爽无遮挡免费看| 欧美成狂野欧美在线观看| 国产亚洲精品综合一区在线观看 | 岛国视频午夜一区免费看| 国内精品久久久久久久电影| 中文亚洲av片在线观看爽| 国产一区二区在线av高清观看| 精品久久蜜臀av无| 亚洲精品在线美女| 欧美丝袜亚洲另类 | 真人做人爱边吃奶动态| 久久国产精品男人的天堂亚洲| a在线观看视频网站| 国产精品野战在线观看| 国产av一区二区精品久久| 麻豆成人午夜福利视频| 曰老女人黄片| 国产黄a三级三级三级人| 狂野欧美激情性xxxx| 亚洲国产高清在线一区二区三 | 中文字幕久久专区| 久久99热这里只有精品18| 亚洲人成电影免费在线| 黑人巨大精品欧美一区二区mp4| 成人免费观看视频高清| 欧美中文综合在线视频| 淫秽高清视频在线观看| 国产精品永久免费网站| 波多野结衣高清无吗| 中亚洲国语对白在线视频| 日日爽夜夜爽网站| 亚洲激情在线av| 精品一区二区三区四区五区乱码| 999久久久精品免费观看国产| 丝袜美腿诱惑在线| 人妻丰满熟妇av一区二区三区| 成人三级黄色视频| 久久久国产成人免费| 国产精品乱码一区二三区的特点| 色综合站精品国产| 久久国产精品人妻蜜桃| av欧美777| 久久香蕉激情| 国产亚洲精品久久久久久毛片| 午夜福利免费观看在线| 欧美一级a爱片免费观看看 | 十八禁人妻一区二区| 99久久99久久久精品蜜桃| 免费av毛片视频| 亚洲欧美激情综合另类| 丝袜在线中文字幕| 1024视频免费在线观看| 欧美日韩瑟瑟在线播放| 在线天堂中文资源库| 亚洲av电影不卡..在线观看| 99久久久亚洲精品蜜臀av| 亚洲一码二码三码区别大吗| 久久久久免费精品人妻一区二区 | 亚洲精品国产一区二区精华液| 久久久久久人人人人人| 大香蕉久久成人网| 老司机在亚洲福利影院| 性色av乱码一区二区三区2| 高潮久久久久久久久久久不卡| 美女免费视频网站| 在线观看日韩欧美| 一边摸一边做爽爽视频免费| 中文字幕另类日韩欧美亚洲嫩草| 看免费av毛片| 宅男免费午夜| 中出人妻视频一区二区| 男人舔女人的私密视频| 国产亚洲av嫩草精品影院| 一本大道久久a久久精品| 老汉色av国产亚洲站长工具| 99热这里只有精品一区 | 亚洲一区二区三区不卡视频| 中文字幕人妻丝袜一区二区| 美女扒开内裤让男人捅视频| 波多野结衣av一区二区av| 一级作爱视频免费观看| 久久精品91蜜桃| 亚洲av熟女| 久久国产亚洲av麻豆专区| 国产高清激情床上av| 欧美色欧美亚洲另类二区| 国产一区二区三区视频了| 亚洲一码二码三码区别大吗| 一本大道久久a久久精品| 亚洲精品美女久久av网站| 国产激情久久老熟女| 国产精品影院久久| 人人妻人人澡人人看| 女同久久另类99精品国产91| 69av精品久久久久久| 亚洲av片天天在线观看| 国产三级黄色录像| 国产成人av激情在线播放| av欧美777| 国产视频内射| 日韩精品青青久久久久久| 国产精华一区二区三区| 欧美成人午夜精品| 久久亚洲真实| 国产成+人综合+亚洲专区| АⅤ资源中文在线天堂| 黄色a级毛片大全视频| av欧美777| 日本撒尿小便嘘嘘汇集6| 国产一区二区三区视频了| 欧美黑人巨大hd| 欧美乱色亚洲激情| 国产精品精品国产色婷婷| 免费观看精品视频网站| 国产成年人精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利欧美成人| 欧美成人免费av一区二区三区| 999久久久国产精品视频| 在线观看午夜福利视频| 久热爱精品视频在线9| 51午夜福利影视在线观看| 人人妻人人看人人澡| 日本精品一区二区三区蜜桃| 国产三级在线视频| 国内精品久久久久久久电影| 一区二区三区高清视频在线| 一本一本综合久久| 成人特级黄色片久久久久久久| 狠狠狠狠99中文字幕| 人人妻人人看人人澡| 亚洲精品美女久久av网站| 一本精品99久久精品77| 午夜免费观看网址| 天堂影院成人在线观看| 亚洲激情在线av| 美女扒开内裤让男人捅视频| 男人舔女人下体高潮全视频| 国产97色在线日韩免费| 91成人精品电影| 可以在线观看的亚洲视频| 亚洲国产高清在线一区二区三 | 欧美绝顶高潮抽搐喷水| 亚洲狠狠婷婷综合久久图片| 亚洲午夜理论影院| 黑丝袜美女国产一区| 精品国内亚洲2022精品成人| 亚洲精品美女久久久久99蜜臀| svipshipincom国产片| 久久国产精品男人的天堂亚洲| 日韩一卡2卡3卡4卡2021年| 老熟妇仑乱视频hdxx| 欧美乱妇无乱码| 又紧又爽又黄一区二区| 色精品久久人妻99蜜桃| 成年版毛片免费区| 免费人成视频x8x8入口观看| 亚洲一区中文字幕在线| 亚洲aⅴ乱码一区二区在线播放 | 久久精品国产亚洲av香蕉五月| 精品久久蜜臀av无| 中文字幕人成人乱码亚洲影| 精品午夜福利视频在线观看一区| 99re在线观看精品视频| 国产v大片淫在线免费观看| 久久人人精品亚洲av| 亚洲欧美激情综合另类| 国产野战对白在线观看| 欧美丝袜亚洲另类 | 国内少妇人妻偷人精品xxx网站 | 男人操女人黄网站| 一区二区三区高清视频在线| 精品熟女少妇八av免费久了| 伊人久久大香线蕉亚洲五| 麻豆国产av国片精品| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久九九精品二区国产 | 亚洲在线自拍视频| 一进一出好大好爽视频| 亚洲国产精品久久男人天堂| 欧美午夜高清在线| 九色国产91popny在线| 久久天堂一区二区三区四区| 精品第一国产精品| 波多野结衣高清作品| 男女那种视频在线观看| 十八禁网站免费在线| 黑人欧美特级aaaaaa片| 免费在线观看完整版高清| 国产久久久一区二区三区| 亚洲精品久久成人aⅴ小说| 亚洲avbb在线观看| 欧美黑人精品巨大| 国产av一区在线观看免费| 久久精品aⅴ一区二区三区四区| 老司机福利观看| 中文字幕人成人乱码亚洲影| svipshipincom国产片| 久久亚洲精品不卡| 久久伊人香网站| 久久狼人影院| 亚洲国产中文字幕在线视频| 久久精品夜夜夜夜夜久久蜜豆 | 淫妇啪啪啪对白视频| 欧美黑人欧美精品刺激| 两个人看的免费小视频| 亚洲熟妇熟女久久| 国产精品电影一区二区三区| 国产一区二区激情短视频| 50天的宝宝边吃奶边哭怎么回事| 在线天堂中文资源库| 99精品在免费线老司机午夜| 视频区欧美日本亚洲| 亚洲av中文字字幕乱码综合 | 不卡av一区二区三区| 51午夜福利影视在线观看| 欧美日韩瑟瑟在线播放| 久久精品91无色码中文字幕| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久国产高清桃花| 少妇 在线观看| 欧美成人一区二区免费高清观看 | 国产乱人伦免费视频| 亚洲七黄色美女视频| 欧美日韩瑟瑟在线播放| 最近最新免费中文字幕在线| 夜夜看夜夜爽夜夜摸| 成人手机av| 黄色 视频免费看| 麻豆国产av国片精品| 国产精品综合久久久久久久免费| 久久这里只有精品19| 欧美日本亚洲视频在线播放| 欧美日韩亚洲国产一区二区在线观看| 两个人免费观看高清视频| 亚洲第一电影网av| 99热只有精品国产| 免费一级毛片在线播放高清视频| or卡值多少钱| 无限看片的www在线观看| 国产黄色小视频在线观看| 超碰成人久久| 老鸭窝网址在线观看| 国产精品98久久久久久宅男小说| 成人手机av| 午夜视频精品福利| 亚洲人成伊人成综合网2020| 97人妻精品一区二区三区麻豆 | 日韩免费av在线播放| 精品无人区乱码1区二区| 国产男靠女视频免费网站| 国产单亲对白刺激| 久久99热这里只有精品18| 久久热在线av| 久久久国产精品麻豆| 一级毛片精品| 亚洲九九香蕉| 午夜激情av网站| 亚洲精品久久国产高清桃花| 99久久久亚洲精品蜜臀av| 久久这里只有精品19| 亚洲国产毛片av蜜桃av| 久久伊人香网站| 岛国视频午夜一区免费看| 国产激情欧美一区二区| 91av网站免费观看| 一二三四在线观看免费中文在| 高清毛片免费观看视频网站| 999久久久国产精品视频| 日本撒尿小便嘘嘘汇集6| 美女免费视频网站| 欧美大码av| 日韩三级视频一区二区三区| 国产熟女xx| 满18在线观看网站| 黄网站色视频无遮挡免费观看| 99精品久久久久人妻精品| 日韩欧美国产一区二区入口| 午夜福利成人在线免费观看| 88av欧美| 99久久99久久久精品蜜桃| 日本a在线网址| 欧美精品亚洲一区二区| 最近在线观看免费完整版| 久久欧美精品欧美久久欧美| 色av中文字幕| 亚洲男人天堂网一区| 成人精品一区二区免费| 国产成人精品久久二区二区91| 亚洲av熟女| 老司机在亚洲福利影院| 香蕉国产在线看| 此物有八面人人有两片| 欧美成人免费av一区二区三区| 亚洲激情在线av| 中文字幕av电影在线播放| 精品国产乱码久久久久久男人| 又黄又爽又免费观看的视频| 成年人黄色毛片网站| 国产精品久久久久久亚洲av鲁大| 精品午夜福利视频在线观看一区| 免费在线观看完整版高清| 一级毛片高清免费大全| 亚洲一码二码三码区别大吗| 国产一级毛片七仙女欲春2 | 国产精品久久久久久精品电影 | av视频在线观看入口| 变态另类成人亚洲欧美熟女| 一级作爱视频免费观看| 国产精品久久久av美女十八| 亚洲欧洲精品一区二区精品久久久| 一二三四社区在线视频社区8| 亚洲国产日韩欧美精品在线观看 | 老熟妇仑乱视频hdxx| 女人被狂操c到高潮| 他把我摸到了高潮在线观看| 成人精品一区二区免费| 最近最新中文字幕大全免费视频| 国产午夜福利久久久久久| 国产成人av教育| 欧美在线黄色| 最近最新中文字幕大全电影3 | 免费观看精品视频网站| 免费女性裸体啪啪无遮挡网站| 麻豆成人av在线观看| 亚洲成av人片免费观看| 亚洲自拍偷在线| 欧美大码av| 国产精品爽爽va在线观看网站 | 俄罗斯特黄特色一大片| 女同久久另类99精品国产91| 非洲黑人性xxxx精品又粗又长| 欧美av亚洲av综合av国产av| 中文在线观看免费www的网站 | 美女 人体艺术 gogo| 在线观看一区二区三区| 看免费av毛片| 欧美成人性av电影在线观看| 色哟哟哟哟哟哟| 欧美中文日本在线观看视频| 亚洲av美国av| 一本综合久久免费| 免费看十八禁软件| 18禁观看日本| 精品久久蜜臀av无| 国产精品二区激情视频| 成人永久免费在线观看视频| 久久久久国内视频| 此物有八面人人有两片| 久久久久精品国产欧美久久久| 中文字幕精品亚洲无线码一区 | 成人午夜高清在线视频 | 久久精品成人免费网站| 亚洲国产高清在线一区二区三 | 99国产精品99久久久久| 91成年电影在线观看| 国产成+人综合+亚洲专区| 夜夜爽天天搞| 欧美三级亚洲精品| 国产麻豆成人av免费视频| 91成人精品电影| 国产成人欧美在线观看| 女人高潮潮喷娇喘18禁视频| 99久久精品国产亚洲精品| 午夜影院日韩av| 可以免费在线观看a视频的电影网站| 日韩成人在线观看一区二区三区| 亚洲成a人片在线一区二区| 午夜福利在线在线| 热99re8久久精品国产| 成人国产一区最新在线观看| 日韩有码中文字幕| 欧美一区二区精品小视频在线| 国产麻豆成人av免费视频| 国内少妇人妻偷人精品xxx网站 | 亚洲精品久久成人aⅴ小说| 亚洲国产精品久久男人天堂| av在线播放免费不卡| 成人三级做爰电影| 亚洲片人在线观看| 国产人伦9x9x在线观看| 国产一卡二卡三卡精品| xxx96com| 国产亚洲欧美在线一区二区| 中文资源天堂在线| 视频在线观看一区二区三区| 嫩草影视91久久| 精品久久久久久成人av| 色综合站精品国产| 国产亚洲av高清不卡| 精品无人区乱码1区二区| 黄色女人牲交| 搡老熟女国产l中国老女人| 国产亚洲av嫩草精品影院| 丁香欧美五月| 日日夜夜操网爽| 啪啪无遮挡十八禁网站| 久久久久亚洲av毛片大全| 国产片内射在线| 国产精品精品国产色婷婷| 精品国产亚洲在线| 成熟少妇高潮喷水视频| 欧美性猛交黑人性爽| 精品少妇一区二区三区视频日本电影| 最新在线观看一区二区三区| 亚洲精品色激情综合| 女性被躁到高潮视频| 18禁黄网站禁片午夜丰满| 在线永久观看黄色视频| 成人永久免费在线观看视频| 久9热在线精品视频| 国产高清激情床上av| 久久精品国产综合久久久| 一本久久中文字幕| 欧美在线黄色| 久久久久久久久中文| 男人的好看免费观看在线视频 | 最近最新中文字幕大全免费视频| 精品电影一区二区在线| 一本精品99久久精品77| 日本a在线网址| 999久久久国产精品视频| 天堂影院成人在线观看| 性色av乱码一区二区三区2| 欧美一区二区精品小视频在线| 欧美国产日韩亚洲一区| 亚洲九九香蕉| 国产成人欧美| 国内久久婷婷六月综合欲色啪| 丝袜美腿诱惑在线| 久久久久久人人人人人| 正在播放国产对白刺激| 51午夜福利影视在线观看| 亚洲 国产 在线| xxx96com| 免费搜索国产男女视频| 2021天堂中文幕一二区在线观 | 叶爱在线成人免费视频播放| 黄色 视频免费看| 伦理电影免费视频| а√天堂www在线а√下载| 日韩精品青青久久久久久| 91老司机精品| 免费人成视频x8x8入口观看| 免费看a级黄色片| 啦啦啦韩国在线观看视频| 麻豆av在线久日| 18禁黄网站禁片免费观看直播| 99国产精品99久久久久| 91国产中文字幕| 亚洲片人在线观看| 无限看片的www在线观看|