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

    浸沒式外輸浮筒的動(dòng)力響應(yīng)影響

    2022-11-01 09:54:50張佳敏竇宏波
    中國海洋平臺(tái) 2022年5期
    關(guān)鍵詞:浮筒系泊阻尼比

    李 英, 張佳敏, 王 坤, 竇宏波

    (1. 天津大學(xué) 水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室, 天津 300072;2. 中海油研究總院有限責(zé)任公司, 北京 100028)

    0 引 言

    單點(diǎn)系泊系統(tǒng)是一種廣泛用于海上平臺(tái)定位作業(yè)的系泊方式,其通常由穿梭油船、系船纜、漂浮軟管、外輸浮筒和系泊線組成,是原油傳輸?shù)闹修D(zhuǎn)站。浮筒在環(huán)境載荷的作用下,與系泊線、穿梭油船之間耦合運(yùn)動(dòng),其運(yùn)動(dòng)響應(yīng)對(duì)安全輸送油氣至關(guān)重要。

    目前,國內(nèi)外許多學(xué)者對(duì)浮筒的運(yùn)動(dòng)響應(yīng)進(jìn)行分析和研究。張磊等[1]對(duì)深水浮筒進(jìn)行水動(dòng)力分析并與試驗(yàn)結(jié)果進(jìn)行對(duì)比,發(fā)現(xiàn)西非海域的風(fēng)浪和涌浪對(duì)浮筒的作用差別很大。SALEM等[2]在頻域內(nèi)將二次阻力/阻尼進(jìn)行線性化處理,將時(shí)域、頻域結(jié)果與試驗(yàn)?zāi)P蜏y試結(jié)果進(jìn)行對(duì)比驗(yàn)證,結(jié)果證明所用的線性化方法對(duì)浮筒縱搖響應(yīng)峰值有很好的估計(jì)效果。文曉東[3]分別對(duì)極限工況和外輸作業(yè)工況下浮筒的運(yùn)動(dòng)響應(yīng)和系泊錨索的張力響應(yīng)進(jìn)行分析預(yù)測,HWANG[4]也進(jìn)行類似的研究。姜季江[5]分析六錨鏈懸鏈?zhǔn)絾吸c(diǎn)系泊浮筒在不同錨鏈和錨點(diǎn)分布方式中的系泊力變化。

    為減小浮筒的垂蕩和縱搖運(yùn)動(dòng)響應(yīng),工程上會(huì)在浮筒上安裝裙板以增大浮筒的垂蕩和縱搖阻尼。RYU等[6]對(duì)深水浮筒垂蕩板的縱蕩、垂蕩和縱搖運(yùn)動(dòng)進(jìn)行預(yù)測,發(fā)現(xiàn)縱搖運(yùn)動(dòng)對(duì)裙板的阻力效應(yīng)特別敏感??登f等[7]提出一種莫里森單元與蝶形單元相組合的方法以計(jì)算浮筒的黏性載荷??涤袨閇8]分別分析帶裙板式浮筒主體尺度的變化、裙板位置和外徑的變化對(duì)水面處外輸浮筒運(yùn)動(dòng)特性的影響,發(fā)現(xiàn)增大浮筒直徑可改善浮筒的垂蕩運(yùn)動(dòng),而對(duì)浮筒的縱蕩運(yùn)動(dòng)幾乎沒有影響。孫世鵬[9]也對(duì)裙板半徑、板厚等因素進(jìn)行類似的研究。另外,浮筒的裙板與Spar平臺(tái)的垂蕩板結(jié)構(gòu)和作用相似,滕斌等[10]對(duì)Spar平臺(tái)垂蕩板的水動(dòng)力系數(shù)進(jìn)行計(jì)算,并與試驗(yàn)結(jié)果進(jìn)行對(duì)比,發(fā)現(xiàn)在真實(shí)流體中垂蕩板的輻射阻尼在總阻尼中所占的百分比較小。

    上述研究都針對(duì)漂浮在水面處的浮筒,目前關(guān)于完全浸沒在水下的外輸浮筒研究較少。不同于水面處的浮筒,浸沒式浮筒由于其浮力大于重力,需要靠系泊線的張力補(bǔ)償才能在一定位置達(dá)到平衡,其水動(dòng)力分析方法也不同于水面處的浮筒。另外,浮筒主尺度變化對(duì)其動(dòng)力響應(yīng)的影響趨勢也與水面處浮筒有所差別。因此,本文基于邊際油田開發(fā)實(shí)際項(xiàng)目中遇到的一個(gè)概念性設(shè)計(jì)的浸沒式外輸浮筒,研究其附加阻尼、幾何形狀對(duì)浮筒水動(dòng)力特性和運(yùn)動(dòng)響應(yīng),以及系泊線張力的影響,進(jìn)而優(yōu)化浸沒式外輸浮筒設(shè)計(jì),為實(shí)際工程應(yīng)用提供指導(dǎo)。

    1 浮筒基本參數(shù)

    所研究的浸沒式浮筒的作業(yè)水深為40 m,浮筒主尺度參數(shù)和相關(guān)信息如表1所示。

    表1 浮筒主尺度信息

    由于浸沒式浮筒浮力大于重力,因此采用張緊式系泊系統(tǒng)對(duì)其進(jìn)行系泊定位,系泊線的布置采用3×2的分組形式,每組系泊線夾角為120°,組內(nèi)相鄰2根系泊線的夾角為4°,布置如圖1所示。每根系泊線由錨鏈和鋼纜兩種材料組成,其主要參數(shù)如表2所示。

    注:圖中編號(hào)ML1~ML6分別代表錨鏈1~錨鏈6圖1 系泊系統(tǒng)布置示例

    表2 系泊線組成

    2 浮筒水動(dòng)力模型

    應(yīng)用ANSYS APDL建立浮筒濕表面模型,浮筒的有限元模型如圖2所示。

    圖2 浮筒有限元模型

    模型網(wǎng)格的劃分對(duì)浮筒的水動(dòng)力響應(yīng)有一定的影響,因此在進(jìn)行水動(dòng)力分析之前,先對(duì)網(wǎng)格劃分的質(zhì)量進(jìn)行評(píng)估。AQWA默認(rèn)使用遠(yuǎn)場法進(jìn)行波浪平均漂移力的計(jì)算,網(wǎng)格劃分質(zhì)量不影響遠(yuǎn)場法計(jì)算結(jié)果;近場法通過對(duì)濕表面進(jìn)行積分來求解平均漂移力,計(jì)算結(jié)果受網(wǎng)格劃分質(zhì)量的影響。圖3給出0°波浪入射角度下分別用遠(yuǎn)場法和近場法計(jì)算的縱蕩方向二階漂移力,二者幾乎沒有差別,表明該模型的網(wǎng)格劃分合理。由于遠(yuǎn)場法只能計(jì)算3個(gè)自由度的平均漂移力,因此最終采用近場法。

    圖3 0°方向縱蕩二階漂移力

    3 黏性阻尼對(duì)浮筒運(yùn)動(dòng)特性的影響

    3.1 水動(dòng)力分析結(jié)果

    AQWA是一款基于勢流理論的計(jì)算軟件,只能考慮輻射阻尼,無法考慮黏性的影響,因此應(yīng)用AQWA-Line進(jìn)行水動(dòng)力分析需要額外考慮黏性阻尼。為研究黏性阻尼對(duì)浸沒式外輸浮筒水動(dòng)力性能的影響,根據(jù)式(1)計(jì)算臨界阻尼Dc,分別考慮5%、8%和10%阻尼比的臨界阻尼作為黏性阻尼。

    (1)

    式中:M為浮筒對(duì)應(yīng)自由度的質(zhì)量,包括慣性質(zhì)量和附加質(zhì)量;K為對(duì)應(yīng)自由度的剛度。

    由于所研究的浮筒完全浸沒在水下,浮力大于重力,需要靠系泊線保持穩(wěn)定,因此為準(zhǔn)確預(yù)報(bào)浮筒真實(shí)的水動(dòng)力特性,需要充分考慮系泊的影響。由于浮筒的對(duì)稱性,僅給出0°~90°波浪入射方向的浮筒水動(dòng)力計(jì)算結(jié)果。阻尼比為10%時(shí)浸沒式浮筒的縱蕩、垂蕩和縱搖位移響應(yīng)幅值算子(Response Amplitude Operator,RAO)曲線如圖4所示。由圖4可知:由于浮筒的對(duì)稱性,浮筒的垂蕩運(yùn)動(dòng)與波浪入射角度無關(guān);由于是張緊式系泊系統(tǒng),浮筒在3個(gè)自由度方向的固有周期都比較小,縱蕩、垂蕩、縱搖方向固有周期分別約5.0 s、5.3 s 和3.5 s,遠(yuǎn)離波浪的主要能量范圍。該浸沒式浮筒的固有周期特性與張力腿平臺(tái)(Tension Leg Platform,TLP)有相似之處。

    圖4 10%阻尼比浮筒位移RAO

    在0°波浪入射方向下,不同阻尼比對(duì)該浸沒式浮筒縱蕩、垂蕩和縱搖方向位移RAO的影響如圖5所示,相應(yīng)的幅值如圖6所示。由圖6可知:隨著阻尼比增大,浮筒縱搖方向RAO響應(yīng)幅值小幅度降低,而縱蕩、垂蕩方向RAO響應(yīng)幅值變化不大。由于浸沒式浮筒采用張緊式系泊,在進(jìn)行水動(dòng)力分析時(shí)考慮了系泊的作用,浮筒本身的運(yùn)動(dòng)響應(yīng)就比較小,采取5%、8%和10%的阻尼比整體來說對(duì)浮筒的運(yùn)動(dòng)響應(yīng)影響不大。

    圖5 不同阻尼比下0°方向浮筒位移RAO

    圖6 不同阻尼比下位移RAO響應(yīng)幅值

    3.2 時(shí)域動(dòng)態(tài)計(jì)算結(jié)果

    3.2.1 浮筒位移結(jié)果

    為研究阻尼比對(duì)浮筒時(shí)域運(yùn)動(dòng)響應(yīng)和系泊線張力的影響,開展浮筒-系泊系統(tǒng)的動(dòng)態(tài)時(shí)域分析,針對(duì)一年一遇環(huán)境條件,選用JONSWAP隨機(jī)波浪譜,模擬海況時(shí)長為3 h,分別得到該浮筒在縱蕩、垂蕩和縱搖等3個(gè)自由度方向的位移時(shí)程曲線,各自由度的時(shí)域位移最大值均出現(xiàn)在6 800~7 000 s。

    當(dāng)波浪入射角度為0°、模擬時(shí)間為6 800~7 000 s 時(shí),浮筒3個(gè)自由度方向的位移時(shí)程曲線如圖7所示。由圖7可知,隨著阻尼比的增大,浮筒在縱蕩、垂蕩和縱搖等3個(gè)自由度方向的時(shí)域位移峰值均減小。當(dāng)波浪入射角為0°、阻尼比分別為5%、8%和10%時(shí),相應(yīng)自由度時(shí)域位移峰值的變化情況如圖8所示。

    圖7 浮筒位移時(shí)程曲線

    圖8 不同阻尼比下浮筒時(shí)域位移峰值

    當(dāng)阻尼比由5%變化至8%時(shí),垂蕩運(yùn)動(dòng)峰值減小30%,縱搖運(yùn)動(dòng)峰值減小18%,垂蕩和縱搖位移運(yùn)動(dòng)響應(yīng)分別成約10倍、6倍變化。阻尼比對(duì)該浸沒式浮筒垂蕩和縱搖方向上時(shí)域位移運(yùn)動(dòng)的影響都比較大,說明該浸沒式浮筒的垂蕩、縱搖運(yùn)動(dòng)對(duì)阻尼比較敏感,可以考慮安裝裙板等結(jié)構(gòu)以增大浮筒的垂蕩和縱搖阻尼,進(jìn)而減小浮筒的垂蕩和縱搖運(yùn)動(dòng)響應(yīng)。

    3.2.2 系泊線張力結(jié)果

    當(dāng)波浪入射角度為0°時(shí),6根系泊纜中ML4受到的張力最大。不同阻尼比下ML4的最大有效張力結(jié)果如圖9所示。 由圖9可知,隨著阻尼比的增大,ML4的最大有效張力減小。與阻尼比為5%的結(jié)果比較,在10%阻尼比時(shí)最大系泊張力下降28%,系統(tǒng)更安全。

    圖9 不同阻尼比下ML4最大有效張力

    4 幾何形狀對(duì)浮筒運(yùn)動(dòng)性能的影響

    為優(yōu)化浮筒設(shè)計(jì),研究浸沒式浮筒的幾何形狀對(duì)浮筒水動(dòng)力性能的影響,在直徑10 m、高3 m的浮筒的基礎(chǔ)上,保持浮筒的整體排水量、重心和慣性半徑不變,建立不同幾何形狀的浮筒有限元模型進(jìn)行水動(dòng)力計(jì)算,在水動(dòng)力分析中統(tǒng)一考慮10%阻尼比。不同幾何形狀浮筒的主尺度如表3所示。

    表3 浮筒主尺度

    4.1 水動(dòng)力分析結(jié)果

    4.1.1 附加質(zhì)量和阻尼

    對(duì)不同幾何形狀的浮筒分別進(jìn)行水動(dòng)力分析,得到帶系泊的浸沒式浮筒在縱蕩、垂蕩和縱搖等3個(gè)自由度方向的附加質(zhì)量、輻射阻尼結(jié)果分別如圖10和圖11所示。

    圖10 浮筒附加質(zhì)量

    圖11 浮筒輻射阻尼

    由圖10和圖11可知:附加質(zhì)量和輻射阻尼的結(jié)果變化趨勢相同,隨著浸沒式浮筒直徑的增大,縱蕩方向的附加質(zhì)量、輻射阻尼都逐漸減小,垂蕩和縱搖方向的附加質(zhì)量、輻射阻尼都逐漸增大。直徑為12 m 的浮筒附加質(zhì)量在縱蕩方向上比直徑為8 m 的浮筒減小約1/2,在垂蕩和縱搖方向上分別增大約3倍和7倍。圖11表明直徑為12 m 的浮筒與直徑為8 m 的浮筒相比,輻射阻尼在縱蕩方向上減小約1/2,在垂蕩和縱搖方向上分別增大約2倍和6倍。其中,浮筒在縱搖方向上的水動(dòng)力特性對(duì)幾何形狀更敏感。

    4.1.2 位移RAO對(duì)比分析

    在0°波浪入射方向下不同幾何形狀對(duì)浸沒式浮筒的縱蕩、垂蕩和縱搖等3個(gè)自由度方向的位移RAO的影響如圖12所示。由圖12可知,隨著浸沒式浮筒直徑的增大,縱蕩和縱搖RAO響應(yīng)幅值均減小,但是垂蕩方向上的RAO響應(yīng)幅值增大。這是由于隨著浮筒直徑的增大,圓柱形浮筒在縱蕩方向上的投影面積逐漸減小,因此受到的波浪作用力也逐漸減小,但是浮筒在底面上的投影面積逐漸增大,因此在垂蕩方向上受到的波浪作用力也逐漸增大。另外,隨著浮筒直徑的增加,其在縱搖方向上受到的波浪力矩雖然增加,但是由于浮筒在縱搖方向上的附加質(zhì)量和阻尼增加得更多,因此,浮筒在縱搖方向上的響應(yīng)幅值反而變小。由圖12(c)可知,直徑為8 m和9 m的浮筒在縱搖方向的RAO響應(yīng)幅值遠(yuǎn)大于直徑為10 m、11 m和12 m的浮筒,直徑8 m浮筒的縱搖RAO幅值約為直徑10 m浮筒的4倍。因此,基于頻域運(yùn)動(dòng)RAO設(shè)計(jì)浮筒,直徑8 m和9 m的浮筒幾何形狀顯然不合適。

    圖12 不同幾何形狀下浮筒位移RAO

    不同幾何形狀下浮筒相應(yīng)自由度固有周期的變化如表4所示,隨著浸沒浮筒直徑的增大,浮筒幾何形狀變得矮胖,浮筒自身的縱蕩、垂蕩、縱搖固有周期逐漸變大,越來越接近波浪的主要能量范圍。

    表4 不同幾何形狀浮筒固有周期 s

    4.2 時(shí)域動(dòng)態(tài)計(jì)算位移結(jié)果

    4.2.1 浮筒位移結(jié)果

    為研究浮筒幾何形狀對(duì)浮筒時(shí)域運(yùn)動(dòng)響應(yīng)和系泊線張力的影響,針對(duì)一年一遇的環(huán)境條件,應(yīng)用JONSWAP建立系泊系統(tǒng)模型,開展浮筒-系泊系統(tǒng)時(shí)域內(nèi)的動(dòng)態(tài)耦合分析。模擬時(shí)長為3 h,在時(shí)域模擬內(nèi)各自由度位移的最大值均出現(xiàn)在6 800~7 000 s。當(dāng)波浪入射角度為0°、模擬時(shí)間為6 800~7 000 s時(shí),浮筒縱蕩、垂蕩和縱搖等3個(gè)自由度方向的位移時(shí)程曲線如圖13所示。

    圖13 浮筒位移時(shí)程曲線

    由圖13可知,隨著浮筒幾何直徑的增大,浮筒在縱蕩、縱搖方向上的位移均增大,不同于浸沒式浮筒水動(dòng)力分析的位移RAO結(jié)果,縱蕩方向和縱搖方向的RAO響應(yīng)幅值均減小。這可能是由于浮筒的固有周期、運(yùn)動(dòng)的耦合效應(yīng)、在頻域水動(dòng)力分析中系泊剛度是固定的而在時(shí)域分析下系泊剛度是變化的以及系泊線本身的阻尼影響等。垂蕩方向RAO響應(yīng)幅值無論是位移RAO還是時(shí)域動(dòng)態(tài)分析結(jié)果都隨著直徑增大而增大。因此,對(duì)于該浸沒式浮筒,頻域內(nèi)的水動(dòng)力分析并不能完全準(zhǔn)確預(yù)報(bào)其水動(dòng)力響應(yīng),非常有必要開展時(shí)域內(nèi)的動(dòng)態(tài)分析。

    4.2.2 系泊線張力結(jié)果

    系泊系統(tǒng)耦合的時(shí)域動(dòng)態(tài)分析表明,當(dāng)波浪入射角度為0°時(shí)ML4上的系泊張力最大。不同幾何形狀下浮筒系泊系統(tǒng)ML4最大有效張力變化情況如圖14所示。由圖14可知,系泊線時(shí)域最大有效張力的變化與時(shí)域分析的浮筒位移變化趨勢一致,隨著浮筒幾何直徑的增大,ML4的最大有效張力增加。該浸沒式圓柱形浮筒的直徑越小,形狀越細(xì)長,浮筒的位移運(yùn)動(dòng)響應(yīng)和系泊線的最大有效張力越小。

    圖14 不同幾何形狀下ML4最大有效張力

    與直徑為8 m的浮筒ML4最大有效張力相比,不同直徑浮筒的ML4最大有效張力的無量綱比值如表5所示。同時(shí)結(jié)合表4,浮筒直徑為8 m時(shí)固有周期最小,更遠(yuǎn)離波浪的主要能量范圍,因此,時(shí)域模擬結(jié)果表明,直徑8 m的浮筒是這5種不同幾何形狀的最優(yōu)設(shè)計(jì),此結(jié)論與基于頻域分析的RAO的結(jié)論正好相反,這表明此類型的浮筒優(yōu)化設(shè)計(jì)要基于時(shí)域耦合分析,不能只考慮水動(dòng)力分析的位移RAO。

    表5 ML4最大有效張力無量綱比值

    5 結(jié) 論

    對(duì)浸沒式浮筒進(jìn)行水動(dòng)力分析和時(shí)域耦合分析,研究不同阻尼比、浮筒不同幾何形狀對(duì)浮筒水動(dòng)力性能和時(shí)域運(yùn)動(dòng)響應(yīng)以及系泊線張力的影響,并且開展浮筒幾何形狀的優(yōu)化設(shè)計(jì),主要得出如下結(jié)論:

    (1) 阻尼比越大,浸沒式浮筒在縱蕩、垂蕩和縱搖方向的運(yùn)動(dòng)RAO幅值越小,但由于考慮了系泊的作用,阻尼比對(duì)運(yùn)動(dòng)RAO的影響幅度很小。隨著阻尼比增大,浮筒時(shí)域運(yùn)動(dòng)響應(yīng)和系泊張力也逐漸減小。浮體的阻尼起到抑制浮體運(yùn)動(dòng)響應(yīng)、減小系泊線最大有效張力的作用,在實(shí)際工程中可考慮在浮筒上安裝阻尼板等結(jié)構(gòu)以增大浮體的阻尼。

    (2) 隨著浸沒式浮筒幾何直徑的增大:浮筒在縱蕩方向的附加質(zhì)量、阻尼減小,在垂蕩和縱搖方向的附加質(zhì)量、阻尼增大;浮筒在縱蕩和縱搖方向的位移RAO幅值減小,在垂蕩方向的位移RAO幅值增大;在3個(gè)自由度方向上的固有周期均增大。

    (3) 隨著浮筒幾何直徑的增大,浮筒越矮胖,時(shí)域動(dòng)態(tài)分析中浮筒在3個(gè)自由度方向的運(yùn)動(dòng)響應(yīng)和系泊線張力越大。對(duì)于該浸沒式浮筒,幾何形狀對(duì)浮筒頻域水動(dòng)力性能和時(shí)域運(yùn)動(dòng)響應(yīng)的影響趨勢不同。在頻域內(nèi),當(dāng)浮筒形狀為直徑8 m、9 m時(shí),浮筒縱搖方向RAO響應(yīng)幅值遠(yuǎn)大于其他形狀,顯然不合適;在時(shí)域內(nèi),浮筒形狀越細(xì)長、直徑越小,浮筒的時(shí)域運(yùn)動(dòng)響應(yīng)和系泊線的張力響應(yīng)越小。因此,浸沒式浮筒的選型設(shè)計(jì)需要綜合考量,在設(shè)計(jì)浸沒式浮筒形狀時(shí),通過水動(dòng)力分析的位移RAO來選型是基礎(chǔ)設(shè)計(jì),另外開展時(shí)域動(dòng)態(tài)分析非常有必要。

    猜你喜歡
    浮筒系泊阻尼比
    基于CFD方法的系泊狀態(tài)FSO舵載荷計(jì)算
    高韌性海洋用三級(jí)系泊鏈鋼的研制
    山東冶金(2022年2期)2022-08-08 01:50:50
    一種基于淺水靜水環(huán)境的浮筒構(gòu)架式水上鉆探平臺(tái)
    深水單點(diǎn)系泊維修中水浮筒保護(hù)性回收與下放工藝的應(yīng)用
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    黏滯阻尼器在時(shí)程分析下的附加有效阻尼比研究
    鋼質(zhì)浮筒在沉船人性化打撈中的應(yīng)用
    中國航海(2018年1期)2018-05-07 11:12:23
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    鋼制浮筒在“世越號(hào)”打撈工程中的應(yīng)用
    結(jié)構(gòu)構(gòu)件阻尼比對(duì)大跨度懸索橋地震響應(yīng)的影響
    中文字幕人妻丝袜制服| 不卡一级毛片| 亚洲精品久久成人aⅴ小说| 亚洲av成人av| 国产黄色免费在线视频| 性欧美人与动物交配| 国产激情欧美一区二区| 女人被躁到高潮嗷嗷叫费观| 久久中文看片网| 亚洲精品中文字幕一二三四区| 在线观看免费午夜福利视频| 黄片播放在线免费| 午夜激情av网站| 国产三级在线视频| 午夜福利免费观看在线| 亚洲色图 男人天堂 中文字幕| 国产单亲对白刺激| 久久国产亚洲av麻豆专区| 在线观看午夜福利视频| www.精华液| 咕卡用的链子| 新久久久久国产一级毛片| 免费高清视频大片| 国产精华一区二区三区| 欧美日韩av久久| 精品国产亚洲在线| 男女下面进入的视频免费午夜 | xxx96com| 国产乱人伦免费视频| 精品久久久久久久毛片微露脸| 91国产中文字幕| 国产激情久久老熟女| 亚洲欧美一区二区三区久久| 亚洲少妇的诱惑av| 天天添夜夜摸| 国产欧美日韩综合在线一区二区| 久久影院123| 色老头精品视频在线观看| 久久欧美精品欧美久久欧美| 成人免费观看视频高清| 很黄的视频免费| 黄片播放在线免费| 悠悠久久av| 亚洲精品一二三| 一个人观看的视频www高清免费观看 | 中国美女看黄片| 大码成人一级视频| 国产精品久久久久成人av| 久久人人97超碰香蕉20202| 人人妻人人爽人人添夜夜欢视频| 青草久久国产| 久久久久久久久免费视频了| 国产99白浆流出| 丰满的人妻完整版| 久久久久九九精品影院| 最近最新中文字幕大全电影3 | 精品久久蜜臀av无| 日韩免费高清中文字幕av| 999精品在线视频| 亚洲视频免费观看视频| 在线观看免费日韩欧美大片| 黑人欧美特级aaaaaa片| 亚洲一区二区三区色噜噜 | 老司机午夜十八禁免费视频| 91九色精品人成在线观看| 女人被躁到高潮嗷嗷叫费观| 黄色视频不卡| 日本wwww免费看| 精品国产乱子伦一区二区三区| 变态另类成人亚洲欧美熟女 | 99国产精品一区二区蜜桃av| 免费久久久久久久精品成人欧美视频| 精品福利观看| 国产主播在线观看一区二区| 亚洲精品中文字幕一二三四区| 性欧美人与动物交配| 99久久99久久久精品蜜桃| ponron亚洲| 中文字幕人妻丝袜一区二区| 欧美日韩福利视频一区二区| 在线观看舔阴道视频| 好男人电影高清在线观看| 国产午夜精品久久久久久| 日韩高清综合在线| 亚洲avbb在线观看| 亚洲色图av天堂| 国产又色又爽无遮挡免费看| 这个男人来自地球电影免费观看| 免费在线观看日本一区| ponron亚洲| 视频区欧美日本亚洲| 黄色片一级片一级黄色片| 窝窝影院91人妻| 国产精品免费视频内射| 纯流量卡能插随身wifi吗| 老熟妇仑乱视频hdxx| 人人澡人人妻人| 国产精品99久久99久久久不卡| 国产精品久久久久成人av| 亚洲av五月六月丁香网| av天堂在线播放| 无遮挡黄片免费观看| 国产成年人精品一区二区 | 国产欧美日韩一区二区三| 午夜亚洲福利在线播放| 超色免费av| 午夜久久久在线观看| 人妻久久中文字幕网| 搡老乐熟女国产| 一级毛片女人18水好多| 免费一级毛片在线播放高清视频 | 在线观看66精品国产| 亚洲国产欧美日韩在线播放| 精品国产一区二区三区四区第35| 亚洲精品美女久久av网站| 亚洲欧美一区二区三区黑人| 99在线视频只有这里精品首页| 中文字幕色久视频| 精品久久久久久久久久免费视频 | 久久国产精品男人的天堂亚洲| 在线看a的网站| 高清av免费在线| 手机成人av网站| 国产av一区二区精品久久| www.www免费av| 国产日韩一区二区三区精品不卡| 亚洲国产欧美日韩在线播放| 成人亚洲精品一区在线观看| 麻豆一二三区av精品| 18禁观看日本| 国产一区二区三区在线臀色熟女 | 久久天堂一区二区三区四区| 国产欧美日韩综合在线一区二区| 精品久久久精品久久久| 另类亚洲欧美激情| 99香蕉大伊视频| 亚洲av成人一区二区三| 亚洲片人在线观看| 精品免费久久久久久久清纯| 久久热在线av| 新久久久久国产一级毛片| 18禁美女被吸乳视频| 一级片'在线观看视频| 97人妻天天添夜夜摸| 他把我摸到了高潮在线观看| 国产av又大| 夜夜夜夜夜久久久久| 国产精品野战在线观看 | 亚洲国产欧美一区二区综合| 91老司机精品| 91成人精品电影| 亚洲精品中文字幕一二三四区| 欧洲精品卡2卡3卡4卡5卡区| 成人av一区二区三区在线看| 国产一区二区三区在线臀色熟女 | 美女高潮喷水抽搐中文字幕| 精品国产一区二区久久| 手机成人av网站| 亚洲激情在线av| 如日韩欧美国产精品一区二区三区| 免费观看精品视频网站| 免费看a级黄色片| 在线看a的网站| 免费av中文字幕在线| 亚洲精品中文字幕在线视频| 黄频高清免费视频| 国产成人av激情在线播放| 亚洲 国产 在线| 中文字幕人妻丝袜制服| 国产aⅴ精品一区二区三区波| 亚洲熟女毛片儿| 久久香蕉国产精品| 法律面前人人平等表现在哪些方面| 日韩中文字幕欧美一区二区| 国内毛片毛片毛片毛片毛片| 国产aⅴ精品一区二区三区波| 999久久久国产精品视频| 一夜夜www| 国产高清videossex| 亚洲国产欧美日韩在线播放| 免费女性裸体啪啪无遮挡网站| 在线播放国产精品三级| 亚洲精品成人av观看孕妇| 黄色视频不卡| 免费av中文字幕在线| 久久香蕉精品热| 国产午夜精品久久久久久| 91麻豆av在线| 天天躁夜夜躁狠狠躁躁| 久久精品91无色码中文字幕| 精品国产国语对白av| 日韩欧美一区视频在线观看| 久久性视频一级片| 又黄又爽又免费观看的视频| 久久久久久大精品| 91大片在线观看| 男人舔女人的私密视频| 视频在线观看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 国产高清激情床上av| 亚洲色图 男人天堂 中文字幕| 久久精品亚洲精品国产色婷小说| 老汉色av国产亚洲站长工具| 亚洲精品国产一区二区精华液| 欧美日韩黄片免| 欧美日韩亚洲综合一区二区三区_| 午夜久久久在线观看| 999久久久精品免费观看国产| 日韩视频一区二区在线观看| 国产精品一区二区免费欧美| 午夜福利欧美成人| 一进一出抽搐动态| 黑人巨大精品欧美一区二区mp4| 美女高潮到喷水免费观看| 熟女少妇亚洲综合色aaa.| 欧美亚洲日本最大视频资源| 日本a在线网址| 久久国产精品人妻蜜桃| 国产成人系列免费观看| 国产欧美日韩综合在线一区二区| 岛国在线观看网站| 淫秽高清视频在线观看| 9色porny在线观看| 一a级毛片在线观看| tocl精华| 国产精品99久久99久久久不卡| 国产三级黄色录像| 天堂影院成人在线观看| 一级a爱片免费观看的视频| 国产精华一区二区三区| 欧美激情久久久久久爽电影 | 午夜精品久久久久久毛片777| 精品高清国产在线一区| 在线观看免费视频网站a站| 午夜免费激情av| 好看av亚洲va欧美ⅴa在| 他把我摸到了高潮在线观看| 性欧美人与动物交配| 亚洲精品国产区一区二| 日韩人妻精品一区2区三区| xxxhd国产人妻xxx| 婷婷六月久久综合丁香| 老司机午夜福利在线观看视频| 激情视频va一区二区三区| 99精国产麻豆久久婷婷| 看免费av毛片| 少妇被粗大的猛进出69影院| 亚洲一区高清亚洲精品| 纯流量卡能插随身wifi吗| 日本a在线网址| 国产精品电影一区二区三区| 黄色女人牲交| 精品无人区乱码1区二区| 麻豆成人av在线观看| 久久精品影院6| 婷婷丁香在线五月| 悠悠久久av| 交换朋友夫妻互换小说| 亚洲 欧美一区二区三区| 精品一品国产午夜福利视频| 免费在线观看黄色视频的| 97人妻天天添夜夜摸| 国产成人精品在线电影| 人人妻人人添人人爽欧美一区卜| 午夜福利一区二区在线看| aaaaa片日本免费| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 亚洲av片天天在线观看| 天堂√8在线中文| 校园春色视频在线观看| 乱人伦中国视频| 国产成人影院久久av| 亚洲一区二区三区不卡视频| 91精品三级在线观看| 成人18禁在线播放| 在线观看一区二区三区| 久久草成人影院| 纯流量卡能插随身wifi吗| 丝袜人妻中文字幕| 成年版毛片免费区| 男女之事视频高清在线观看| 亚洲午夜理论影院| 国产成人精品久久二区二区91| 中文字幕高清在线视频| 欧美日本中文国产一区发布| 国产精品一区二区在线不卡| 国产亚洲精品第一综合不卡| 免费在线观看亚洲国产| 日本免费一区二区三区高清不卡 | 最新美女视频免费是黄的| 美女大奶头视频| 日韩有码中文字幕| 麻豆成人av在线观看| 性色av乱码一区二区三区2| 中文字幕av电影在线播放| 午夜福利一区二区在线看| 精品无人区乱码1区二区| 国产一区二区三区视频了| 一级片免费观看大全| 亚洲午夜理论影院| 成人免费观看视频高清| 脱女人内裤的视频| 长腿黑丝高跟| 制服人妻中文乱码| 久久精品国产清高在天天线| 亚洲精品国产一区二区精华液| 神马国产精品三级电影在线观看 | 91成人精品电影| 欧美日韩福利视频一区二区| 女人被躁到高潮嗷嗷叫费观| 久久久国产成人免费| 久久久久国产精品人妻aⅴ院| 一个人免费在线观看的高清视频| 看片在线看免费视频| 一进一出好大好爽视频| www.精华液| 国产在线观看jvid| 欧美中文日本在线观看视频| 国产男靠女视频免费网站| 久久久国产成人免费| 中文亚洲av片在线观看爽| 久久久久久久精品吃奶| 两性午夜刺激爽爽歪歪视频在线观看 | 国产激情久久老熟女| 最好的美女福利视频网| 日本wwww免费看| 麻豆久久精品国产亚洲av | 亚洲精品一区av在线观看| 中文欧美无线码| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| netflix在线观看网站| 黄色视频不卡| 日日夜夜操网爽| 国产欧美日韩一区二区三| 亚洲国产中文字幕在线视频| 欧美在线黄色| 中文字幕色久视频| 国产一区在线观看成人免费| 精品久久蜜臀av无| 男女下面进入的视频免费午夜 | 亚洲国产精品一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 久久久国产精品麻豆| 夜夜躁狠狠躁天天躁| 亚洲第一青青草原| 国产aⅴ精品一区二区三区波| av欧美777| 婷婷丁香在线五月| 国产不卡一卡二| 久久久国产成人精品二区 | 黄色丝袜av网址大全| 午夜福利欧美成人| 亚洲国产中文字幕在线视频| 男女做爰动态图高潮gif福利片 | 一级毛片精品| 亚洲熟妇中文字幕五十中出 | 久久久国产精品麻豆| 人妻丰满熟妇av一区二区三区| 18禁国产床啪视频网站| 黄色a级毛片大全视频| 在线永久观看黄色视频| 亚洲国产中文字幕在线视频| 国产精品国产av在线观看| 精品午夜福利视频在线观看一区| 波多野结衣av一区二区av| 欧美日韩福利视频一区二区| 亚洲av成人不卡在线观看播放网| 99热国产这里只有精品6| 亚洲精品国产区一区二| 女警被强在线播放| 亚洲av成人不卡在线观看播放网| 亚洲欧美激情综合另类| 村上凉子中文字幕在线| 热99re8久久精品国产| 三上悠亚av全集在线观看| 成人黄色视频免费在线看| 亚洲中文字幕日韩| 99riav亚洲国产免费| 麻豆成人av在线观看| a在线观看视频网站| 国产精品免费视频内射| 这个男人来自地球电影免费观看| 欧美日韩黄片免| 91国产中文字幕| 亚洲成人免费电影在线观看| 欧美日韩视频精品一区| 黄色 视频免费看| 亚洲色图 男人天堂 中文字幕| 麻豆久久精品国产亚洲av | 国产在线观看jvid| 高清av免费在线| 在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 80岁老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 国产成年人精品一区二区 | 久久人妻熟女aⅴ| 亚洲专区字幕在线| 69av精品久久久久久| 久久精品亚洲精品国产色婷小说| 亚洲一码二码三码区别大吗| 国产精品免费视频内射| 十八禁人妻一区二区| 少妇 在线观看| 国产成人影院久久av| 一边摸一边做爽爽视频免费| 一区二区三区精品91| 一边摸一边做爽爽视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 国产 在线| 午夜老司机福利片| 日本欧美视频一区| 亚洲精品国产一区二区精华液| 日韩欧美一区视频在线观看| 亚洲在线自拍视频| 免费高清视频大片| 国产成人精品久久二区二区免费| 日本a在线网址| 亚洲国产精品合色在线| 十八禁人妻一区二区| 精品国产超薄肉色丝袜足j| 欧美乱妇无乱码| 国产精品综合久久久久久久免费 | 校园春色视频在线观看| 少妇被粗大的猛进出69影院| 久久久久久久久免费视频了| 淫秽高清视频在线观看| 成人永久免费在线观看视频| 满18在线观看网站| 一边摸一边做爽爽视频免费| 后天国语完整版免费观看| 黄色 视频免费看| x7x7x7水蜜桃| 美女大奶头视频| 亚洲avbb在线观看| 男女床上黄色一级片免费看| 看片在线看免费视频| 精品国产乱码久久久久久男人| 色婷婷久久久亚洲欧美| 1024视频免费在线观看| 丝袜在线中文字幕| 很黄的视频免费| 亚洲美女黄片视频| www.熟女人妻精品国产| 岛国视频午夜一区免费看| 免费在线观看日本一区| 91av网站免费观看| 老司机在亚洲福利影院| 一区二区三区精品91| 亚洲专区中文字幕在线| 一级毛片精品| 久久精品人人爽人人爽视色| 无限看片的www在线观看| 国产成人精品无人区| 777久久人妻少妇嫩草av网站| 久久久久九九精品影院| 天天添夜夜摸| 国产精品99久久99久久久不卡| a在线观看视频网站| 999久久久精品免费观看国产| 精品国产一区二区久久| 天天影视国产精品| 国产aⅴ精品一区二区三区波| 国产一区二区三区综合在线观看| 日本三级黄在线观看| 18禁国产床啪视频网站| 国产精品野战在线观看 | 又紧又爽又黄一区二区| 成人黄色视频免费在线看| 国产深夜福利视频在线观看| 一级片'在线观看视频| 免费在线观看完整版高清| 99国产精品99久久久久| 日韩av在线大香蕉| 999久久久精品免费观看国产| 很黄的视频免费| 美女午夜性视频免费| 成人三级做爰电影| 亚洲精品一二三| 黄色毛片三级朝国网站| 天堂中文最新版在线下载| 国产精品1区2区在线观看.| tocl精华| 国产三级在线视频| 在线播放国产精品三级| videosex国产| 亚洲国产精品合色在线| tocl精华| 久久久国产欧美日韩av| 久久伊人香网站| 亚洲精品美女久久av网站| 黑人操中国人逼视频| 亚洲aⅴ乱码一区二区在线播放 | 人妻久久中文字幕网| 久久中文看片网| 国产99久久九九免费精品| 99精国产麻豆久久婷婷| 亚洲人成网站在线播放欧美日韩| 国产成人一区二区三区免费视频网站| 午夜a级毛片| 国产精品一区二区免费欧美| 亚洲性夜色夜夜综合| 侵犯人妻中文字幕一二三四区| 中文欧美无线码| 午夜免费成人在线视频| 高清在线国产一区| 亚洲中文av在线| 日韩精品中文字幕看吧| 青草久久国产| 国产亚洲精品第一综合不卡| 最近最新中文字幕大全电影3 | 丁香欧美五月| xxx96com| 午夜亚洲福利在线播放| 中文亚洲av片在线观看爽| 免费不卡黄色视频| 久久精品aⅴ一区二区三区四区| 日日爽夜夜爽网站| 国产成人精品无人区| 精品国产超薄肉色丝袜足j| 最新在线观看一区二区三区| 久久精品国产亚洲av香蕉五月| 人人妻人人添人人爽欧美一区卜| 久久精品国产亚洲av高清一级| 亚洲成人免费电影在线观看| 欧美黄色淫秽网站| 少妇被粗大的猛进出69影院| 国产成人免费无遮挡视频| 黄频高清免费视频| 国产乱人伦免费视频| 99久久国产精品久久久| 精品国产一区二区三区四区第35| 日韩大尺度精品在线看网址 | 精品午夜福利视频在线观看一区| 亚洲av美国av| 一本大道久久a久久精品| 日韩高清综合在线| av中文乱码字幕在线| 国产aⅴ精品一区二区三区波| 露出奶头的视频| 变态另类成人亚洲欧美熟女 | 美女福利国产在线| 人成视频在线观看免费观看| 亚洲 国产 在线| 亚洲精品国产一区二区精华液| 欧美日本中文国产一区发布| 成年版毛片免费区| 亚洲人成电影免费在线| 80岁老熟妇乱子伦牲交| 中文字幕高清在线视频| 99riav亚洲国产免费| av中文乱码字幕在线| 色尼玛亚洲综合影院| 嫩草影院精品99| 久热爱精品视频在线9| 精品久久久久久电影网| 久久精品国产亚洲av香蕉五月| 美女福利国产在线| a在线观看视频网站| 亚洲伊人色综图| 最近最新中文字幕大全电影3 | 在线观看一区二区三区激情| 人人澡人人妻人| 高清黄色对白视频在线免费看| 国产一区二区三区视频了| 久久久国产一区二区| 国产一区二区三区综合在线观看| 色老头精品视频在线观看| 国产国语露脸激情在线看| 男女做爰动态图高潮gif福利片 | 欧美成人午夜精品| 亚洲欧美激情综合另类| 精品国产美女av久久久久小说| 欧美乱码精品一区二区三区| 18禁裸乳无遮挡免费网站照片 | 国产亚洲精品一区二区www| 十八禁人妻一区二区| 无遮挡黄片免费观看| 黄色a级毛片大全视频| 免费人成视频x8x8入口观看| 老司机靠b影院| 欧美中文综合在线视频| 亚洲中文av在线| 侵犯人妻中文字幕一二三四区| 成人国语在线视频| 亚洲第一青青草原| 久久久久国产精品人妻aⅴ院| 亚洲一卡2卡3卡4卡5卡精品中文| 久热这里只有精品99| 午夜日韩欧美国产| 在线观看免费日韩欧美大片| 精品一区二区三卡| 中文字幕高清在线视频| 国产色视频综合| 免费少妇av软件| 丝袜美足系列| 国产黄色免费在线视频| 免费少妇av软件| 国产精品野战在线观看 | av中文乱码字幕在线| 国产精品一区二区免费欧美| 国产成人精品久久二区二区91| 久久人人97超碰香蕉20202| 久久热在线av| 亚洲七黄色美女视频| 国产深夜福利视频在线观看| 亚洲精品国产色婷婷电影| 亚洲av第一区精品v没综合| 成人av一区二区三区在线看| 天堂影院成人在线观看| 欧美日韩乱码在线| 国产高清视频在线播放一区| 中文字幕最新亚洲高清|