• <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ù)值模擬
    女性被躁到高潮视频| 水蜜桃什么品种好| 99热全是精品| 久久狼人影院| 国产av精品麻豆| videosex国产| 久久精品亚洲av国产电影网| 国产精品 欧美亚洲| 人妻人人澡人人爽人人| 十八禁高潮呻吟视频| 看免费成人av毛片| 日本色播在线视频| 久久久久久久精品精品| 亚洲精品第二区| 十八禁高潮呻吟视频| 女人久久www免费人成看片| 欧美日韩成人在线一区二区| 69精品国产乱码久久久| 99热国产这里只有精品6| 日韩欧美一区视频在线观看| 精品一区二区三卡| 国产午夜精品一二区理论片| 国产伦理片在线播放av一区| av.在线天堂| 国产男女超爽视频在线观看| 日韩人妻精品一区2区三区| 国产伦人伦偷精品视频| 汤姆久久久久久久影院中文字幕| 777久久人妻少妇嫩草av网站| 欧美日韩精品网址| 热re99久久精品国产66热6| 一区二区三区激情视频| 久久精品国产a三级三级三级| 国产精品久久久av美女十八| 久久毛片免费看一区二区三区| 99国产精品免费福利视频| 人人妻人人澡人人爽人人夜夜| 男女之事视频高清在线观看 | 欧美日韩国产mv在线观看视频| 国产 精品1| 男女无遮挡免费网站观看| 久久久国产精品麻豆| 最黄视频免费看| 日韩成人av中文字幕在线观看| 妹子高潮喷水视频| 欧美黑人欧美精品刺激| 欧美精品一区二区免费开放| 黄色怎么调成土黄色| 久久久久久免费高清国产稀缺| 久久天堂一区二区三区四区| 麻豆av在线久日| 亚洲欧洲日产国产| 麻豆乱淫一区二区| 亚洲视频免费观看视频| 色94色欧美一区二区| 亚洲国产欧美一区二区综合| 欧美日韩国产mv在线观看视频| 亚洲婷婷狠狠爱综合网| 人妻人人澡人人爽人人| 亚洲精品成人av观看孕妇| 午夜影院在线不卡| 大香蕉久久网| 欧美亚洲 丝袜 人妻 在线| 中国三级夫妇交换| 久久久久人妻精品一区果冻| 女性被躁到高潮视频| 在线观看一区二区三区激情| 久久女婷五月综合色啪小说| 又大又爽又粗| 日本猛色少妇xxxxx猛交久久| 国产在视频线精品| 免费看av在线观看网站| 免费看不卡的av| 亚洲七黄色美女视频| 婷婷色综合大香蕉| 国产精品三级大全| 精品少妇内射三级| 一二三四在线观看免费中文在| 免费黄频网站在线观看国产| 男女床上黄色一级片免费看| 黑丝袜美女国产一区| 亚洲一区二区三区欧美精品| 亚洲美女视频黄频| 大码成人一级视频| 在线观看国产h片| 亚洲欧美一区二区三区久久| 亚洲欧美日韩另类电影网站| 亚洲图色成人| 免费在线观看完整版高清| 叶爱在线成人免费视频播放| 亚洲国产欧美一区二区综合| 99久久精品国产亚洲精品| 中国国产av一级| 亚洲综合精品二区| 亚洲精品国产区一区二| av视频免费观看在线观看| 99国产综合亚洲精品| 欧美精品一区二区免费开放| 国产男人的电影天堂91| 交换朋友夫妻互换小说| 欧美黑人欧美精品刺激| 国产成人精品在线电影| 国产成人精品无人区| 两性夫妻黄色片| 精品一区在线观看国产| 国产在视频线精品| 美女视频免费永久观看网站| 日韩视频在线欧美| 黑丝袜美女国产一区| 日本色播在线视频| 成人影院久久| 不卡av一区二区三区| 午夜精品国产一区二区电影| 国产av一区二区精品久久| 久久久久精品人妻al黑| 亚洲人成电影观看| 欧美中文综合在线视频| 亚洲欧美精品综合一区二区三区| 久久av网站| 美女主播在线视频| 免费在线观看视频国产中文字幕亚洲 | 又大又黄又爽视频免费| 热99国产精品久久久久久7| 日韩 欧美 亚洲 中文字幕| 国产女主播在线喷水免费视频网站| 男女高潮啪啪啪动态图| 中文字幕制服av| 夫妻性生交免费视频一级片| 熟女少妇亚洲综合色aaa.| 麻豆乱淫一区二区| 久久国产精品大桥未久av| 国产男人的电影天堂91| 天堂俺去俺来也www色官网| 啦啦啦 在线观看视频| 午夜老司机福利片| 欧美 亚洲 国产 日韩一| 午夜福利一区二区在线看| 精品福利永久在线观看| 成人18禁高潮啪啪吃奶动态图| 飞空精品影院首页| 黄频高清免费视频| 久久综合国产亚洲精品| 十八禁高潮呻吟视频| 波多野结衣一区麻豆| 中文字幕高清在线视频| 久久精品熟女亚洲av麻豆精品| 又大又爽又粗| av国产久精品久网站免费入址| 亚洲国产看品久久| 欧美人与性动交α欧美精品济南到| 丰满少妇做爰视频| 热99久久久久精品小说推荐| 亚洲一卡2卡3卡4卡5卡精品中文| 国产老妇伦熟女老妇高清| 国产又色又爽无遮挡免| 亚洲伊人久久精品综合| 麻豆乱淫一区二区| 一级毛片黄色毛片免费观看视频| 久久久国产一区二区| 久久人人爽人人片av| 午夜老司机福利片| av女优亚洲男人天堂| 国产成人一区二区在线| 女人精品久久久久毛片| 色婷婷久久久亚洲欧美| 国产毛片在线视频| 国产成人欧美在线观看 | 80岁老熟妇乱子伦牲交| 又大又黄又爽视频免费| 国产 精品1| 老汉色av国产亚洲站长工具| 亚洲av福利一区| 黄网站色视频无遮挡免费观看| 国产成人精品福利久久| 国产激情久久老熟女| 黄色 视频免费看| 99国产综合亚洲精品| 男女边吃奶边做爰视频| 乱人伦中国视频| 亚洲欧美一区二区三区黑人| 久久久久久久精品精品| 国产免费又黄又爽又色| 久久久久精品久久久久真实原创| 久热爱精品视频在线9| 亚洲国产日韩一区二区| 一级片'在线观看视频| 激情视频va一区二区三区| 国产男女超爽视频在线观看| 国产精品一区二区在线不卡| 伊人久久大香线蕉亚洲五| 国产免费福利视频在线观看| 日日爽夜夜爽网站| 国产又色又爽无遮挡免| 日韩大片免费观看网站| 国产精品久久久久久精品古装| 欧美最新免费一区二区三区| 亚洲情色 制服丝袜| 伊人久久国产一区二区| 亚洲av国产av综合av卡| 欧美乱码精品一区二区三区| 精品少妇黑人巨大在线播放| 亚洲成国产人片在线观看| 麻豆精品久久久久久蜜桃| 亚洲成人av在线免费| 欧美少妇被猛烈插入视频| 中文字幕亚洲精品专区| 久久久亚洲精品成人影院| 国产av国产精品国产| 亚洲第一av免费看| 亚洲欧美成人综合另类久久久| 亚洲精品日本国产第一区| 99国产精品免费福利视频| 大码成人一级视频| www.自偷自拍.com| 国产精品成人在线| 老司机在亚洲福利影院| 国产av精品麻豆| 亚洲国产精品一区三区| 中文字幕制服av| 性高湖久久久久久久久免费观看| 如日韩欧美国产精品一区二区三区| 2018国产大陆天天弄谢| 精品久久久精品久久久| 男人舔女人的私密视频| 捣出白浆h1v1| 免费高清在线观看视频在线观看| 国产精品一区二区在线不卡| 久久久久国产精品人妻一区二区| 久久精品熟女亚洲av麻豆精品| 国产成人精品无人区| 精品午夜福利在线看| 色婷婷久久久亚洲欧美| 国产日韩欧美在线精品| 一级片免费观看大全| 成人漫画全彩无遮挡| 天天躁狠狠躁夜夜躁狠狠躁| 看非洲黑人一级黄片| 日韩 亚洲 欧美在线| 天天添夜夜摸| 久久青草综合色| 久久精品熟女亚洲av麻豆精品| 又大又黄又爽视频免费| 亚洲成人一二三区av| 国产免费福利视频在线观看| av国产久精品久网站免费入址| 国产伦理片在线播放av一区| 多毛熟女@视频| 亚洲国产精品国产精品| 看非洲黑人一级黄片| 久久久久国产精品人妻一区二区| 欧美日韩综合久久久久久| 在线观看国产h片| 亚洲第一av免费看| 操美女的视频在线观看| 免费黄网站久久成人精品| 午夜福利网站1000一区二区三区| 丰满少妇做爰视频| 看十八女毛片水多多多| 久久精品久久精品一区二区三区| 国产伦理片在线播放av一区| 一级毛片电影观看| 久久婷婷青草| 久久久久视频综合| 亚洲av在线观看美女高潮| 色视频在线一区二区三区| 啦啦啦中文免费视频观看日本| 国产视频首页在线观看| 久久久久精品性色| 精品国产一区二区三区四区第35| 夫妻性生交免费视频一级片| 亚洲一卡2卡3卡4卡5卡精品中文| av在线播放精品| 国产乱人偷精品视频| 国产精品.久久久| 99久久精品国产亚洲精品| 国产精品熟女久久久久浪| 久久精品aⅴ一区二区三区四区| 久久精品亚洲av国产电影网| 久久国产亚洲av麻豆专区| 乱人伦中国视频| 欧美xxⅹ黑人| 久久ye,这里只有精品| 久久久久久免费高清国产稀缺| www.av在线官网国产| 久久这里只有精品19| 我要看黄色一级片免费的| 在线观看免费视频网站a站| 十八禁人妻一区二区| 亚洲 欧美一区二区三区| 性少妇av在线| 最近手机中文字幕大全| 精品视频人人做人人爽| 两性夫妻黄色片| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 久久久亚洲精品成人影院| 观看美女的网站| 久久久久久久大尺度免费视频| 国产探花极品一区二区| 国产欧美日韩综合在线一区二区| 黄色毛片三级朝国网站| 国产日韩欧美亚洲二区| 久久精品国产亚洲av高清一级| 欧美精品av麻豆av| 飞空精品影院首页| 香蕉丝袜av| 亚洲精品国产一区二区精华液| 亚洲欧美激情在线| 精品酒店卫生间| 亚洲精品乱久久久久久| 欧美成人午夜精品| 色网站视频免费| 激情五月婷婷亚洲| 精品国产超薄肉色丝袜足j| 成人亚洲欧美一区二区av| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成人国产一区在线观看 | 亚洲情色 制服丝袜| 人人妻,人人澡人人爽秒播 | 日韩精品免费视频一区二区三区| 中文天堂在线官网| 国产精品女同一区二区软件| 在线观看www视频免费| 久久久久久久久久久免费av| 国产男人的电影天堂91| 亚洲熟女毛片儿| 久久久国产精品麻豆| 高清av免费在线| 韩国高清视频一区二区三区| videosex国产| 日本猛色少妇xxxxx猛交久久| 老汉色∧v一级毛片| 免费观看人在逋| 久久热在线av| 精品亚洲乱码少妇综合久久| 久久久久久久久久久久大奶| 国产av国产精品国产| avwww免费| 欧美日韩福利视频一区二区| 亚洲欧美中文字幕日韩二区| 90打野战视频偷拍视频| 美女扒开内裤让男人捅视频| 老熟女久久久| 操美女的视频在线观看| 免费黄频网站在线观看国产| 精品午夜福利在线看| 波多野结衣av一区二区av| 七月丁香在线播放| av在线老鸭窝| 精品久久蜜臀av无| 国产免费又黄又爽又色| www.精华液| 日韩电影二区| 久久久精品94久久精品| 亚洲精品一二三| 国产精品免费大片| 国产精品国产av在线观看| 国产男女内射视频| 国产午夜精品一二区理论片| 亚洲av电影在线观看一区二区三区| 蜜桃在线观看..| 男女国产视频网站| 亚洲国产看品久久| bbb黄色大片| 国产精品女同一区二区软件| 日本av手机在线免费观看| 精品亚洲乱码少妇综合久久| 亚洲中文av在线| 美女扒开内裤让男人捅视频| 超碰97精品在线观看| 一级毛片黄色毛片免费观看视频| 久久久久人妻精品一区果冻| 高清在线视频一区二区三区| 又大又爽又粗| 亚洲综合精品二区| 亚洲成av片中文字幕在线观看| 一本色道久久久久久精品综合| 亚洲精品在线美女| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 国产淫语在线视频| 亚洲欧美成人综合另类久久久| 水蜜桃什么品种好| 欧美日韩综合久久久久久| 9色porny在线观看| 男女无遮挡免费网站观看| 亚洲国产av新网站| 超碰97精品在线观看| 国产女主播在线喷水免费视频网站| av国产久精品久网站免费入址| 99国产综合亚洲精品| 九九爱精品视频在线观看| 成人手机av| 国产精品免费视频内射| 亚洲精品久久久久久婷婷小说| 日韩中文字幕欧美一区二区 | 操美女的视频在线观看| 久久国产精品男人的天堂亚洲| av网站在线播放免费| 久热这里只有精品99| 亚洲欧洲精品一区二区精品久久久 | 亚洲av电影在线进入| 久久ye,这里只有精品| 最近中文字幕高清免费大全6| 久久精品国产综合久久久| 久久精品熟女亚洲av麻豆精品| 别揉我奶头~嗯~啊~动态视频 | 亚洲av电影在线进入| www.熟女人妻精品国产| 亚洲av在线观看美女高潮| 久久精品亚洲av国产电影网| 久久久久网色| 90打野战视频偷拍视频| 丰满迷人的少妇在线观看| 一级爰片在线观看| 99热国产这里只有精品6| 建设人人有责人人尽责人人享有的| 黄频高清免费视频| 少妇猛男粗大的猛烈进出视频| 少妇被粗大的猛进出69影院| 操美女的视频在线观看| 如何舔出高潮| 亚洲一区二区三区欧美精品| 国产成人a∨麻豆精品| 黄片小视频在线播放| 国产精品亚洲av一区麻豆 | 国产精品免费大片| kizo精华| 看免费av毛片| 国产成人精品在线电影| 中文字幕精品免费在线观看视频| 一级毛片黄色毛片免费观看视频| 成年人午夜在线观看视频| 日韩av不卡免费在线播放| 性色av一级| 国产极品粉嫩免费观看在线| 极品少妇高潮喷水抽搐| 各种免费的搞黄视频| 中文字幕人妻丝袜一区二区 | 亚洲成人av在线免费| 国产精品熟女久久久久浪| 久久天躁狠狠躁夜夜2o2o | 捣出白浆h1v1| 大片免费播放器 马上看| 黄频高清免费视频| 日韩 亚洲 欧美在线| 99热网站在线观看| 精品亚洲成国产av| 国产亚洲午夜精品一区二区久久| 亚洲久久久国产精品| 午夜日韩欧美国产| 人成视频在线观看免费观看| 亚洲三区欧美一区| 国产精品秋霞免费鲁丝片| 免费黄网站久久成人精品| 亚洲精品乱久久久久久| 亚洲视频免费观看视频| 黄色毛片三级朝国网站| 国产免费福利视频在线观看| 国产精品一二三区在线看| 国产极品粉嫩免费观看在线| 巨乳人妻的诱惑在线观看| 激情五月婷婷亚洲| 五月开心婷婷网| 国产在视频线精品| 婷婷色麻豆天堂久久| 婷婷成人精品国产| 中国三级夫妇交换| 少妇被粗大的猛进出69影院| 日日摸夜夜添夜夜爱| 国产精品av久久久久免费| 成年人午夜在线观看视频| 久久精品亚洲av国产电影网| 免费看av在线观看网站| 香蕉国产在线看| 天堂俺去俺来也www色官网| 午夜激情av网站| 制服丝袜香蕉在线| 亚洲欧美一区二区三区久久| 黄色一级大片看看| 久久久久久久精品精品| 国产毛片在线视频| 在线观看三级黄色| 一级爰片在线观看| 视频区图区小说| 婷婷成人精品国产| 精品免费久久久久久久清纯 | 熟女av电影| 国产精品免费大片| 男女无遮挡免费网站观看| 久久这里只有精品19| 狂野欧美激情性bbbbbb| 欧美亚洲 丝袜 人妻 在线| 婷婷成人精品国产| 午夜日本视频在线| 精品少妇一区二区三区视频日本电影 | 中文字幕另类日韩欧美亚洲嫩草| 99久久99久久久精品蜜桃| 看十八女毛片水多多多| 午夜91福利影院| 你懂的网址亚洲精品在线观看| 9191精品国产免费久久| 欧美激情高清一区二区三区 | 亚洲一码二码三码区别大吗| 欧美日韩精品网址| av在线观看视频网站免费| 亚洲精品第二区| 国产精品偷伦视频观看了| 飞空精品影院首页| 日日撸夜夜添| 免费黄频网站在线观看国产| 日韩中文字幕视频在线看片| 亚洲第一av免费看| √禁漫天堂资源中文www| 日韩一区二区视频免费看| 精品午夜福利在线看| 精品国产一区二区三区四区第35| 又粗又硬又长又爽又黄的视频| 成人午夜精彩视频在线观看| 欧美少妇被猛烈插入视频| videos熟女内射| 日本色播在线视频| 18禁动态无遮挡网站| 亚洲av成人精品一二三区| 国产在线一区二区三区精| 国产一区二区在线观看av| 国产一区有黄有色的免费视频| 各种免费的搞黄视频| 亚洲第一av免费看| 夫妻午夜视频| 91精品三级在线观看| 99久久综合免费| 夫妻性生交免费视频一级片| 嫩草影院入口| 色网站视频免费| 国产在视频线精品| 国产成人精品久久二区二区91 | 午夜福利乱码中文字幕| av天堂久久9| 欧美激情 高清一区二区三区| 9色porny在线观看| 99re6热这里在线精品视频| 中文字幕av电影在线播放| 性高湖久久久久久久久免费观看| www.精华液| 午夜91福利影院| 国产一区二区三区综合在线观看| 91老司机精品| 亚洲国产日韩一区二区| 高清不卡的av网站| 精品亚洲成a人片在线观看| 亚洲av日韩精品久久久久久密 | 亚洲欧洲日产国产| 久久国产精品男人的天堂亚洲| 夫妻午夜视频| 欧美日本中文国产一区发布| 校园人妻丝袜中文字幕| 国产免费现黄频在线看| 日韩欧美精品免费久久| 纵有疾风起免费观看全集完整版| 热re99久久精品国产66热6| 亚洲欧洲国产日韩| 国产一区有黄有色的免费视频| av天堂久久9| 90打野战视频偷拍视频| 国产精品熟女久久久久浪| 美女脱内裤让男人舔精品视频| 黄色视频在线播放观看不卡| 免费观看a级毛片全部| 国产在线一区二区三区精| 在线观看三级黄色| 免费看不卡的av| 亚洲,一卡二卡三卡| 精品国产一区二区三区久久久樱花| www日本在线高清视频| 美女大奶头黄色视频| 青春草国产在线视频| 男人操女人黄网站| 国产深夜福利视频在线观看| kizo精华| 日日撸夜夜添| 青春草视频在线免费观看| 久久人人爽av亚洲精品天堂| 美女脱内裤让男人舔精品视频| 最近2019中文字幕mv第一页| 日本av免费视频播放| 免费黄网站久久成人精品| 亚洲一码二码三码区别大吗| 亚洲国产av新网站| 18禁观看日本| 99热国产这里只有精品6| 黑人巨大精品欧美一区二区蜜桃| 最近最新中文字幕免费大全7| 一个人免费看片子| 黑人猛操日本美女一级片| 国产成人一区二区在线| 热re99久久国产66热| 欧美xxⅹ黑人| 欧美成人午夜精品| 又大又黄又爽视频免费| 午夜影院在线不卡| 国产人伦9x9x在线观看| 精品免费久久久久久久清纯 | 国产精品久久久av美女十八| 啦啦啦视频在线资源免费观看| 男女边摸边吃奶| videos熟女内射| 高清不卡的av网站| 亚洲一区二区三区欧美精品| 亚洲欧美色中文字幕在线| 丁香六月欧美| 免费人妻精品一区二区三区视频| 国产日韩欧美在线精品| 亚洲婷婷狠狠爱综合网| 黄片无遮挡物在线观看|