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

    下?lián)舯┝鞒隽鞫畏欠€(wěn)態(tài)風場的大氣邊界層風洞模擬

    2022-01-21 10:12:48鐘永力晏致濤楊小剛
    實驗流體力學 2021年6期
    關(guān)鍵詞:暴流時變風場

    鐘永力,晏致濤,,*,李 妍,楊小剛,蔣 森

    1.重慶科技學院 建筑工程學院,重慶 401331;2.重慶大學 土木工程學院,重慶 400045

    0 引 言

    下?lián)舯┝魇且环N在雷暴天氣中由強下沉氣流猛烈沖擊地面而形成并沿地表傳播的近地面短時破壞性強風[1]。在雷暴天氣中,發(fā)生微下?lián)舯┝鞯母怕士蛇_60%~70%,產(chǎn)生的近地面強風最大風速可超過60 m/s。在世界范圍內(nèi),下?lián)舯┝髟斐闪舜罅抗こ探Y(jié)構(gòu)物的破壞[2]。

    下?lián)舯┝骶哂休^強的突發(fā)性,較難記錄實際的下?lián)舯┝魇录?,因此,風洞試驗是研究下?lián)舯┝鞯闹饕侄巍嶒炇已芯肯聯(lián)舯┝髌毡椴捎脹_擊射流裝置。早期,沖擊射流裝置主要用于模擬平穩(wěn)風場特性,由于水平放置,一定程度上會受到地面的影響。之后,通過將噴嘴懸掛于足夠高處,盡量排除了地面的影響[3]。此類試驗在后續(xù)研究中得到了進一步改進,形成了移動式壁面射流裝置[4-5]。同時,可以通過控制沖擊射流閥門的開啟速度來實現(xiàn)非穩(wěn)態(tài)的下?lián)舯┝髂M[6]。國內(nèi)對下?lián)舯┝鞯奈锢砟M相對較晚。浙江大學[7]與中國空氣動力研究與發(fā)展中心[8]建成了類似的動量沖擊射流裝置。

    從物理模型角度而言,沖擊射流是模擬下?lián)舯┝鞯倪壿嬒嗨颇P?。然而,由于實驗室條件的限制,很難采用沖擊射流對下?lián)舯┝鏖_展大尺度、較快平移速度的試驗研究。即便能夠快速移動,試驗縮尺比也較小,需要較大的試驗空間。目前,國內(nèi)外文獻中最大的噴嘴直徑為4.5 m[9],對于如此大的設(shè)備,模擬下?lián)舯┝鞯钠絼樱ㄋ俣瓤蛇_水平風速的1/3)是非常困難的。

    由于下?lián)舯┝鞒隽鞫蔚钠骄L速剖面符合典型的平面壁面射流特征,Lin 等[10]提出了一種新的物理模型,采用壁面射流裝置研究了二維平面內(nèi)的下?lián)舯┝鞒隽鞫蔚牧鲌鎏卣鳎⒃诔R?guī)邊界層風洞中對下?lián)舯┝鞒隽鞫芜M行了模擬[11],實現(xiàn)了傳統(tǒng)風洞試驗中的幾何縮尺(1∶100~1∶250)。但前期的研究尚未考慮下?lián)舯┝髦欣妆┰频钠絼幼饔?。段旻等[12]通過一個縮聚裝置在邊界層風洞中得到了壁面射流,模擬了下?lián)舯┝髌骄L速剖面;辛亞兵等[13]采用類似裝置研究了穩(wěn)態(tài)下?lián)舯┝髯饔孟碌拇罂邕B續(xù)剛構(gòu)橋風振響應(yīng)。

    現(xiàn)有下?lián)舯┝黠L場模擬試驗裝置大部分僅模擬下?lián)舯┝鞣€(wěn)態(tài)風場,而對下?lián)舯┝鞯姆欠€(wěn)態(tài)風場模擬較少。實際上,空間層面的下?lián)舯┝魇侨S的,產(chǎn)生的渦系是環(huán)狀擴散的,而采用帶協(xié)同流的壁面射流模型僅能模擬二維平面流動情況。為了評估采用二維壁面射流模擬三維下?lián)舯┝鞒隽鞫物L場的準確性,Lin 等[10]對理想下?lián)舯┝鞒隽鞫蔚亩S性假設(shè)進行了驗證,研究表明平面壁面射流與徑向沖擊射流在橫風向的誤差完全可以忽略;即使存在較大的分流角度,陣風前端產(chǎn)生的正面曲率效應(yīng)對下?lián)舯┝髦行木€(移動路徑)流場的影響也非常微弱。總體來看,壁面射流模型與下?lián)舯┝鞯牧鲃訖C制并不完全相同,但在造成結(jié)構(gòu)物破環(huán)的主要區(qū)域—出流段—采用壁面射流模型完全能夠模擬下?lián)舯┝鞒隽鞫物L場的主要特征。

    根據(jù)Joint Airport Weather Studies(JAWS)的實測總結(jié)[14],下?lián)舯┝髌茐牡拇蟾怕蕝^(qū)域主要位于出流段,占整個下?lián)舯┝黠L場的90%以上,如圖1所示(圖中um表示局部最大風速)。因此,在常規(guī)邊界層風洞中采用壁面射流方法進行風洞試驗,可以實現(xiàn)更大的幾何縮尺比。本文在參考已有研究的基礎(chǔ)上,設(shè)計制作了基于大氣邊界層風洞的非穩(wěn)態(tài)風場模擬裝置。在常規(guī)邊界層風洞的基礎(chǔ)上,通過增加壁面射流噴口及風機,在噴口處布置氣動閥門,實現(xiàn)對邊界層風洞中下?lián)舯┝鞒隽鞫畏欠€(wěn)態(tài)風場的模擬。

    圖1 典型微下?lián)舯┝魇疽鈭D[14]Fig.1 Schematic diagram of micro-downburst[14]

    1 試驗概況

    1.1 非穩(wěn)態(tài)風場模擬裝置

    試驗在重慶大學直流式教學風洞中進行。通過改進文獻[10]提出的壁面射流方法并對風洞進行改裝,使直流式邊界層風洞具備了壁面射流模擬功能,如圖2所示。風洞試驗段尺寸為2.4 m×1.8 m×15.0 m(寬×高×長)。壁面射流裝置加裝于試驗段入口位置,通過4 個千斤頂與支架連接,可實現(xiàn)壁面射流裝置的上升與下降:當進行壁面射流試驗時,以千斤頂將壁面射流裝置上升至與試驗段底面平齊;試驗完成后,降下壁面射流裝置,可以常規(guī)使用邊界層風洞。壁面射流裝置分為動力段、過渡段、穩(wěn)流段、回轉(zhuǎn)段和射流噴口。動力段采用3 臺風機并聯(lián)安裝;回轉(zhuǎn)段采用對數(shù)螺旋線設(shè)計,盡量減小由于風場轉(zhuǎn)向?qū)е碌娘L速損耗;壁面射流噴口的高度b=60 mm,寬度與邊界層風洞基本一致,最大出流速度為45 m/s。

    圖2 壁面射流裝置圖Fig.2 Modifications to boundary layer wind tunnel

    大多數(shù)下?lián)舯┝魍瑫r伴隨雷暴云層的平動,這種平動會導(dǎo)致雷暴云層前端下?lián)舯┝黠L場速度增大,增大至峰值后迅速衰減。下?lián)舯┝鞯倪@種時變特征與常規(guī)邊界層穩(wěn)態(tài)風場有著極大的不同,因此,對于平動型下?lián)舯┝?,需要考慮其非穩(wěn)態(tài)風場特性。通過在壁面射流噴口段安裝快開閥門,可以形成具有下?lián)舯┝黠L速剖面的突變風場,實現(xiàn)時變下?lián)舯┝黠L場的試驗?zāi)M;通過閥門控制壁面射流的流量大小,實現(xiàn)噴口風速的非線性脈動沖擊模擬。閥門快開由伸縮氣缸完成,并增加手動節(jié)流閥以調(diào)節(jié)氣缸運動速度,實現(xiàn)閥門開/閉時間可調(diào)(1~10 s 范圍內(nèi)),如圖3所示。試驗風場采用眼鏡蛇三維脈動風速采集系統(tǒng)(Series 100 Cobra Probe)進行采集,采樣頻率為512 Hz。

    圖3 快開閥門裝置示意圖Fig.3 Schematic diagram of quick opening valve device

    2 風場模擬結(jié)果

    2.1 穩(wěn)態(tài)風場結(jié)果

    當閥門保持全開、噴口出流速度為28.30 m/s 時,壁面射流不同順流向位置的平均風速剖面如圖4所示,湍流度剖面如圖5所示。圖中,x為順流向位置,y為豎直高度方向位置(坐標系如圖2所示)。可以看出,在測量高度范圍內(nèi),基本沒有回流的出現(xiàn),隨著順流向距離的增大,風速剖面下部的速度逐漸減小,而上部速度逐漸增大,轉(zhuǎn)折區(qū)域約在200~300 mm范圍內(nèi)。壁面射流湍流度剖面則呈現(xiàn)出明顯的雙峰特征(即內(nèi)層近壁面峰值與外層峰值),且隨著順流向距離的增大,外層峰值出現(xiàn)的豎向高度位置不斷上移。

    圖4 不同順流向位置平均風速剖面Fig.4 Mean velocity profiles at various streamwise locations

    圖5 不同順流向位置湍流度剖面Fig.5 Turbulent intensity profiles at various streamwise locations

    為得到壁面射流的自相似剖面,通常采用特征長度與特征速度對其進行無量綱處理[15]。壁面射流的無量綱速度剖面如圖6所示??梢钥闯觯诿嫔淞黠L洞試驗結(jié)果表現(xiàn)出較好的自相似性,僅在x=20b處有一定偏差,這是由于此位置壁面射流轉(zhuǎn)捩沒有完成,尚未進入完全發(fā)展階段。風洞試驗與Eriksson等[16]的試驗結(jié)果、三種典型下?lián)舯┝黠L速剖面經(jīng)驗?zāi)P蚚17-19]非常吻合。

    圖6 無量綱平均風速剖面Fig.6 Normalized mean velocity profiles

    2.2 非穩(wěn)態(tài)風場結(jié)果

    實際上,下?lián)舯┝黠L場中各位置的風速與距離下?lián)舯┝鳑_擊地面點的相對位置、下?lián)舯┝骼妆┰频钠絼铀俣?、下沉氣流的出流速度和高度以及下?lián)舯┝鞯某隽鲝姸扔嘘P(guān),這些因素同樣具有較強的隨機性。因此,本文僅通過控制快開閥門的轉(zhuǎn)動速度和轉(zhuǎn)動角度,重點模擬下?lián)舯┝麟S時間的變化及其豎向風速剖面,未考慮同一下?lián)舯┝魇录胁煌恢玫娘L速分布特征。

    通過手動節(jié)流閥控制伸縮氣缸氣壓大小,調(diào)整閥門開關(guān)時間。如圖7所示,瞬態(tài)風場采用2 個眼鏡蛇探頭同時進行測量,測量的豎向位置分別為30、50、80、120、200、300 和500 mm,一個探頭固定于30 mm處,另一個探頭則上下移動。在測量不同高度的風速時,數(shù)據(jù)采集時間并不同步,因此,通過統(tǒng)一最大風速對應(yīng)的時間點,采用固定探頭對不同高度的風速進行時間同步,從而更加合理地分析豎向風場特征[20]。下?lián)舯┝鞯膶崪y記錄較少,本文選擇文獻中可用的安德魯斯空軍基地(Andrews AFB)下?lián)舯┝鲗崪y記錄作為目標事件[1],從而得到典型的下?lián)舯┝黠L速時程。

    圖7 非穩(wěn)態(tài)風場測量的探頭布置Fig.7 Probe arrangement for unsteady wind field measurement

    任意高度處下?lián)舯┝鞯捻樍飨蝻L速可以分解為確定的時變平均分量V(z,t)和 脈動分量v′(z,t)[21]:

    其中,平均風速隨時間連續(xù)變化且能夠反映原始風速的變化趨勢。Su 等[22]采用不同的方法研究時變平均風速與脈動風速的分離,建議采用高階的離散小波變換(DWT)與總體經(jīng)驗?zāi)B(tài)分解(EEMD)方法提取時變平均風速,就能夠較好地反映原始風速的變化趨勢。因此,本文采取高階的DWT 方法(采用10 階小波)進行時變平均風速的分離。當閥門開啟時間為3 s、壁面射流出流速度為33.06 m/s 時,x=100b處移動探頭在不同高度時(50、80 和120 mm)固定探頭的時變平均風速如圖8所示。圖例中的2 個高度分別表示固定探頭和移動探頭所在的豎向高度位置。

    圖8 不同測量高度風速同步Fig.8 Wind speed synchronization of fixed probe data (unsynchronized and synchronized)

    可以看出,由壁面射流裝置及閥門開關(guān)產(chǎn)生的出流速度在每次試驗中都存在一定差異,這主要是由電機啟動扭矩、制動性能以及手動節(jié)流閥中氣壓大小的不確定性所導(dǎo)致。因此,在測量風速時,移動探頭在每個高度處進行6 次重復(fù)試驗,挑選較為接近的幾組數(shù)據(jù)中的一組數(shù)據(jù)作為最終結(jié)果。壁面射流出流速度為33.06 m/s、主風洞產(chǎn)出的環(huán)境風速為6.60 m/s 時,x=100b處不同高度的風速時程以及時變平均風速和脈動風速的分離如圖9所示。

    從圖9可以看出,不同高度處的風速時程都呈現(xiàn)出2 個峰值的特征,且隨著高度增大,最大平均風速減小,說明采用本文方法能有效模擬典型下?lián)舯┝鞒隽鞫物L速的時變特征。最大時變平均風速豎向剖面如圖10所示(圖中陰影區(qū)域為Hjelmfelt[14]統(tǒng)計的實測數(shù)據(jù)區(qū)間),可以看出:豎向風速剖面表現(xiàn)出下?lián)舯┝鞯湫惋L速剖面的基本特征,與Hjelmfelt 的實測數(shù)據(jù)較為吻合;但由于測點較少,且僅是取時變平均風速的最大值,曲線并不非常平滑。

    圖9 不同高度同步后的風速時程Fig.9 Synchronized velocity time histories at various heights

    圖10 最大時變平均風速豎向剖面Fig.10 Vertical profile of maximum time-varying mean velocity

    下?lián)舯┝髯畲箫L速發(fā)生高度約為0.02~0.05D(D為下?lián)舯┝鞒隽髦睆剑23-24]。從圖4和10 可以看出,順流向x=100b處的最大風速大約發(fā)生在豎向高度80 mm 處,可以估計出該位置的等效直徑約在2.67~1.60 m 范圍內(nèi)。Abd-Elaal 等通過大量參數(shù)分析,估計出Andrews AFB 下?lián)舯┝鞯某隽髦睆郊s為660 m[24]。因此,本文的壁面射流裝置在x=100b處模擬Andrews AFB 下?lián)舯┝鞯目s尺比約為1∶247~1∶412。出流速度為33.06 m/s 時,最大時變平均風速為15.70 m/s,Andrews AFB 下?lián)舯┝髯畲笏矔r風速約為66.82 m/s[1],對其移動平均后得到的最大風速為54.21 m/s,如圖11所示,可以得到速度縮尺比為1∶3.67,需要模擬的時間縮尺比為1∶67.30~1∶112.26。從圖11可以看出,風速從開始增大直至峰值所需時間約為105 s,因此需要閥門打開的時間為0.94~1.56 s。取閥門打開時間為1.50 s 進行模擬,模擬結(jié)果如圖12所示。從圖中可以看出:風洞試驗結(jié)果和實際的下?lián)舯┝骶哂休^好的吻合度,特別是在第一個峰值處,而在第二個峰值處較實際風場偏大。

    圖11 Andrews AFB 下?lián)舯┝饕苿悠骄L速時程[1]Fig.11 Moving-averaged wind speed history of Andrews AFB downburst[1]

    圖12 風洞試驗結(jié)果與目標下?lián)舯┝黠L場時變平均風速對比Fig.12 Comparison of simulated and recorded time-varying mean velocity

    不同順流向位置的最大風速高度基本不變,因此,通過調(diào)節(jié)壁面射流的出流速度控制速度縮尺比,根據(jù)閥門開啟的速度控制時間縮尺比,能夠有效得到目標風場的時變風速。

    順流向x=100b處的湍流度豎向剖面如圖13所示。由于缺乏下?lián)舯┝鞯膶崪y湍流度豎向剖面,本文采用中國《建筑結(jié)構(gòu)荷載規(guī)范》(GB 50009—2012)[25]中的大氣邊界層風場與試驗結(jié)果進行對比。下?lián)舯┝黠L速具有時變特征,因此湍流度的計算并沒有固定方法。本文湍流度采用特征時間內(nèi)脈動風速的均方根值與最大平均風速的比值。其中,定義特征時間為下?lián)舯┝鳑_擊地面后、風速從開始急劇增大至衰減為最大風速一半所需的時間[26]。

    圖13 湍流度豎向剖面Fig.13 Vertical profiles of turbulent intensity

    從圖13可以看出,風洞試驗結(jié)果在近地面與規(guī)范值較為一致,當豎向高度超過最大風速所在高度后(0.05D,30 m),湍流度迅速增大,這是由于壁面射流風場與邊界層風場發(fā)生混合,渦旋的產(chǎn)生與耗散加速,這與穩(wěn)態(tài)壁面射流試驗中得到的湍流度剖面非常類似。

    風洞試驗x=100b、豎向高度80 mm 處的特征時間內(nèi)脈動風速功率譜如圖14所示,圖中同時給出von Karman 譜進行對比。根據(jù)Kolmogorov 湍流能譜假設(shè),在慣性子區(qū)內(nèi),湍流的能量將由最小尺度的渦耗散,這種渦的尺度稱為Kolmogorov 長度尺度,而慣性子區(qū)內(nèi)的湍流能譜滿足“?5/3 定律”。從圖中可以看出,特征時間內(nèi)脈動風速功率譜與30 m 處的von Karman 譜在5 Hz 之后吻合,并且符合“?5/3 定律”。

    圖14 80 mm 高度處特征時間內(nèi)脈動風速功率譜Fig.14 PSD of residual fluctuating velocity during characteristic time at 80 mm height

    從上述分析可以看出,采用壁面射流裝置能夠較好地實現(xiàn)下?lián)舯┝鞒隽鞫物L場的模擬,主要包括豎向風速剖面以及時變風速時程。對于湍流度,由于下?lián)舯┝鞯耐牧鞫热狈崪y資料,無法進行對比研究,本文僅研究了光滑壁面一種情況,尚未考慮粗糙元的設(shè)置。總體而言,采用壁面射流裝置模擬下?lián)舯┝鞣欠€(wěn)態(tài)風場是較為有效的。

    3 結(jié) 論

    基于平面壁面射流模型,在邊界層風洞中模擬了穩(wěn)態(tài)及非穩(wěn)態(tài)下?lián)舯┝鞒隽鞫物L場,結(jié)果表明:

    1)穩(wěn)態(tài)壁面射流湍流度剖面呈現(xiàn)出明顯的雙峰特征(內(nèi)層近壁面峰值與外層峰值),且隨著順流向距離的增大,外層峰值出現(xiàn)的豎向高度位置不斷上移;而隨著射流出流速度的增大,湍流度呈增大趨勢,外層峰值出現(xiàn)的豎向高度位置也上移。

    2)與中國規(guī)范湍流度相比,非穩(wěn)態(tài)壁面射流湍流度剖面在最大風速高度以下較為接近,而在最大風速高度之上,非穩(wěn)態(tài)壁面射流湍流度較規(guī)范值偏大。

    3)采用壁面射流裝置產(chǎn)生的下?lián)舯┝鞣欠€(wěn)態(tài)出流段風速功率譜與von Karman 譜較為吻合,且符合“–5/3 定律”。

    4)基于平面壁面射流方法的下?lián)舯┝髂M試驗裝置在邊界層風洞中可實現(xiàn)非穩(wěn)態(tài)風場模擬,為下?lián)舯┝魉斤L場作用下的結(jié)構(gòu)響應(yīng)研究打下基礎(chǔ)。

    猜你喜歡
    暴流時變風場
    江蘇典型下?lián)舯┝黠L暴結(jié)構(gòu)特征統(tǒng)計分析*
    氣象學報(2022年4期)2022-09-06 01:59:50
    基于FLUENT的下?lián)舯┝魅S風場建模
    基于壁面射流的下?lián)舯┝鞣欠€(wěn)態(tài)風場大渦模擬
    一種改進的多普勒天氣雷達下?lián)舯┝髯R別算法
    “最美風場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    基于時變Copula的股票市場相關(guān)性分析
    智富時代(2017年4期)2017-04-27 17:08:47
    煙氣輪機復(fù)合故障時變退化特征提取
    側(cè)向風場中無人機的飛行研究
    基于MEP法的在役橋梁時變可靠度研究
    風場條件下LPG 瞬時泄漏擴散的數(shù)值模擬
    久久人妻福利社区极品人妻图片| 亚洲精品中文字幕在线视频| 欧美国产精品一级二级三级| 久久精品亚洲av国产电影网| 亚洲性夜色夜夜综合| 女人久久www免费人成看片| 久久精品aⅴ一区二区三区四区| 精品视频人人做人人爽| 久久中文看片网| www.999成人在线观看| 性高湖久久久久久久久免费观看| 视频区欧美日本亚洲| www.999成人在线观看| av有码第一页| 精品一品国产午夜福利视频| 涩涩av久久男人的天堂| 搡老岳熟女国产| 午夜免费成人在线视频| 韩国精品一区二区三区| 日韩制服丝袜自拍偷拍| 成人国产一区最新在线观看| 欧美+亚洲+日韩+国产| 99久久人妻综合| 国产在线一区二区三区精| 精品亚洲成a人片在线观看| 日本91视频免费播放| 91精品三级在线观看| 日韩视频一区二区在线观看| 欧美日本中文国产一区发布| 亚洲精品乱久久久久久| 国产成人av教育| 两性夫妻黄色片| 国产精品免费视频内射| 一级毛片女人18水好多| 国产片内射在线| 国产日韩一区二区三区精品不卡| 国产一区二区 视频在线| 亚洲情色 制服丝袜| 亚洲中文日韩欧美视频| 国产老妇伦熟女老妇高清| 亚洲三区欧美一区| 97精品久久久久久久久久精品| 麻豆乱淫一区二区| 国产精品一二三区在线看| 免费在线观看影片大全网站| 亚洲中文字幕日韩| 国产成人欧美| 久久人人97超碰香蕉20202| 女人高潮潮喷娇喘18禁视频| 亚洲第一av免费看| 欧美老熟妇乱子伦牲交| 国产伦理片在线播放av一区| 捣出白浆h1v1| 午夜精品久久久久久毛片777| 精品久久久久久电影网| 国产男人的电影天堂91| 久久热在线av| 欧美乱码精品一区二区三区| www日本在线高清视频| 丝袜喷水一区| 欧美在线一区亚洲| 搡老乐熟女国产| 一级a爱视频在线免费观看| 国产精品久久久av美女十八| 美女大奶头黄色视频| 美女扒开内裤让男人捅视频| 一二三四在线观看免费中文在| 男女无遮挡免费网站观看| 精品免费久久久久久久清纯 | 人人妻人人澡人人爽人人夜夜| 最黄视频免费看| 久久99热这里只频精品6学生| 美女视频免费永久观看网站| 老汉色av国产亚洲站长工具| 性少妇av在线| 一级片'在线观看视频| 各种免费的搞黄视频| 亚洲精品久久久久久婷婷小说| 后天国语完整版免费观看| 一区二区av电影网| 黑人猛操日本美女一级片| 日日爽夜夜爽网站| 精品福利观看| 黑人巨大精品欧美一区二区蜜桃| avwww免费| 一二三四社区在线视频社区8| 免费观看人在逋| 汤姆久久久久久久影院中文字幕| 国产一区二区在线观看av| 大片电影免费在线观看免费| 制服人妻中文乱码| 久久国产亚洲av麻豆专区| 夜夜骑夜夜射夜夜干| 99国产精品99久久久久| 国产1区2区3区精品| 丝袜美足系列| 老司机在亚洲福利影院| 老熟妇仑乱视频hdxx| 国产精品香港三级国产av潘金莲| 女性生殖器流出的白浆| av有码第一页| 中亚洲国语对白在线视频| 国产熟女午夜一区二区三区| 国产男人的电影天堂91| 天堂8中文在线网| 各种免费的搞黄视频| 亚洲欧洲日产国产| 91麻豆av在线| 脱女人内裤的视频| 手机成人av网站| 国产精品香港三级国产av潘金莲| 如日韩欧美国产精品一区二区三区| 精品一区二区三区av网在线观看 | 侵犯人妻中文字幕一二三四区| 日韩制服丝袜自拍偷拍| 国产亚洲av片在线观看秒播厂| 国产成人系列免费观看| 久久久欧美国产精品| 国产精品一区二区在线观看99| 满18在线观看网站| 69精品国产乱码久久久| 成人三级做爰电影| 午夜福利视频精品| 在线观看www视频免费| 亚洲av成人不卡在线观看播放网 | 国产一区二区在线观看av| 十分钟在线观看高清视频www| 国产精品欧美亚洲77777| 50天的宝宝边吃奶边哭怎么回事| 热re99久久国产66热| 九色亚洲精品在线播放| 少妇粗大呻吟视频| 新久久久久国产一级毛片| 欧美午夜高清在线| 女人久久www免费人成看片| 国产黄频视频在线观看| av片东京热男人的天堂| 亚洲国产精品一区三区| 欧美人与性动交α欧美软件| 欧美亚洲日本最大视频资源| 搡老岳熟女国产| 美女高潮到喷水免费观看| 国产成人av教育| 欧美性长视频在线观看| 亚洲精品美女久久av网站| 午夜福利在线观看吧| 精品久久久久久电影网| 亚洲国产欧美在线一区| 久久热在线av| 97人妻天天添夜夜摸| 91字幕亚洲| 欧美黑人欧美精品刺激| 欧美xxⅹ黑人| 精品国产一区二区久久| 一区二区三区四区激情视频| 婷婷色av中文字幕| 美女中出高潮动态图| 老汉色av国产亚洲站长工具| 俄罗斯特黄特色一大片| 一区二区三区精品91| 午夜老司机福利片| www.熟女人妻精品国产| 人妻人人澡人人爽人人| 少妇粗大呻吟视频| 亚洲视频免费观看视频| 久久这里只有精品19| 国产精品一区二区免费欧美 | 国产亚洲精品久久久久5区| 男人操女人黄网站| 午夜成年电影在线免费观看| 美国免费a级毛片| 欧美 日韩 精品 国产| 考比视频在线观看| 纵有疾风起免费观看全集完整版| 三级毛片av免费| 高潮久久久久久久久久久不卡| 亚洲少妇的诱惑av| 桃红色精品国产亚洲av| 夜夜骑夜夜射夜夜干| 国产91精品成人一区二区三区 | 窝窝影院91人妻| 日本91视频免费播放| 久久久久国内视频| 搡老熟女国产l中国老女人| 国产精品久久久人人做人人爽| 国产av又大| 午夜日韩欧美国产| 蜜桃在线观看..| 亚洲人成电影免费在线| 亚洲国产欧美网| 999精品在线视频| 久久精品熟女亚洲av麻豆精品| 国产伦理片在线播放av一区| 日韩三级视频一区二区三区| 99国产精品一区二区蜜桃av | 亚洲av成人不卡在线观看播放网 | 久久人人爽av亚洲精品天堂| 久久久久国产精品人妻一区二区| 涩涩av久久男人的天堂| 美女高潮喷水抽搐中文字幕| 少妇裸体淫交视频免费看高清 | a级毛片黄视频| 91成人精品电影| 黄色视频,在线免费观看| 欧美成狂野欧美在线观看| 欧美少妇被猛烈插入视频| 在线十欧美十亚洲十日本专区| 亚洲精品一区蜜桃| 黄色视频不卡| 久久精品成人免费网站| 成人三级做爰电影| 老司机深夜福利视频在线观看 | 国产区一区二久久| 99国产精品一区二区三区| 国产野战对白在线观看| 国产av又大| 色婷婷久久久亚洲欧美| 欧美日韩中文字幕国产精品一区二区三区 | 啦啦啦 在线观看视频| 成人国语在线视频| 欧美 日韩 精品 国产| 男女下面插进去视频免费观看| 天堂8中文在线网| 搡老岳熟女国产| 最黄视频免费看| 在线十欧美十亚洲十日本专区| 老汉色av国产亚洲站长工具| 正在播放国产对白刺激| 成人av一区二区三区在线看 | 女警被强在线播放| 成人三级做爰电影| 欧美精品啪啪一区二区三区 | 午夜福利在线免费观看网站| 国产熟女午夜一区二区三区| 我的亚洲天堂| 美女视频免费永久观看网站| www.av在线官网国产| 精品国产乱码久久久久久小说| 91精品伊人久久大香线蕉| 12—13女人毛片做爰片一| 中文字幕人妻丝袜一区二区| 91国产中文字幕| 亚洲国产精品一区三区| 色婷婷久久久亚洲欧美| 乱人伦中国视频| 亚洲精品粉嫩美女一区| 欧美精品一区二区免费开放| 久久影院123| 久久久久久免费高清国产稀缺| 亚洲久久久国产精品| 亚洲人成电影免费在线| 国产成人精品久久二区二区免费| www.自偷自拍.com| 视频区欧美日本亚洲| 久久久久久亚洲精品国产蜜桃av| 国产成+人综合+亚洲专区| 亚洲第一av免费看| 国产精品久久久久久人妻精品电影 | av超薄肉色丝袜交足视频| 淫妇啪啪啪对白视频 | 99精国产麻豆久久婷婷| 亚洲国产av新网站| 搡老乐熟女国产| 丰满少妇做爰视频| av天堂在线播放| 狠狠狠狠99中文字幕| 精品一区二区三区四区五区乱码| 日日夜夜操网爽| 自线自在国产av| 丝袜美腿诱惑在线| 亚洲精品国产一区二区精华液| 日韩一区二区三区影片| a级片在线免费高清观看视频| 日韩人妻精品一区2区三区| av一本久久久久| 国产亚洲一区二区精品| 国产又色又爽无遮挡免| 高清黄色对白视频在线免费看| svipshipincom国产片| 亚洲精品一卡2卡三卡4卡5卡 | 国产黄色免费在线视频| 伦理电影免费视频| 欧美 日韩 精品 国产| 俄罗斯特黄特色一大片| 精品欧美一区二区三区在线| 秋霞在线观看毛片| 国产一区二区激情短视频 | 日韩制服丝袜自拍偷拍| 精品免费久久久久久久清纯 | 欧美日韩中文字幕国产精品一区二区三区 | 大型av网站在线播放| 啦啦啦中文免费视频观看日本| 国产不卡av网站在线观看| 在线av久久热| 日本91视频免费播放| 人人妻人人澡人人看| 国产精品一二三区在线看| 亚洲黑人精品在线| 欧美国产精品va在线观看不卡| 久久女婷五月综合色啪小说| 热re99久久国产66热| 老汉色av国产亚洲站长工具| 黄色怎么调成土黄色| 人妻 亚洲 视频| 最新在线观看一区二区三区| 涩涩av久久男人的天堂| 夜夜夜夜夜久久久久| 天天躁夜夜躁狠狠躁躁| 国产极品粉嫩免费观看在线| 欧美 日韩 精品 国产| a在线观看视频网站| a 毛片基地| 亚洲欧美一区二区三区黑人| 欧美性长视频在线观看| 亚洲七黄色美女视频| 亚洲国产欧美网| 19禁男女啪啪无遮挡网站| 丰满人妻熟妇乱又伦精品不卡| 国产真人三级小视频在线观看| 女人精品久久久久毛片| 日本91视频免费播放| 国产成人一区二区三区免费视频网站| 在线亚洲精品国产二区图片欧美| 桃红色精品国产亚洲av| 激情视频va一区二区三区| 老鸭窝网址在线观看| 亚洲自偷自拍图片 自拍| 午夜福利在线免费观看网站| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| 国产av精品麻豆| 一区在线观看完整版| 国产亚洲精品一区二区www | 黑人巨大精品欧美一区二区mp4| 久久热在线av| 丝袜在线中文字幕| 国产男女超爽视频在线观看| 99香蕉大伊视频| 王馨瑶露胸无遮挡在线观看| 欧美av亚洲av综合av国产av| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美网| 考比视频在线观看| 人人妻人人澡人人看| 精品福利观看| 欧美黑人欧美精品刺激| 首页视频小说图片口味搜索| 国产精品国产三级国产专区5o| 亚洲精品美女久久av网站| 十分钟在线观看高清视频www| 久久国产精品影院| 一个人免费在线观看的高清视频 | 欧美精品亚洲一区二区| 五月开心婷婷网| 欧美中文综合在线视频| 精品亚洲成国产av| 一区二区av电影网| 午夜免费成人在线视频| 一本色道久久久久久精品综合| 欧美+亚洲+日韩+国产| 人人妻人人澡人人看| 天天躁夜夜躁狠狠躁躁| 少妇人妻久久综合中文| 少妇精品久久久久久久| 99热网站在线观看| 妹子高潮喷水视频| 少妇裸体淫交视频免费看高清 | 日韩欧美一区二区三区在线观看 | 老司机深夜福利视频在线观看 | 成人18禁高潮啪啪吃奶动态图| 女人高潮潮喷娇喘18禁视频| 一级,二级,三级黄色视频| 啦啦啦在线免费观看视频4| 午夜福利视频在线观看免费| 黄片大片在线免费观看| 99久久综合免费| 男人添女人高潮全过程视频| 免费看十八禁软件| 久久国产亚洲av麻豆专区| 国产淫语在线视频| av福利片在线| 男人爽女人下面视频在线观看| 十八禁网站免费在线| 久久久久视频综合| 妹子高潮喷水视频| 亚洲七黄色美女视频| 国精品久久久久久国模美| 人妻 亚洲 视频| 两人在一起打扑克的视频| 性色av一级| 免费女性裸体啪啪无遮挡网站| av在线播放精品| 午夜福利一区二区在线看| 国产一区二区在线观看av| 韩国高清视频一区二区三区| 99国产精品一区二区蜜桃av | 午夜老司机福利片| 在线观看人妻少妇| 欧美日韩视频精品一区| 下体分泌物呈黄色| 精品一品国产午夜福利视频| 亚洲天堂av无毛| 欧美 亚洲 国产 日韩一| 久久亚洲国产成人精品v| 午夜老司机福利片| 欧美性长视频在线观看| 精品视频人人做人人爽| 嫁个100分男人电影在线观看| 欧美性长视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品国产综合久久久| 日韩 欧美 亚洲 中文字幕| 国产欧美日韩一区二区精品| 亚洲成人手机| 蜜桃在线观看..| 一本—道久久a久久精品蜜桃钙片| 国产精品一二三区在线看| 亚洲欧美一区二区三区黑人| 91精品国产国语对白视频| 黄片小视频在线播放| 亚洲avbb在线观看| 日韩欧美国产一区二区入口| 汤姆久久久久久久影院中文字幕| 国产精品熟女久久久久浪| 美女高潮到喷水免费观看| 国产日韩一区二区三区精品不卡| 国产淫语在线视频| 夜夜夜夜夜久久久久| 成人av一区二区三区在线看 | 9热在线视频观看99| 久久中文字幕一级| 窝窝影院91人妻| 欧美黄色片欧美黄色片| 国产激情久久老熟女| 在线观看一区二区三区激情| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产区一区二| 满18在线观看网站| 啦啦啦在线免费观看视频4| 日韩欧美一区视频在线观看| 日韩制服骚丝袜av| 久久人人爽人人片av| 69av精品久久久久久 | 久久国产精品人妻蜜桃| 国产精品久久久av美女十八| 三级毛片av免费| 午夜视频精品福利| a在线观看视频网站| 老司机在亚洲福利影院| 久久人妻熟女aⅴ| 精品国产乱子伦一区二区三区 | 亚洲精品国产一区二区精华液| 亚洲av成人不卡在线观看播放网 | 国产不卡av网站在线观看| 国产熟女午夜一区二区三区| 精品福利观看| 美女脱内裤让男人舔精品视频| 成人18禁高潮啪啪吃奶动态图| av欧美777| 日韩中文字幕视频在线看片| 日本五十路高清| 男人操女人黄网站| 一区二区三区激情视频| 香蕉丝袜av| 欧美亚洲 丝袜 人妻 在线| 国产精品99久久99久久久不卡| 欧美97在线视频| 2018国产大陆天天弄谢| 精品人妻一区二区三区麻豆| 最近最新免费中文字幕在线| 又紧又爽又黄一区二区| 国产成人精品无人区| 国产精品久久久久成人av| 免费女性裸体啪啪无遮挡网站| netflix在线观看网站| 亚洲精品粉嫩美女一区| 91精品国产国语对白视频| av有码第一页| 精品国内亚洲2022精品成人 | e午夜精品久久久久久久| 乱人伦中国视频| 久久这里只有精品19| av线在线观看网站| 午夜老司机福利片| 91字幕亚洲| 最黄视频免费看| 欧美 日韩 精品 国产| 男人操女人黄网站| 超碰97精品在线观看| 在线观看免费日韩欧美大片| 亚洲欧美精品综合一区二区三区| 天天操日日干夜夜撸| 老司机影院成人| 男女之事视频高清在线观看| 国产免费一区二区三区四区乱码| 国产一区二区激情短视频 | 99国产精品免费福利视频| 欧美精品av麻豆av| 国产av一区二区精品久久| 法律面前人人平等表现在哪些方面 | 亚洲人成77777在线视频| xxxhd国产人妻xxx| 午夜福利,免费看| 成人18禁高潮啪啪吃奶动态图| 久久天堂一区二区三区四区| 色综合欧美亚洲国产小说| 亚洲成人免费av在线播放| 热99re8久久精品国产| 亚洲精品国产av蜜桃| 美女主播在线视频| 免费在线观看影片大全网站| 男人操女人黄网站| 悠悠久久av| 精品国产乱码久久久久久小说| 三上悠亚av全集在线观看| 汤姆久久久久久久影院中文字幕| 一级a爱视频在线免费观看| a级毛片在线看网站| 亚洲专区国产一区二区| 这个男人来自地球电影免费观看| 一二三四在线观看免费中文在| 久久免费观看电影| 亚洲欧美日韩另类电影网站| 亚洲精品日韩在线中文字幕| 国产精品国产三级国产专区5o| 亚洲国产欧美一区二区综合| a级毛片黄视频| 人成视频在线观看免费观看| 国产精品久久久久成人av| 国产免费福利视频在线观看| 一区二区三区乱码不卡18| 免费高清在线观看日韩| 天天躁日日躁夜夜躁夜夜| 中文字幕另类日韩欧美亚洲嫩草| 国产免费一区二区三区四区乱码| 国产精品自产拍在线观看55亚洲 | 男女之事视频高清在线观看| 午夜福利影视在线免费观看| 中文字幕人妻丝袜制服| 日韩免费高清中文字幕av| 国产亚洲精品久久久久5区| 免费在线观看视频国产中文字幕亚洲 | 操出白浆在线播放| 亚洲av成人不卡在线观看播放网 | 亚洲全国av大片| 国产麻豆69| 自线自在国产av| e午夜精品久久久久久久| 一个人免费看片子| 精品欧美一区二区三区在线| 精品国产一区二区三区四区第35| 亚洲国产欧美在线一区| 欧美黑人精品巨大| 制服人妻中文乱码| 国产精品国产三级国产专区5o| 日本a在线网址| 免费在线观看视频国产中文字幕亚洲 | 一进一出抽搐动态| 国产欧美日韩一区二区三区在线| 国产国语露脸激情在线看| 自线自在国产av| 亚洲av男天堂| 免费少妇av软件| 首页视频小说图片口味搜索| 久久久国产欧美日韩av| 人妻久久中文字幕网| 久久亚洲精品不卡| 男女国产视频网站| 亚洲视频免费观看视频| 两个人看的免费小视频| 欧美黑人欧美精品刺激| 看免费av毛片| 午夜福利免费观看在线| 日韩免费高清中文字幕av| 亚洲一码二码三码区别大吗| 国产精品免费视频内射| 嫁个100分男人电影在线观看| 悠悠久久av| 色婷婷av一区二区三区视频| 亚洲第一av免费看| 正在播放国产对白刺激| 一二三四社区在线视频社区8| 久久国产亚洲av麻豆专区| 久久中文字幕一级| 热99国产精品久久久久久7| 满18在线观看网站| 黑人巨大精品欧美一区二区mp4| 最新的欧美精品一区二区| 激情视频va一区二区三区| 欧美日韩视频精品一区| 自线自在国产av| 少妇的丰满在线观看| 午夜福利乱码中文字幕| 黄色a级毛片大全视频| 午夜福利在线免费观看网站| 欧美国产精品va在线观看不卡| 中国美女看黄片| 亚洲精品在线美女| bbb黄色大片| 巨乳人妻的诱惑在线观看| 亚洲欧洲精品一区二区精品久久久| 91精品三级在线观看| 99九九在线精品视频| 亚洲欧美精品自产自拍| 亚洲欧美激情在线| 日韩有码中文字幕| 日韩欧美国产一区二区入口| av又黄又爽大尺度在线免费看| 国产精品国产av在线观看| 黄色怎么调成土黄色| 国产av精品麻豆| 国产精品99久久99久久久不卡| 超色免费av| 一区二区三区激情视频| 国产精品秋霞免费鲁丝片|