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

    螺旋側板的導程對VIV影響的數(shù)值模擬

    2016-01-08 05:41:05王亞非孫麗萍吳梓鑫
    關鍵詞:導程數(shù)值模擬

    王亞非 孫麗萍 吳梓鑫

    (哈爾濱工程大學深海工程技術研究中心1) 哈爾濱 150001) (江蘇科技大學船舶與海洋工程學院2) 鎮(zhèn)江 212000)

    螺旋側板的導程對VIV影響的數(shù)值模擬

    王亞非1)孫麗萍1)吳梓鑫2)

    (哈爾濱工程大學深海工程技術研究中心1)哈爾濱150001)(江蘇科技大學船舶與海洋工程學院2)鎮(zhèn)江212000)

    摘要:海洋立管在來流影響下極易發(fā)生渦激振動(VIV),VIV會對立管產(chǎn)生嚴重的疲勞破壞作用,導致立管很快失效.之前有大量的研究主要集中在如何采用有效措施抑制其影響,在工程中螺旋側板的使用居多,螺旋側板在設計中有很多參數(shù),如螺距、側板高度和覆蓋率等.文中基于Fluent軟件對不同導程螺旋側板的立管的尾流場進行數(shù)值模擬.結果發(fā)現(xiàn),螺旋側板可以明顯的抑制渦激振動的影響.同時得到:在研究范圍內,螺旋側板隨導程增大控制旋渦發(fā)放效果越好,從而抑制渦激振動效果越佳;將數(shù)值模擬結果和文獻的實驗結果相比較,發(fā)現(xiàn)二者吻合良好,驗證了數(shù)值模擬的可靠性.

    關鍵詞:螺旋側板;導程;VIV;數(shù)值模擬

    王亞非(1989- ):男,碩士生,主要研究領域為海洋立管

    0引言

    海洋立管是連接海底井口與海洋平臺或鉆井船的重要構件,既可以傳輸油氣,又可以用于鉆探生產(chǎn).VIV會對其產(chǎn)生嚴重的結構疲勞破壞作用[1],減少了立管的使用壽命.目前人們研制出了各種渦激抑制裝置,主要有螺旋側板(helicalstrakes)、導流板(splitterplate)、屏蔽裝置(shroud),以及減振器(fairings)等,由于加工及安裝簡單、生產(chǎn)成本低等原因,實際工程中大多采用螺旋側板作為抑制設備,實驗也證明了其有效的抑制效果[2-4].螺旋側板的設計參數(shù)很多,主要是螺距、側板高度[5-6].

    Korkischko等[7]在水動力實驗室中通過改變螺旋側板的幾何參數(shù)研究了對VIV的抑制效果,最后得到結論:螺旋側板可以改變旋渦脫落的頻率和降低立管的振動幅值.T.Zhou等[8]在UniversityofWesternAustralia風洞中用安裝了3個導程P=10 D和高度H=0.12 D螺旋側板立管進行實驗,實驗結果發(fā)現(xiàn)裸管在約化速度5~8.5范圍內有明顯的自鎖現(xiàn)象,而螺旋側板可以很好的抑制旋渦脫落同時立管沒有發(fā)生自鎖現(xiàn)象.Lee等[9]在OsakaUniversity的實驗水池用模型研究發(fā)現(xiàn)當螺旋側板的P=10 D和H=0.1 D時,側板可以明顯的減小VIV對柔性立管的水動力作用,升力大小對螺旋側板高度比較敏感.

    但之前的研究多數(shù)集中于實驗方法,較少的涉及對螺旋側板的導程和VIV抑制作用的關系的進行控制變量的數(shù)值分析.故本文的目的是運用Fluent對螺旋側板的導程大小和VIV作用之間關系進行數(shù)值模擬,重點研究了立管尾流場中的旋渦脫落和立管所受的水動力隨導程變化情況.在數(shù)值模擬中,使用高度H=0.08 D(D=0.02m)和導程P分別為5D,8D,15D,17D的螺旋側板進行研究,將結果和文獻的實驗數(shù)據(jù)對比,發(fā)現(xiàn)二者吻合良好,表明了數(shù)值模擬的可靠性.

    1數(shù)學模型

    1.1流體力學控制方程

    當流體運動處于湍流狀態(tài)時,由非穩(wěn)態(tài)的連續(xù)方程和N-S方程而求解得到.

    (1)

    式中:F為質量力;ρ為密度;t為時間;p為壓力;μ為流體的動力粘性系數(shù).

    (2)

    1.2結構動力學方程

    圓柱體在運動的流體中因為所受的渦激力產(chǎn)生流向和橫向的運動,其運動方程為

    (3)

    (4)

    式中:m為單位長度圓柱體的質量;c為結構阻尼系數(shù);k為彈簧剛度系數(shù);Fx(t)為流向的阻力;Fy(t)為橫向的升力.

    1.3LES原理

    (5)

    式中:u為速度矢量在X,Y,Z方向的分量;p為壓力;ρ為密度;τij為亞格子雷諾應力;μ為流體的動力粘性系數(shù).

    2數(shù)值模擬

    參考之前的實驗研究和模擬計算[10],該算例流場區(qū)域采用長30D,寬20D,高15D的長方體,D=20mm,D為立管模型直徑,上游為10D,下游為20D,兩邊寬各為10D.考慮到不同區(qū)域網(wǎng)格對流場變化的影響,以立管模型為中心,周圍8D范圍內網(wǎng)格加密.與立管壁面相鄰區(qū)域也是流動狀態(tài)變化較為劇烈的區(qū)域,對整個流場的影響較大,也需進行加密處理.光滑立管模型一共分為12個區(qū)域,均采用結構網(wǎng)格;帶螺旋側板立管模型一共分為9個區(qū)域,進口、出口以及壁面區(qū)各分為3個區(qū)域,周圍八個區(qū)域采用結構網(wǎng)格,中間區(qū)域采用非結構網(wǎng)格,具體平面見圖1,立體見圖2.圖3為帶螺旋側板立管的表邊網(wǎng)格,圖4為帶不同導程螺旋側板立管的物理模型.

    圖1 二維模型平面圖

    圖2 三維計算區(qū)域網(wǎng)格劃分示意圖

    圖4 帶螺旋側板立管的物理模型

    圖3帶螺旋側板立管表面網(wǎng)格劃分

    流體計算域的設置對于流體計算域所進行的簡化設定為:

    1) 參考物理模型中,立管系統(tǒng)完全處于自由液面以下的水體中心,不考慮空氣影響,采取單相模擬.

    2) 流體計算域的左側設為速度入口(velocityinlet);側面和上下面均采用無滑移條件(wall),右側設為自由出流邊界(outflow).

    3結果與討論

    3.1光滑立管

    將網(wǎng)格劃分信息輸入Fluent,在Re=5 000時使用LES模型進行求解.得到拖拽力系數(shù)Cd和升力系數(shù)Cl曲線圖.

    可以求出拖拽力系數(shù)約為1.27,與文獻[11]的模擬結果1.26吻合良好.運用MATLAB軟件通過編程做出頻率分析圖(見圖5),得到該數(shù)值模擬情況下主頻為2.560 8Hz,由此算得斯特勞哈數(shù)為0.205,與文獻[12]中的計算結果0.225相接近,說明了計算模型的可靠性,因而可以將所用區(qū)域及其他設置模擬帶螺旋側板的立管,并可以將該模擬結果作為帶螺旋側板立管的模擬結果的對比.

    圖5 頻率分析圖

    3.2帶不同導程的螺旋側板立管

    為了和光滑立管相對比,同時參考文獻[8-9]實驗研究,此算例采取側板高為0.08D,選取不同導程(5D,8D,15D,17D)進行模擬計算,記為H0.08D-P5,H0.08D-P8,H0.08D-P15,以及H0.08D-P17四種螺旋側板.在Fluent軟件中帶螺旋側板的立管參數(shù)設置與光滑立管相同.通過數(shù)值模擬得到4種螺旋側板與光滑立管的拖拽力系數(shù)圖.

    圖6 四種螺旋側板立管與光滑立管的Cd曲線

    從圖6可以清晰的看出加上側板高為0.08D的螺旋側板后,隨著螺旋側板導程的變化,旋渦瀉放的頻率發(fā)生很大變化,拖拽力系數(shù)總體上介于0.9~1.05之間,相比于光滑立管有了顯著地下降,即加上螺旋側板后立管所受拖拽力減小了,立管所受的損壞也就減輕了.除此而外,還可以看出加上螺旋側板后,拖拽力系數(shù)的振幅有了明顯的減小,同樣說明了立管受力的減弱.通過數(shù)據(jù)處理,得到表1,針對拖拽力系數(shù)分析導程的影響情況.

    表1 光滑立管和相同側板高不同導程螺旋側板拖拽力系數(shù)

    分析表1可知,加上螺旋側板后立管的拖拽力系數(shù)平均減少了20%之多,極差不及原先光滑立管的10%,通過數(shù)據(jù)量化的說明了螺旋側板減少對立管損傷的優(yōu)越性.對比不同導程的螺旋側板,可以看到當螺旋側板高度一定時,在5D~15D導程之內,拖拽力系數(shù)及其極差隨著導程增大而變大,而導程為17D時,拖拽力系數(shù)與其極差小于15D導程值.

    雖然入口邊界來流是單向的,但是在繞流過螺旋側板后,旋渦流由端部向尾流移動,流體呈現(xiàn)出復雜的三維特性.文中選用Z向渦量云圖比較其尾流場.

    圖7為光滑立管在Re=5 000時流場內部渦瀉放Z向渦量云圖,可以看出光滑立管在整個流場中呈現(xiàn)出較為規(guī)則的渦街發(fā)放形態(tài).圖8~11為0.08D側板高時不同導程下Z向渦量情況,可以看出螺旋側板影響了渦的發(fā)放,并且隨著導程的增大,立管后方渦的瀉放越加紊亂,渦街已不能有規(guī)律的發(fā)放,表明螺旋側板很好的抑制了渦激振動的影響.同時也證實了螺旋側板導程越大,對尾流場影響越大,進而說明了能夠更好地抑制渦激振動.

    圖7 光滑立管Z向渦量云圖

    圖8 H0.08-P5立管Z向渦量云圖

    圖9 H0.08-P8立管Z向渦量云圖

    圖10 H0.08-P15立管Z向渦量云圖

    圖11 H0.08-P17立管Z向渦量云圖

    4結束語

    本文利用大渦模擬數(shù)值計算了三維情況下雷諾數(shù)Re=5 000時光滑立管和帶有螺旋側板的立管,分別通過選用不同導程研究了螺旋側板對尾流場的影響.通過拖拽力系數(shù)、升力系數(shù)以及Z向渦量云圖的對比分析,得出結論:在模擬范圍內螺旋側板隨導程增大明顯的抑制了渦激振動對立管的影響,其中導程P=17D綜合效果最佳,且水動力系數(shù)減小幅度達到了23.93%.將數(shù)值模擬和文獻[3]的實驗結果(水動力減額為25%)相對比,發(fā)現(xiàn)二者的水動力系數(shù)減幅基本接近,同時和DNV-RP-F204推薦導程值(17.5D)相符合.

    本文在研究螺旋側板抑制VIV時僅考慮了導程對其影響,但實際主要影響因素還有螺旋擦板的高度,建議以后在相關的研究中考慮高度的作用效果,這樣可以更好的探究螺旋側板抑制渦激振動的本質.

    參 考 文 獻

    [1]HOWELLSH,BOWMANJ.Drillingriser/wellsysteminteractionindeepwater,harshenvironments[R].2HOffshoreEngineeringLimited,1997.

    [2]NAOAKIM,AKIYOSHIB,KOJIT.Experimentalstudyofaflexibleriseroveredwithhelicalstrakes[C].TheThirteenthInternatinalOffshoreAndPolarEngineeringConference,2003.

    [3]HUANGShan.VIVsuppressionofatwo-degree-of-freedomcircularcylinderanddragreductionofafixedcircularcylinderbytheuseofhelicalgrooves[J].JournalofFluidsandStructures,2011,27:1124-1133.

    [4]BLUMBERGM,TELLIERE.Experimentalevaluationofvortexinducedvibrationresponseofstrakedpipesintandemarrangements[C].OMAE2012-83772.

    [5]沙勇,曹靜,張恩勇,等.抑制渦激振動的螺旋列板設計參數(shù)研究[J]. 海洋工程,2013, 31(1): 43-48.

    [6]魏澤,吳乘勝,倪陽.基于正交設計和CFD模擬的Spar平臺螺旋側板水動力優(yōu)化設計研究[J]. 船舶力學, 2013, 17(10): 1133-1139.

    [7]KORKISCHKOI,MENEGHINIJR.Experimentalinvestigationofflow-inducedvibrationonisolatedandtandemcircularcylindersfittedwithstrakes[J].JournalofFluidsandStructures,2010,26:611-625.

    [8]ZHOUT,MOHDSF.RAZALI,Z,etal.Onthestudyofvortex-inducedvibrationofacylinderwithhelicalstrakes[J].JournalofFluidsandStructures,2011,27:903-917.

    [9]LEEKQ,AMINUDINA,NAOMIK,etal.InvestigationontheeffectivenessofhelicalstrakesinsuppressingVIVofflexibleriser[J].AppliedOceanResearch,2014 ,44:82-91.

    [10]TEZDUYARTE,SHIHR.Numericalexperimentsondownstreamboundaryofflowpastcylinder[J].J.Engrg.Mech.1991,117:854-871.

    [11]詹昊,李萬平,方秦漢,等.不同雷諾數(shù)下圓柱繞流仿真計算[J].武漢理工大學學報,2008,30(12):129-132.

    [12]WILLIAMSONCHK.Vortexdynamicsinthecylinderwake[J].AnnualReviewofFluidMechanics,1996,28:477-539.

    NumericalSimulationontheEffectivenessofthe

    中圖法分類號:U661.1

    doi:10.3963/j.issn.2095-3844.2015.01.040

    收稿日期:2014-11-16

    PitchofHelicalStrakesinSuppressingVIV
    WANGYafei1)SUNLiping1)WUZixin2)
    (Deepwater Engineering Technology Research Center,

    Harbin Engineering University, Harbin 150001, China)1)

    (College of Navel Architecture and Marine Engineering,

    Jiangsu University of Science and Technology, Zhenjiang 212000, China)2)

    Abstract:Due to the serious fatigue damage of VIV (vortex-induced vibration) to the riser, it is very likely to failure soon. Much corresponding research has been done to deal with how to suppress the effect recently, and the helical strakes is widely used in engineering now. Many parameters have certain effects on its design in practice, including the height, pitch and coverage scale. The main aim of this paper is to study the wake field of a riser with different pitch by numerical simulation based on the Fluent software. In conclusion, the VIV can be obviously suppressed by the helical strakes, and the effectiveness of inhibiting vortex shedding is better with increasing the pitch of helical strakes. Besides, the numerical simulation has a good agreement with the experiment result of reference, which demonstrates the accuracy of research approach.

    Key words:helical strakes; pitch; VIV; numerical simulation

    猜你喜歡
    導程數(shù)值模擬
    雙螺母滾珠絲杠副導程誤差與摩擦力矩關系研究
    車削變齒厚蝸桿時交換齒輪齒數(shù)的計算方法及車削時應注意的幾個問題
    論四頭雙導程蝸桿車削掛輪的選配
    論四頭雙導程蝸桿車削掛輪的選配
    普通車床車削多線螺紋的工藝探究
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    尾随美女入室| 国产成年人精品一区二区| 国产一级毛片七仙女欲春2| 欧美三级亚洲精品| 国产午夜福利久久久久久| 免费无遮挡裸体视频| 成人亚洲精品一区在线观看 | 欧美高清成人免费视频www| 天堂影院成人在线观看| av黄色大香蕉| 中文在线观看免费www的网站| 2018国产大陆天天弄谢| 蜜臀久久99精品久久宅男| 日韩欧美 国产精品| 女人久久www免费人成看片| 国产激情偷乱视频一区二区| 晚上一个人看的免费电影| 午夜激情欧美在线| 亚洲综合精品二区| 日日干狠狠操夜夜爽| 久久久久久久午夜电影| 99久久中文字幕三级久久日本| 欧美日韩在线观看h| 国产精品久久久久久久电影| 午夜免费观看性视频| 国产一区有黄有色的免费视频 | 免费观看的影片在线观看| 精品久久久久久久久av| 99热网站在线观看| 国产精品无大码| 在线观看av片永久免费下载| 日本免费在线观看一区| 国产精品一区二区三区四区免费观看| 日韩电影二区| 又爽又黄a免费视频| 亚洲av电影不卡..在线观看| 麻豆乱淫一区二区| 国产精品一区www在线观看| 一级二级三级毛片免费看| 街头女战士在线观看网站| 小蜜桃在线观看免费完整版高清| 亚洲图色成人| 欧美3d第一页| 2021少妇久久久久久久久久久| av卡一久久| 色综合站精品国产| 亚州av有码| 在线观看一区二区三区| 天天躁夜夜躁狠狠久久av| 高清午夜精品一区二区三区| 不卡视频在线观看欧美| 久久久久久久久中文| 青春草视频在线免费观看| 亚洲在线自拍视频| 精品久久久久久成人av| 女人被狂操c到高潮| 亚洲四区av| 亚洲精品影视一区二区三区av| 超碰av人人做人人爽久久| 欧美性猛交╳xxx乱大交人| 日韩 亚洲 欧美在线| 亚洲伊人久久精品综合| 日韩中字成人| 国产精品无大码| 免费无遮挡裸体视频| 国产探花在线观看一区二区| 老师上课跳d突然被开到最大视频| 蜜桃久久精品国产亚洲av| 亚洲色图av天堂| 成人综合一区亚洲| 国产黄片视频在线免费观看| 亚洲精品成人av观看孕妇| 黑人高潮一二区| 99视频精品全部免费 在线| 亚洲人成网站在线观看播放| 日本色播在线视频| 亚洲美女搞黄在线观看| 婷婷色麻豆天堂久久| 高清视频免费观看一区二区 | 亚洲av成人精品一区久久| 蜜桃久久精品国产亚洲av| 熟妇人妻不卡中文字幕| 午夜福利成人在线免费观看| 久久精品国产鲁丝片午夜精品| 国产精品99久久久久久久久| 国产黄片美女视频| 亚洲国产精品专区欧美| 你懂的网址亚洲精品在线观看| 亚洲精品亚洲一区二区| 色综合亚洲欧美另类图片| 综合色av麻豆| 成年人午夜在线观看视频 | 高清视频免费观看一区二区 | 亚洲综合精品二区| 国产成人精品一,二区| 国产av码专区亚洲av| 一个人看视频在线观看www免费| 汤姆久久久久久久影院中文字幕 | 男女啪啪激烈高潮av片| 欧美精品一区二区大全| 久久久午夜欧美精品| av国产免费在线观看| 99热网站在线观看| 有码 亚洲区| 国精品久久久久久国模美| 欧美zozozo另类| 2021少妇久久久久久久久久久| 国语对白做爰xxxⅹ性视频网站| 老司机影院成人| 97精品久久久久久久久久精品| 国产亚洲一区二区精品| 国产精品伦人一区二区| 老女人水多毛片| av专区在线播放| 亚洲av中文字字幕乱码综合| 欧美+日韩+精品| av在线蜜桃| 91av网一区二区| 亚洲18禁久久av| 亚洲av一区综合| 午夜福利视频1000在线观看| 最近的中文字幕免费完整| 蜜桃亚洲精品一区二区三区| 最近最新中文字幕大全电影3| 日韩电影二区| 国产成人一区二区在线| 国产男女超爽视频在线观看| 国产乱人视频| 亚洲精品成人久久久久久| 日韩大片免费观看网站| 日韩av免费高清视频| 久久97久久精品| 成人av在线播放网站| 99久国产av精品| 久久久久久久亚洲中文字幕| av在线亚洲专区| 久久久久久久久久人人人人人人| 国产高潮美女av| 久久精品夜夜夜夜夜久久蜜豆| 又粗又硬又长又爽又黄的视频| 成人特级av手机在线观看| 卡戴珊不雅视频在线播放| 97精品久久久久久久久久精品| 在现免费观看毛片| 成人鲁丝片一二三区免费| 别揉我奶头 嗯啊视频| 一级毛片aaaaaa免费看小| 女人久久www免费人成看片| 国产精品久久久久久久电影| 中文字幕av成人在线电影| 黄片wwwwww| 国产亚洲91精品色在线| 日韩av在线大香蕉| 亚洲人与动物交配视频| 直男gayav资源| 天堂av国产一区二区熟女人妻| 黄色一级大片看看| 免费观看在线日韩| 男人舔女人下体高潮全视频| 国国产精品蜜臀av免费| videossex国产| 青春草国产在线视频| 成人综合一区亚洲| 国产精品一区二区三区四区免费观看| 男女啪啪激烈高潮av片| 我的女老师完整版在线观看| 嫩草影院新地址| 国产精品美女特级片免费视频播放器| 亚洲不卡免费看| 可以在线观看毛片的网站| 欧美不卡视频在线免费观看| 国内少妇人妻偷人精品xxx网站| 男人舔女人下体高潮全视频| 综合色av麻豆| 毛片女人毛片| 欧美极品一区二区三区四区| 97超碰精品成人国产| 日本-黄色视频高清免费观看| 精品人妻视频免费看| 亚洲av成人av| 免费观看无遮挡的男女| 能在线免费观看的黄片| 丝瓜视频免费看黄片| 午夜视频国产福利| 三级男女做爰猛烈吃奶摸视频| 中文欧美无线码| 国产黄a三级三级三级人| 久久午夜福利片| 婷婷色av中文字幕| 真实男女啪啪啪动态图| 大香蕉久久网| 婷婷色麻豆天堂久久| 在线免费观看的www视频| 91在线精品国自产拍蜜月| 久久人人爽人人片av| 丝袜喷水一区| 成人高潮视频无遮挡免费网站| 亚洲激情五月婷婷啪啪| 国产真实伦视频高清在线观看| 2018国产大陆天天弄谢| 国产成人a∨麻豆精品| 建设人人有责人人尽责人人享有的 | 国产黄色小视频在线观看| 最近手机中文字幕大全| 国产精品久久久久久久久免| 欧美97在线视频| 精品不卡国产一区二区三区| 日韩欧美精品v在线| 国产视频首页在线观看| 久久久久网色| 禁无遮挡网站| 永久免费av网站大全| 精品久久久久久久末码| 国产成人午夜福利电影在线观看| 国产精品人妻久久久影院| 亚洲欧美精品专区久久| 日本午夜av视频| 久久久精品94久久精品| 国产老妇伦熟女老妇高清| 国产高清三级在线| 国产高清不卡午夜福利| 精品国产三级普通话版| 亚洲精品久久午夜乱码| 国产精品国产三级国产专区5o| 国产在线男女| 国产免费一级a男人的天堂| 男女下面进入的视频免费午夜| av在线天堂中文字幕| 一级黄片播放器| 日日啪夜夜爽| 99久久九九国产精品国产免费| 2021天堂中文幕一二区在线观| 色哟哟·www| 五月天丁香电影| 亚洲在线自拍视频| 九九在线视频观看精品| 国产三级在线视频| 听说在线观看完整版免费高清| 水蜜桃什么品种好| 亚洲三级黄色毛片| av.在线天堂| 久久久成人免费电影| 久久99热这里只频精品6学生| a级毛色黄片| 欧美97在线视频| 国产成人精品久久久久久| 亚洲av男天堂| 丰满少妇做爰视频| 最近最新中文字幕大全电影3| 久久久久久久久久人人人人人人| 国产成人91sexporn| 国产精品久久视频播放| 直男gayav资源| 波野结衣二区三区在线| 久久精品国产自在天天线| 美女大奶头视频| 九色成人免费人妻av| 简卡轻食公司| a级一级毛片免费在线观看| 成人毛片60女人毛片免费| 国产精品一区二区三区四区久久| 日韩欧美精品免费久久| 蜜桃亚洲精品一区二区三区| 小蜜桃在线观看免费完整版高清| 国产免费福利视频在线观看| 五月玫瑰六月丁香| 禁无遮挡网站| 亚洲人成网站高清观看| av在线蜜桃| 一级片'在线观看视频| 美女cb高潮喷水在线观看| 免费黄网站久久成人精品| 日本熟妇午夜| 亚洲真实伦在线观看| 国产一区二区亚洲精品在线观看| 男人爽女人下面视频在线观看| 联通29元200g的流量卡| 日韩欧美 国产精品| 男女边吃奶边做爰视频| a级毛色黄片| 一个人看视频在线观看www免费| 精品久久久久久久久亚洲| 欧美zozozo另类| 免费黄网站久久成人精品| 国产精品.久久久| 少妇丰满av| 18禁裸乳无遮挡免费网站照片| 亚洲av成人精品一区久久| 日韩一本色道免费dvd| 日韩强制内射视频| 国产麻豆成人av免费视频| 国产亚洲精品av在线| 国产三级在线视频| 淫秽高清视频在线观看| 最近2019中文字幕mv第一页| 午夜免费男女啪啪视频观看| 亚洲成人精品中文字幕电影| 日日干狠狠操夜夜爽| 精品99又大又爽又粗少妇毛片| 久久久久久久久久人人人人人人| 亚洲av.av天堂| 一区二区三区免费毛片| 观看免费一级毛片| 午夜福利高清视频| 一级毛片电影观看| 亚洲精品第二区| 亚洲在久久综合| 少妇被粗大猛烈的视频| 超碰av人人做人人爽久久| 啦啦啦啦在线视频资源| av在线亚洲专区| 亚洲av成人精品一二三区| 欧美xxxx黑人xx丫x性爽| 丝袜喷水一区| 亚洲精品aⅴ在线观看| or卡值多少钱| 久久草成人影院| 一个人看的www免费观看视频| 七月丁香在线播放| 国产淫片久久久久久久久| 午夜免费观看性视频| 精品人妻偷拍中文字幕| av一本久久久久| 国产黄色视频一区二区在线观看| 精品久久久久久久人妻蜜臀av| 国产亚洲5aaaaa淫片| 免费观看性生交大片5| 久久久久久久久久久丰满| 国产一级毛片七仙女欲春2| www.色视频.com| 欧美性猛交╳xxx乱大交人| 久久精品久久久久久久性| 国产高清三级在线| 肉色欧美久久久久久久蜜桃 | 亚洲av成人精品一二三区| 最后的刺客免费高清国语| 天堂av国产一区二区熟女人妻| 精华霜和精华液先用哪个| 一区二区三区乱码不卡18| 午夜精品一区二区三区免费看| 欧美精品一区二区大全| 久久久精品欧美日韩精品| 蜜桃亚洲精品一区二区三区| .国产精品久久| 一边亲一边摸免费视频| 国产乱来视频区| 亚洲成色77777| 国国产精品蜜臀av免费| 亚洲精品国产av蜜桃| 99视频精品全部免费 在线| 免费不卡的大黄色大毛片视频在线观看 | 欧美3d第一页| 最近中文字幕2019免费版| 色哟哟·www| 亚洲精品456在线播放app| 欧美日韩综合久久久久久| 亚洲自偷自拍三级| 午夜激情福利司机影院| 岛国毛片在线播放| 一级毛片 在线播放| 99久国产av精品国产电影| 秋霞伦理黄片| 国产 一区精品| 69av精品久久久久久| 一级爰片在线观看| 国产高潮美女av| 国产高清有码在线观看视频| 亚洲美女搞黄在线观看| 1000部很黄的大片| 观看美女的网站| 久久精品国产鲁丝片午夜精品| 国产午夜福利久久久久久| 午夜激情久久久久久久| 麻豆av噜噜一区二区三区| 网址你懂的国产日韩在线| 精品久久久久久久久av| 极品教师在线视频| 欧美三级亚洲精品| 乱人视频在线观看| 国产综合懂色| 只有这里有精品99| 能在线免费看毛片的网站| 成人欧美大片| 国产不卡一卡二| 午夜精品在线福利| 国产伦一二天堂av在线观看| 国产 亚洲一区二区三区 | 一级毛片我不卡| 少妇的逼水好多| 又大又黄又爽视频免费| 久久久久九九精品影院| 亚洲丝袜综合中文字幕| 五月伊人婷婷丁香| 久久国内精品自在自线图片| av一本久久久久| 伊人久久国产一区二区| 婷婷色麻豆天堂久久| videossex国产| 欧美高清成人免费视频www| 亚洲成人中文字幕在线播放| 国产亚洲av嫩草精品影院| 免费av观看视频| 美女国产视频在线观看| 秋霞在线观看毛片| 日本免费在线观看一区| 校园人妻丝袜中文字幕| 国产色爽女视频免费观看| 91aial.com中文字幕在线观看| 久久午夜福利片| 欧美高清成人免费视频www| 久久久成人免费电影| 久久久久性生活片| 一个人看的www免费观看视频| 麻豆国产97在线/欧美| 精品人妻一区二区三区麻豆| 汤姆久久久久久久影院中文字幕 | 亚洲美女搞黄在线观看| 日韩欧美 国产精品| 少妇猛男粗大的猛烈进出视频 | 国产伦精品一区二区三区四那| 日韩中字成人| 能在线免费观看的黄片| www.av在线官网国产| 国产av不卡久久| 欧美成人精品欧美一级黄| 国产免费一级a男人的天堂| 久久6这里有精品| 欧美变态另类bdsm刘玥| 成人国产麻豆网| 最近的中文字幕免费完整| 亚洲成色77777| 日韩强制内射视频| 精品久久久久久久久久久久久| 自拍偷自拍亚洲精品老妇| 国产亚洲精品久久久com| 成人高潮视频无遮挡免费网站| 精品久久国产蜜桃| 又爽又黄a免费视频| 欧美日韩在线观看h| 熟女人妻精品中文字幕| 女人久久www免费人成看片| 亚洲av成人av| 午夜福利成人在线免费观看| 日韩精品青青久久久久久| 毛片女人毛片| 黄色一级大片看看| 青春草亚洲视频在线观看| 日韩视频在线欧美| 草草在线视频免费看| 秋霞在线观看毛片| 日本-黄色视频高清免费观看| 69av精品久久久久久| 中国美白少妇内射xxxbb| 亚洲欧洲日产国产| 国产 一区 欧美 日韩| 26uuu在线亚洲综合色| 麻豆精品久久久久久蜜桃| 精品熟女少妇av免费看| 欧美zozozo另类| 日日摸夜夜添夜夜添av毛片| 亚洲一级一片aⅴ在线观看| 国产v大片淫在线免费观看| 亚洲精品乱码久久久久久按摩| 丰满少妇做爰视频| www.av在线官网国产| 男女边吃奶边做爰视频| 国产综合懂色| 亚洲av成人精品一二三区| 亚洲精品456在线播放app| 亚洲av国产av综合av卡| 亚洲国产最新在线播放| 99re6热这里在线精品视频| 舔av片在线| 久久这里有精品视频免费| 亚洲av免费在线观看| 高清毛片免费看| 精华霜和精华液先用哪个| 一级片'在线观看视频| 精品久久久精品久久久| 国内精品一区二区在线观看| 国产高清国产精品国产三级 | 九色成人免费人妻av| 精品人妻一区二区三区麻豆| 日本色播在线视频| 亚洲成人久久爱视频| 欧美三级亚洲精品| av国产久精品久网站免费入址| 国产成人freesex在线| 舔av片在线| 欧美一级a爱片免费观看看| 亚洲精品成人久久久久久| 成年免费大片在线观看| 91午夜精品亚洲一区二区三区| 国产精品嫩草影院av在线观看| 老司机影院毛片| 亚洲欧美日韩卡通动漫| 99re6热这里在线精品视频| 青春草亚洲视频在线观看| 国产在线男女| 有码 亚洲区| 欧美潮喷喷水| 亚洲精品乱久久久久久| 高清在线视频一区二区三区| 男女下面进入的视频免费午夜| 天堂俺去俺来也www色官网 | 国产成人福利小说| 中文在线观看免费www的网站| 色综合亚洲欧美另类图片| 午夜精品国产一区二区电影 | 国产高潮美女av| 国产精品久久久久久精品电影小说 | 熟妇人妻不卡中文字幕| 免费播放大片免费观看视频在线观看| 亚洲精品乱码久久久v下载方式| 色综合亚洲欧美另类图片| 国产 一区 欧美 日韩| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品美女特级片免费视频播放器| 亚洲精品aⅴ在线观看| 精品欧美国产一区二区三| 女人久久www免费人成看片| 日韩av不卡免费在线播放| 久久久色成人| 亚洲电影在线观看av| 日韩欧美精品免费久久| 最近中文字幕2019免费版| 一级毛片我不卡| 床上黄色一级片| 亚洲av成人精品一区久久| 亚州av有码| 免费少妇av软件| 精品久久久久久久久av| 成人一区二区视频在线观看| 一级a做视频免费观看| 亚洲人成网站在线观看播放| 亚洲国产欧美人成| 看黄色毛片网站| 18禁在线播放成人免费| 国产精品99久久久久久久久| 天美传媒精品一区二区| 搞女人的毛片| 69人妻影院| 人妻制服诱惑在线中文字幕| 亚洲aⅴ乱码一区二区在线播放| 少妇高潮的动态图| 精品久久久久久久末码| 九九在线视频观看精品| 国产在视频线精品| 亚洲欧美精品自产自拍| 国产精品熟女久久久久浪| 亚洲av日韩在线播放| 国产精品一区www在线观看| 久久精品人妻少妇| 麻豆精品久久久久久蜜桃| 自拍偷自拍亚洲精品老妇| 亚洲精品第二区| 亚洲自拍偷在线| 国产精品蜜桃在线观看| 亚洲欧美成人综合另类久久久| 日本三级黄在线观看| 一级毛片 在线播放| 精品久久久噜噜| 有码 亚洲区| 国产精品不卡视频一区二区| 嫩草影院精品99| 亚洲在久久综合| 国产成人a∨麻豆精品| 最近中文字幕2019免费版| 男插女下体视频免费在线播放| 天堂中文最新版在线下载 | 亚洲乱码一区二区免费版| 久久久久久久大尺度免费视频| 成人国产麻豆网| 中文在线观看免费www的网站| 欧美精品国产亚洲| 国产又色又爽无遮挡免| 综合色丁香网| 亚洲成色77777| 高清欧美精品videossex| av在线蜜桃| 非洲黑人性xxxx精品又粗又长| 日韩欧美精品v在线| 国产久久久一区二区三区| 丰满人妻一区二区三区视频av| 国产成人a区在线观看| 少妇的逼好多水| 欧美变态另类bdsm刘玥| av福利片在线观看| 午夜福利视频1000在线观看| 男女啪啪激烈高潮av片| 狠狠精品人妻久久久久久综合| 亚洲av福利一区| av在线亚洲专区| 能在线免费看毛片的网站| 中文字幕av成人在线电影| 18禁动态无遮挡网站| 国产成人精品婷婷| 色5月婷婷丁香| videossex国产| 边亲边吃奶的免费视频| a级毛色黄片| 亚洲国产精品专区欧美| 赤兔流量卡办理| 69av精品久久久久久| 我的老师免费观看完整版| 婷婷色av中文字幕| 日韩欧美三级三区| 97超碰精品成人国产| 一区二区三区乱码不卡18| 亚洲av中文av极速乱| 99热这里只有精品一区| 国产黄a三级三级三级人| 日韩,欧美,国产一区二区三区| 精品欧美国产一区二区三| 国产一区二区三区av在线|