亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        作用在二維浮體結(jié)構(gòu)上的風(fēng)荷載研究

        2020-07-25 13:51:18黃海龍杜齊魯左其華
        水運(yùn)工程 2020年7期
        關(guān)鍵詞:時(shí)距浮體恒定

        王 紅,黃海龍,2,杜齊魯,左其華,2

        (1.南京水利科學(xué)研究院,江蘇 南京 210029;2.港口航道泥沙工程交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210029;3.南京交通職業(yè)技術(shù)學(xué)院,江蘇 南京 211188)

        風(fēng)是水運(yùn)工程的重要?jiǎng)恿σ刂?,?duì)海岸工程結(jié)構(gòu)的穩(wěn)定和越浪量[1]、港口工程系泊船體的系纜力和擠靠力起著重要作用。理論上,作用在結(jié)構(gòu)上的風(fēng)或流體的壓強(qiáng)p可表示為:

        (1)

        式中:c為阻力系數(shù),與多種因素有關(guān);ρ、u分別為空氣(流體)的密度、速度。在此基礎(chǔ)上參照國(guó)際相關(guān)規(guī)范,《港口工程荷載規(guī)范》[2]給出了我國(guó)港口工程船舶風(fēng)作用荷載的計(jì)算公式:

        (2)

        式中:i為船舶受力方向,取1時(shí)為橫向x,取2時(shí)為縱向y;Ci為系數(shù),i=1時(shí)Cx=73.6×10-5,i=2時(shí)Cy=49.0×10-5;ζi、ζz為風(fēng)壓在水平方向i、垂直方向z的不均勻系數(shù);Ai為船體水面以上受風(fēng)面積。

        需要指出的是,這里的計(jì)算方法多由恒定風(fēng)作用下的理論計(jì)算與風(fēng)洞試驗(yàn)相結(jié)合得到。這樣就存在兩個(gè)問(wèn)題:1)風(fēng)洞試驗(yàn)和設(shè)計(jì)荷載計(jì)算使用的風(fēng)速是恒定的,而天然風(fēng)速是隨機(jī)的、非恒定的,即作用在船舶上的風(fēng)作用力不僅僅在空間上不均勻,在時(shí)域上也是不均勻的。隨機(jī)風(fēng)和恒定風(fēng)對(duì)船舶荷載的影響差異如何?2)雖然天然風(fēng)過(guò)程是隨機(jī)的,但限于設(shè)計(jì)條件,現(xiàn)行工程設(shè)計(jì)應(yīng)用中主要采用恒定風(fēng)。然而不同時(shí)距的平均分析將產(chǎn)生不同的設(shè)計(jì)風(fēng)速,由此得到的恒定風(fēng)速與隨機(jī)風(fēng)作用力的各特征值之間的關(guān)系如何?

        針對(duì)以上問(wèn)題,本文進(jìn)行靜水條件下的水槽試驗(yàn),測(cè)量典型二維浮體結(jié)構(gòu)不同狀況下的動(dòng)力荷載,以期通過(guò)初步分析,為今后我國(guó)港口工程風(fēng)荷載可靠度設(shè)計(jì)和研究提供參考。

        1 試驗(yàn)條件和方法

        試驗(yàn)在波浪水槽中進(jìn)行,如圖1所示。水槽長(zhǎng)40 m、寬0.8 m、深1.0 m。水槽一端配有消浪緩坡格柵。風(fēng)的模擬在加蓋水槽中進(jìn)行,造風(fēng)系統(tǒng)包括1個(gè)脈動(dòng)風(fēng)機(jī)和1個(gè)恒定風(fēng)機(jī)。脈動(dòng)風(fēng)機(jī)正對(duì)風(fēng)道口,恒定風(fēng)機(jī)鑲嵌在風(fēng)道的側(cè)壁,由風(fēng)罩將風(fēng)導(dǎo)入風(fēng)道中。試驗(yàn)風(fēng)速由兩個(gè)風(fēng)機(jī)調(diào)節(jié)完成,風(fēng)過(guò)程由恒定風(fēng)機(jī)給定1個(gè)恒定風(fēng)速,再疊加1個(gè)脈動(dòng)風(fēng)速,即實(shí)際試驗(yàn)風(fēng)速為:

        u=U+u′

        (3)

        式中:U為恒定(平均)風(fēng)速;u′為某一特定波譜產(chǎn)生的均值為零的隨機(jī)風(fēng)速。

        圖1 模型布置立面(單位:mm)

        風(fēng)速測(cè)量采用WindSonic 超聲風(fēng)速風(fēng)向儀,風(fēng)速測(cè)量范圍為0~60 m/s,精度2%,分辨率0.01 m/s,最大輸入頻率為4 Hz。模擬的隨機(jī)風(fēng)速過(guò)程見(jiàn)圖2[3]。對(duì)于1個(gè)風(fēng)速過(guò)程,采用不同時(shí)距進(jìn)行平均,可得到一系列不同的特征平均風(fēng)速。如采用1 min作為時(shí)距,則可以得到一系列風(fēng)速值,這些風(fēng)速值中的最大值,稱為1 min平均最大風(fēng)速。通常,時(shí)距越長(zhǎng)得到的平均風(fēng)速值最大值越小,時(shí)距越短得到的平均風(fēng)速最大值越大。由近海處測(cè)得資料表明:3 s陣風(fēng)風(fēng)速為1.3倍的10 min平均風(fēng)速或1.37倍的1 h平均風(fēng)速;由海岸或陸地測(cè)得資料表明:3 s陣風(fēng)風(fēng)速為1.5倍的1 h平均風(fēng)速。式(3)中的風(fēng)速U,稱為1個(gè)風(fēng)速記錄的全時(shí)距平均風(fēng)速。

        離岸荷載由拉力計(jì)測(cè)量,向岸荷載由壓力計(jì)測(cè)量,測(cè)力點(diǎn)布置如圖1所示。試驗(yàn)風(fēng)速和作用力按重力相似換算為原型。

        圖2 隨機(jī)風(fēng)速u實(shí)測(cè)過(guò)程線

        對(duì)船舶模型進(jìn)行概化處理,浮體橫截面見(jiàn)圖3。LNG船等大多數(shù)船舶均可以此斷面概化,其幾何比尺約為1:80,相應(yīng)的時(shí)間比尺約為1:9??紤]3種不同深度的吃水,即船體高出靜水面的高度分別為11.0、13.0及14.6 cm,相當(dāng)于原型的8.8、10.4及11.7 m,也分別對(duì)應(yīng)系泊浮體滿載、半載、壓載狀態(tài)。

        圖3 試驗(yàn)浮體橫截面(單位:m)

        2 結(jié)果分析

        2.1 恒定風(fēng)的作用

        浮體在3種不同吃水(滿載、半載及壓載)條件下的恒定離岸風(fēng)對(duì)浮體的作用力見(jiàn)圖4。針對(duì)3種不同裝載狀態(tài)的浮體,將實(shí)測(cè)試驗(yàn)結(jié)果與式(1)及規(guī)范公式計(jì)算出的荷載進(jìn)行比較,結(jié)果顯示,現(xiàn)行規(guī)范結(jié)果都比理論公式結(jié)果略大,即有一定的安全余度,但都比實(shí)測(cè)值小。在滿載狀態(tài)下,風(fēng)對(duì)系泊浮體的作用比規(guī)范計(jì)算值大10%左右;在半載狀態(tài)下,比規(guī)范計(jì)算值大約9%;在壓載狀態(tài)時(shí),比規(guī)范計(jì)算值大約11%。

        圖4 恒定風(fēng)作用力與風(fēng)速變化關(guān)系

        這種系統(tǒng)性偏差應(yīng)與其作用機(jī)制有關(guān),并主要由于風(fēng)作用的靜水下邊界與風(fēng)道試驗(yàn)的差異產(chǎn)生。在靜水條件下,風(fēng)對(duì)浮體作用荷載可表示為:

        F=F1+F2

        (4)

        式中:F1為風(fēng)直接作用荷載,可以理解為傳統(tǒng)風(fēng)道試驗(yàn)或理論計(jì)算部分;F2為風(fēng)間接作用荷載,主要由風(fēng)引起的流體表面風(fēng)吹流和局部風(fēng)成浪兩部分組成。這兩部分是理論計(jì)算和風(fēng)道試驗(yàn)中無(wú)法表現(xiàn)的,由于水流密度比空氣密度大得多,由式(1)可知,其量級(jí)雖小,但會(huì)產(chǎn)生一定的荷載。同時(shí),隨著港口建設(shè)的發(fā)展,港池水域面積越來(lái)越大,這兩部分的影響可能也會(huì)越來(lái)越大,應(yīng)在風(fēng)荷載設(shè)計(jì)中引起足夠的重視。

        2.2 隨機(jī)風(fēng)的作用力

        結(jié)合文獻(xiàn)[4],本文選用API譜[5]和Davenport譜[6]兩種隨機(jī)風(fēng)譜。前者常用于海向來(lái)風(fēng),后者常用于岸向來(lái)風(fēng)。

        2.2.1API譜隨機(jī)風(fēng)作用力

        API譜隨機(jī)風(fēng)作用力不同特征值與風(fēng)速的關(guān)系如圖5所示(圖中恒定風(fēng)風(fēng)速指全時(shí)段平均風(fēng)速)。API隨機(jī)風(fēng)作用下,不同載態(tài)作用力相對(duì)恒定風(fēng)作用均有提高,以1/3系纜力特征值來(lái)看,壓載時(shí)提高60%~140%,半載時(shí)提高40%~80%,滿載時(shí)提高40%~94%;以1/10系纜力特征值來(lái)看,壓載時(shí)提高80%~190%,半載時(shí)提高60%~80%,滿載時(shí)提高67%~120%。滿載時(shí)浮體所受風(fēng)最大作用力是恒定風(fēng)的1.9~2.1倍,半載時(shí)為1.7~1.8倍,壓載時(shí)約為2倍。

        圖5 API譜隨機(jī)風(fēng)3種狀態(tài)的作用力特征值

        2.2.2Davenport譜隨機(jī)風(fēng)作用力

        壓載、半載、滿載狀態(tài),Davenport譜隨機(jī)風(fēng)作用下,作用力各特征值與風(fēng)速之間的關(guān)系見(jiàn)圖6。同等風(fēng)速下,Davenport隨機(jī)風(fēng)作用下的系纜力與恒定風(fēng)下的系纜力關(guān)系與API隨機(jī)作用時(shí)變化相似。

        Davenport隨機(jī)風(fēng)相對(duì)恒定風(fēng)作用均有提高,以1/3系纜力特征值來(lái)看,壓載時(shí)提高20%~50%,半載和滿載時(shí)均提高10%~70%;以1/10系纜力特征值來(lái)看,壓載時(shí)提高30%~70%,半載時(shí)提高20%~70%,滿載時(shí)提高20%~80%。

        圖6 Davenport譜隨機(jī)風(fēng)3種狀態(tài)的作用力特征值

        Davenport譜特性隨機(jī)風(fēng)場(chǎng)中,滿載時(shí)船體所受風(fēng)最大作用力是恒定風(fēng)的1.7倍,半載時(shí)為1.3~1.5倍,壓載時(shí)為1.6~1.7倍。

        由圖6可以看出,當(dāng)采用全時(shí)段平均風(fēng)速作為恒定風(fēng)速作用在浮體上的作用力是相同的,但明顯小于API譜和Davenport譜隨機(jī)風(fēng)作用所得結(jié)果,這表明采用全時(shí)段平均風(fēng)速作為浮式結(jié)構(gòu)的設(shè)計(jì)風(fēng)速是不合適的。

        2.2.3不同風(fēng)速譜作用荷載比較

        壓載狀態(tài)時(shí),API譜與Davenport譜隨機(jī)風(fēng)作用荷載特征值比值R見(jiàn)表1。

        表1 API譜和Davenport譜隨機(jī)風(fēng)作用荷載特征值比值

        由表1可知,不同的風(fēng)速譜所得的作用力特征值是不同的,且API譜所得結(jié)果整體上大于Davenport譜所得結(jié)果,其差異隨特征值的加大有增大的趨勢(shì)。這種差異主要是由于與Davenport譜相比,API譜中包含更多低頻部分的風(fēng)能量。對(duì)于浮體結(jié)構(gòu)來(lái)講,低頻荷載會(huì)比高頻的影響更大些。由此,可以建議對(duì)于船舶結(jié)構(gòu)風(fēng)的影響宜采用API譜進(jìn)行估算。

        2.3 不同時(shí)距風(fēng)速特征荷載分析

        工程設(shè)計(jì)中究竟用何種時(shí)距平均(如3 s、1、3、5、10 、15 min和1 h等的特征風(fēng)速值)來(lái)進(jìn)行設(shè)計(jì)尚未形成統(tǒng)一的意見(jiàn)。不同國(guó)家的規(guī)范有不同設(shè)計(jì)風(fēng)速確定方法,如美國(guó)土木工程師協(xié)會(huì)ASCE、澳大利亞新風(fēng)險(xiǎn)管理AS/NZ、國(guó)際標(biāo)準(zhǔn)化組織 ISO等標(biāo)準(zhǔn)對(duì)房屋結(jié)構(gòu)基本風(fēng)速取 3 s的陣風(fēng)速度;日本建筑學(xué)會(huì)AIJ、 歐洲EU 標(biāo)準(zhǔn)等取10 min 平均值;加拿大建筑規(guī)范NBCC取1 h平均值。在海上,日本采用的是瞬時(shí)最大風(fēng)速,英國(guó)和澳大利亞采用3 s風(fēng)速,加拿大采用1 h平均風(fēng)速,美國(guó)的時(shí)距不是一個(gè)定值,其值在0.5~1.0 min之間。用全時(shí)段平均風(fēng)速作為設(shè)計(jì)風(fēng)速偏于危險(xiǎn)。

        根據(jù)我國(guó)的《港口工程荷載規(guī)范》,作用在港口工程結(jié)構(gòu)上的風(fēng)荷載風(fēng)速取“重現(xiàn)期50 a一遇10 min平均最大風(fēng)速”,在其附錄E中規(guī)定作用在船舶上的風(fēng)荷載計(jì)算,“對(duì)沒(méi)有特殊要求的港口,可按9級(jí)風(fēng)即風(fēng)速20.8~24.4 m/s考慮”。這是考慮現(xiàn)有風(fēng)速資料以及相關(guān)規(guī)范常以10 min作為標(biāo)準(zhǔn)時(shí)距,而其他時(shí)距的風(fēng)速資料難以獲得等因素,且研究結(jié)果不足以確定合適的值而采取經(jīng)驗(yàn)的方法。這對(duì)于風(fēng)作用力的可靠度設(shè)計(jì)帶來(lái)很多不確定性。

        由于風(fēng)的隨機(jī)性,也帶來(lái)了作用力的不確定性,可以由不同的特征值,如平均值、1/3值、1/10值和最大值來(lái)表示。

        2.3.1API譜不同時(shí)距的風(fēng)速作用力

        滿載狀態(tài)下,API譜特性隨機(jī)風(fēng)不同時(shí)距下,作用力特征值與不同時(shí)距風(fēng)速的變化關(guān)系見(jiàn)圖7。以1/3特征值為基準(zhǔn)時(shí),時(shí)距采用1 min平均風(fēng)速,恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以1/10特征值為基準(zhǔn)時(shí),時(shí)距采用10 ~30 s平均風(fēng)速,恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以最大峰值為基準(zhǔn)時(shí),當(dāng)時(shí)距采用10 s時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng)。同時(shí)可以看到,時(shí)距采用10~60 min時(shí),作用力的各特征值差別不大,主要是此時(shí)距的風(fēng)速趨于穩(wěn)定值。

        圖7 滿載狀態(tài)API譜隨機(jī)風(fēng)作用力特征值與不同時(shí)距的風(fēng)速關(guān)系

        壓載狀態(tài)下,API譜特性隨機(jī)風(fēng)不同時(shí)距下,作用力特征值與風(fēng)速變化關(guān)系見(jiàn)圖8。以1/3特征值為基準(zhǔn)時(shí),時(shí)距采用30~60 s時(shí),相同風(fēng)速下,恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以1/10特征值為基準(zhǔn)時(shí),時(shí)距采用10~30 s時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);對(duì)于作用力最大峰值,時(shí)距只有采用10 s時(shí),恒定風(fēng)的作用效果才與隨機(jī)風(fēng)的作用效果相當(dāng)。

        圖8 壓載狀態(tài)API譜隨機(jī)風(fēng)作用力特征值與不同時(shí)距的風(fēng)速關(guān)系

        試驗(yàn)結(jié)果還表明,半載狀態(tài)下,以1/3特征值為基準(zhǔn),時(shí)距采用1 min時(shí),同等風(fēng)速下,恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以1/10特征值為基準(zhǔn),時(shí)距采用10~30 s時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);對(duì)于作用力最大峰值,當(dāng)時(shí)距采用10 s時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng)。

        2.3.2Davenport譜不同時(shí)距的作用力分析

        滿載狀態(tài)下,Davenport譜特性隨機(jī)風(fēng)不同時(shí)距下,作用力特征值與不同時(shí)距平均風(fēng)速變化關(guān)系見(jiàn)圖9。以1/3特征值為基準(zhǔn),時(shí)距采用2 min時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以1/10特征值為基準(zhǔn),時(shí)距采用30~60 s時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);對(duì)于作用力最大峰值,時(shí)距采用10 s時(shí),恒定風(fēng)的作用效果才與隨機(jī)風(fēng)的作用效果相當(dāng)。

        圖9 滿載狀態(tài)Davenport譜隨機(jī)風(fēng)作用力特征值與不同時(shí)距風(fēng)速關(guān)系

        壓載狀態(tài)下,Davenport譜特性隨機(jī)風(fēng)不同時(shí)距下,作用力特征值與不同時(shí)距平均風(fēng)速變化關(guān)系為:以1/3特征值為基準(zhǔn),時(shí)距采用2 min時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以1/10特征值為基準(zhǔn),時(shí)距采用1 min時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以作用力最大峰值為基準(zhǔn),時(shí)距采用10 s時(shí),兩者作用相當(dāng)。

        試驗(yàn)結(jié)果還表明,在半載狀態(tài)下,以1/3特征值為基準(zhǔn),時(shí)距采用1~2 min時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以1/10特征值為基準(zhǔn)時(shí),當(dāng)時(shí)距采用30 s時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng);以最大峰值為基準(zhǔn),當(dāng)時(shí)距采用10 s時(shí),恒定風(fēng)的作用效果與隨機(jī)風(fēng)的作用效果相當(dāng)。

        針對(duì)本文試驗(yàn)狀況,Davenport譜特性隨機(jī)風(fēng)比API譜特性隨機(jī)風(fēng)對(duì)時(shí)距選擇的波動(dòng)性小。對(duì)3種載態(tài)下的作用力特征值分析,所得的時(shí)距選取值與特征值關(guān)系是一致的,比如,研究1/3作用力特征值宜采用1~2 min時(shí)距風(fēng)速;研究1/10作用力特征值宜采用30~60 s時(shí)距風(fēng)速;研究最大作用力宜采用10 s時(shí)距風(fēng)速。

        由不同載態(tài)浮體作用力特征值對(duì)比分析看,采用某一固定時(shí)距風(fēng)速平均最大值作為恒定風(fēng)速來(lái)模擬浮體實(shí)際受到的風(fēng)力可能是不妥當(dāng)?shù)?,時(shí)距越大,得到的作用力越小。如恒定風(fēng)采用10 min平均風(fēng)速,其作用力往往偏小,對(duì)工程設(shè)計(jì)來(lái)說(shuō)可能偏于危險(xiǎn)。使用恒定風(fēng)作為結(jié)構(gòu)作用的設(shè)計(jì)值,除應(yīng)根據(jù)設(shè)計(jì)的具體情況所需的作用力特征值選取合適的時(shí)距平均值外,還應(yīng)考慮浮體的載態(tài)、譜特性、隨機(jī)風(fēng)的類別等因素。

        3 結(jié)論

        1) 現(xiàn)有的風(fēng)作用力規(guī)范計(jì)算所得荷載與理論值相比尚有一定的裕度,但忽略了由于水邊界存在而產(chǎn)生的間接風(fēng)作用力的影響,恒定風(fēng)作用力仍偏小10%左右。

        2) 恒定風(fēng)不考慮風(fēng)作用的脈動(dòng)性,從而低估了作用力對(duì)工程的影響。由于天然風(fēng)是隨機(jī)的,因此在工程實(shí)踐中,應(yīng)考慮隨機(jī)風(fēng)脈動(dòng)作用產(chǎn)生的荷載。

        3)將API譜隨機(jī)風(fēng)和Davenport譜隨機(jī)風(fēng)相比,前者引起的系泊浮體運(yùn)動(dòng)大于后者,并且二者相差可能高達(dá)30%~70%。建議海岸工程浮體設(shè)計(jì)和研究中隨機(jī)風(fēng)譜取包含較多低頻能量的API類譜。

        4)從對(duì)不同裝載狀態(tài)時(shí)3種作用力特征值的分析比較看,用恒定風(fēng)來(lái)模擬隨機(jī)風(fēng)對(duì)浮體的作用,不能簡(jiǎn)單地認(rèn)為或采用某一個(gè)固定的風(fēng)速時(shí)距,風(fēng)速時(shí)距應(yīng)該是一個(gè)變量。在研究風(fēng)對(duì)浮體的作用時(shí),應(yīng)根據(jù)不同設(shè)計(jì)所需求的特征作用力,采用不同的風(fēng)速時(shí)距。

        猜你喜歡
        時(shí)距浮體恒定
        浮體結(jié)構(gòu)沉浮過(guò)程周圍水流特性研究
        趣味性進(jìn)度提示設(shè)計(jì)對(duì)時(shí)距知覺(jué)的影響
        物探船硬浮體陣列自擴(kuò)變量分析與應(yīng)用
        超大型浮體結(jié)構(gòu)碰撞損傷研究
        花花世界
        學(xué)習(xí)障礙學(xué)生與普通學(xué)生短時(shí)距估計(jì)的比較
        有限流動(dòng)水域浮體受力及側(cè)傾研究
        回溯式時(shí)距估計(jì)的計(jì)時(shí)機(jī)制*
        漫畫(huà)十萬(wàn)個(gè)為什么
        標(biāo)量計(jì)時(shí)模型的影響因素及發(fā)展*
        精品国产乱子伦一区二区三| 另类亚洲欧美精品久久不卡| 亚洲国产成a人v在线观看| 久久亚洲国产高清av一级 | 免费a级毛片无码a∨免费软件| 免费无码中文字幕A级毛片| 久久国产精品懂色av| 国产一区二区视频免费在| 午夜精品久久久久久99热| 久久成年片色大黄全免费网站| 亚洲一区二区三区在线观看| 亚洲成人av在线第一页| 在线观看在线观看一区二区三区| 一本久久综合亚洲鲁鲁五月天| 中文字幕av一区中文字幕天堂| 亚洲福利视频一区| 亚洲一道一本快点视频| 亚洲一二三四区免费视频| 亚洲欧美精品suv| 在线免费观看韩国a视频| 手机在线观看亚洲av| 四季极品偷拍一区二区三区视频| 18禁裸男晨勃露j毛网站| 五月天精品视频在线观看| 中文字幕亚洲精品码专区| 久久久麻豆精亚洲av麻花| 亚洲日韩一区精品射精| 日韩精品区欧美在线一区| 亚洲一区二区三区视频免费 | 五十路熟女一区二区三区| 青青草针对华人超碰在线| 穿着白丝啪啪的av网站| 亚洲成a v人片在线观看| 国产在线不卡免费播放| 国家一级内射高清视频| 最近免费中文字幕中文高清6| 一二三四在线观看免费视频| 国产精品黄色片在线观看| 久久精品网站免费观看| 久久777国产线看观看精品| a级毛片在线观看|