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

    錯(cuò)列雙圓柱渦激振動(dòng)干涉數(shù)值模擬*

    2022-10-13 11:45:46王曉凱趙德光
    石油機(jī)械 2022年9期
    關(guān)鍵詞:渦激來流振幅

    張 睿 婁 敏 王曉凱 趙德光

    (1.中國石油大學(xué)(華東)石油工程學(xué)院 2.海洋石油工程(青島)有限公司)

    0 引 言

    隨著海洋油氣的開發(fā),海洋立管作為海上油氣開發(fā)中最為關(guān)鍵和易損的部件,其渦激振動(dòng)問題受到廣泛關(guān)注[1]。在實(shí)際工程中,一座海洋平臺(tái)通常會(huì)連接多根海洋立管,立管與立管會(huì)形成多圓柱系統(tǒng)。不同于單根圓柱的渦激振動(dòng),由于存在立管間尾流干涉的影響,多圓柱系統(tǒng)產(chǎn)生的渦激振動(dòng)特性更為復(fù)雜[2-3]。

    杜曉慶等[4]通過數(shù)值模擬對(duì)串列圓柱中上游圓柱的運(yùn)動(dòng)狀態(tài)對(duì)下游圓柱渦激振動(dòng)的影響進(jìn)行了研究,發(fā)現(xiàn)上游圓柱的振動(dòng)會(huì)提高下游圓柱的起振速度和橫流向最大振幅,并且兩類串列圓柱在振動(dòng)過程中的間距比和下游圓柱發(fā)生渦激振動(dòng)的機(jī)制都有所不同。王曉凱等[5]在尾流干涉下考慮流固耦合作用對(duì)小尺寸串聯(lián)圓柱進(jìn)行仿真研究,發(fā)現(xiàn)串聯(lián)圓柱的尾渦脫落模式隨約化速度增大會(huì)發(fā)生顯著變化,進(jìn)而影響圓柱的運(yùn)動(dòng)響應(yīng)和受力情況。林凌霄等[6]基于CFD方法,對(duì)不同間距比下的并列雙圓柱繞流進(jìn)行了數(shù)值模擬,模擬結(jié)果表明,隨著間距比的改變,兩圓柱尾流會(huì)呈現(xiàn)不同的形態(tài)。F.J.HUERA-HUARTE等[7]通過試驗(yàn)研究了高雷諾數(shù)下柔性并聯(lián)雙圓柱的流致振動(dòng),發(fā)現(xiàn)當(dāng)間隙比大于2.5時(shí),兩圓柱之間的干擾耦合影響較小。葉澤華等[8]在改進(jìn)CBS(Characteristic Based Split)有限元求解流場的基礎(chǔ)上,通過改進(jìn)的流-固耦合有限元解法計(jì)算渦激振動(dòng)問題,模擬了低雷諾數(shù)下兩錯(cuò)列圓柱渦激振動(dòng)過程,發(fā)現(xiàn)來流角度改變對(duì)于下游圓柱的流場、受力及振動(dòng)特性有著顯著影響。徐萬海等[9]采用模態(tài)分析法將試驗(yàn)測量應(yīng)變轉(zhuǎn)化為位移,開展了3種排布方式下雙柔性圓柱流激振動(dòng)的模型試驗(yàn),試驗(yàn)結(jié)果表明,當(dāng)柔性圓柱交錯(cuò)排布時(shí),下游圓柱橫流向和順流向的控制頻率接近,順流向的位移顯著增大。

    目前,對(duì)雙圓柱系統(tǒng)的研究主要集中在固定圓柱繞流特性及串、并聯(lián)基礎(chǔ)布置的渦激振動(dòng)研究上。在工程實(shí)際中,兩圓柱錯(cuò)列排布相比串聯(lián)或并聯(lián)情況更為常見,然而關(guān)于錯(cuò)列雙圓柱渦激振動(dòng)特性的研究還相對(duì)較少[10]。為進(jìn)一步探究影響錯(cuò)列雙圓柱渦激振動(dòng)特性的敏感間距及排列角度,筆者基于重疊網(wǎng)格(Overset Mesh)技術(shù)進(jìn)行網(wǎng)格劃分,將振動(dòng)模型簡化為雙自由度系統(tǒng),建立二維流固耦合數(shù)值模型,使用四階龍格庫塔法求解振動(dòng)方程,編譯用戶自定義函數(shù)(UDF)并結(jié)合Fluent動(dòng)網(wǎng)格技術(shù)進(jìn)行求解,在對(duì)所建模型進(jìn)行驗(yàn)證的基礎(chǔ)上,對(duì)錯(cuò)列雙圓柱開展渦激振動(dòng)響應(yīng)研究,研究間距比、來流角度和約化速度對(duì)錯(cuò)列雙圓柱渦激振動(dòng)特性的影響。所得結(jié)果可為實(shí)際海洋立管的布置提供參考。

    1 數(shù)值方法

    1.1 控制方程

    連續(xù)性方程和動(dòng)量方程:

    ?u=0

    (1)

    (2)

    式中:u表示流體的速度,m/s;p表示壓力,Pa;Re表示雷諾數(shù);t表示時(shí)間,s。

    將圓柱振動(dòng)模型簡化為二維彈簧-質(zhì)量-阻尼系統(tǒng),振動(dòng)模型控制方程如下:

    (3)

    (4)

    式中:x、y分別為圓柱順流向和橫流向位移,m;m為單位長度質(zhì)量,kg/m;c為阻尼系數(shù),N/(m/s);k為剛度系數(shù),N/m;FD(t)、FL(t)分別為順流向和橫橫向的流體力,N。

    1.2 流固耦合實(shí)現(xiàn)

    針對(duì)立管的渦激振動(dòng)問題,本文采用動(dòng)網(wǎng)格技術(shù)實(shí)現(xiàn)流固耦合,具體如下:通過Fluent對(duì)圓柱的運(yùn)動(dòng)和圓柱周圍流場變化進(jìn)行模擬,提取圓柱受到的升力和阻力系數(shù);在UDF中采用四階龍格庫塔法對(duì)二維振動(dòng)控制方程進(jìn)行離散并求解,得到每個(gè)時(shí)間步內(nèi)對(duì)應(yīng)的立管速度和位移。

    1.3 模型參數(shù)

    本文模擬所取圓柱直徑D為0.02 m,質(zhì)量比m*為2.55,阻尼比ξ為0.003,固有頻率為1.72 Hz。兩圓柱錯(cuò)列排布,如圖1所示。圓柱間距S為4D,來流角度α=45°,兩圓柱圓心到流場中軸線的距離相同,流場域左側(cè)距離上游圓柱圓心10D,流場域右側(cè)距離上游圓柱圓心30D。

    圖1 圓柱排布方式圖Fig.1 Arrangement of cylinders

    通過Overset Mesh技術(shù)進(jìn)行流場模型的網(wǎng)格劃分[5],以保證較好的網(wǎng)格質(zhì)量,網(wǎng)格劃分結(jié)果如圖2所示。網(wǎng)格類型均采用結(jié)構(gòu)化四邊形網(wǎng)格,對(duì)尾流區(qū)域的網(wǎng)格進(jìn)行加密處理,其他區(qū)域逐漸增大網(wǎng)格密度,以此合理控制網(wǎng)格總量,以便節(jié)約計(jì)算時(shí)間。

    圖2 網(wǎng)格劃分示意圖Fig.2 Schematic diagram of mesh division

    計(jì)算流場域入口采用均勻速度邊界,出口選擇壓力邊界,上下兩側(cè)面邊界為對(duì)稱邊界條件,圓柱表面采用無滑移邊界條件。通過湍流強(qiáng)度和水力直徑定義湍流,計(jì)算可得湍流強(qiáng)度和水力直徑;湍流模型采用SSTk-ω模型,對(duì)流項(xiàng)采用二階迎風(fēng)離散格式,時(shí)間項(xiàng)采用全隱式積分;通過Coupled算法實(shí)現(xiàn)控制方程中速度與壓力的耦合。根據(jù)網(wǎng)格相對(duì)運(yùn)動(dòng)不超過劃分網(wǎng)格的最小尺寸的要求,選取迭代時(shí)間步長為0.01 s。

    2 模型驗(yàn)證

    2.1 圓柱繞流驗(yàn)證

    對(duì)雷諾數(shù)Re=200、間距比S/D=3時(shí)的串聯(lián)雙圓柱體繞流進(jìn)行數(shù)值模擬,將所得結(jié)果與文獻(xiàn)[11-13]的模擬結(jié)果進(jìn)行對(duì)比,結(jié)果如表1所示。從表1可以看出,本文圓柱繞流模擬所得上游阻力系數(shù)Cd1、下游阻力系數(shù)Cd2和斯特勞哈爾數(shù)St與文獻(xiàn)中數(shù)值結(jié)果相近,由于網(wǎng)格、湍流模型以及數(shù)值處理的方法不同,存在一定誤差。這說明本文所采用的網(wǎng)格劃分方法、計(jì)算模型設(shè)置的參數(shù)以及模擬結(jié)果正確。

    表1 雙圓柱繞流模擬結(jié)果對(duì)比Table 1 Comparison of simulation results of flow around double cylinders

    2.2 單圓柱渦激振動(dòng)模擬驗(yàn)證

    由于圓柱流固耦合研究涉及到流場網(wǎng)格的動(dòng)態(tài)變化,與圓柱繞流相比,對(duì)網(wǎng)格質(zhì)量、計(jì)算參數(shù)選取及自定義函數(shù)的準(zhǔn)確性等有著嚴(yán)格的要求,所以需要對(duì)圓柱模型進(jìn)一步開展模擬驗(yàn)證。

    計(jì)算模型采用前文的網(wǎng)格劃分方法及數(shù)值求解格式要求,嚴(yán)格按照文獻(xiàn)[14-15]中的參數(shù)進(jìn)行數(shù)值模擬,將本文模擬結(jié)果與文獻(xiàn)結(jié)果進(jìn)行對(duì)比,結(jié)果如圖3所示。

    圖3 橫流向振幅隨約化速度的變化關(guān)系Fig.3 Variation of transverse amplitude with reduced velocity

    采用振幅均方根與圓柱直徑之比(y/D)對(duì)橫流向振動(dòng)響應(yīng)進(jìn)行描述。從圖3可以看出,本文模擬結(jié)果與參考文獻(xiàn)結(jié)果相似,確定了自定義函數(shù)UDF的準(zhǔn)確性以及流固耦合計(jì)算模型的可靠性,為后文錯(cuò)列圓柱渦激振動(dòng)數(shù)值模擬研究奠定了基礎(chǔ)。

    3 計(jì)算結(jié)果

    基于上述數(shù)值模型,選擇無量綱量約化速度Ur=1.0~15.0,即雷諾數(shù)Re=680~10 320。首先對(duì)間距比S/D=3.0、4.0及5.0的錯(cuò)列排布雙圓柱立管渦激振動(dòng)特性開展研究;同時(shí)研究來流角度α=15°、45°及75°工況下錯(cuò)列雙圓柱渦激振動(dòng)特性。

    3.1 錯(cuò)列圓柱間距比研究

    3.1.1 圓柱振動(dòng)頻率與幅值

    保持來流角度α=45°不變,在不同間距比下,錯(cuò)列圓柱的漩渦脫落頻率fs與固有頻率fn的比值λ隨約化速度的變化如圖4所示。由圖4可知,在低約化速度范圍內(nèi),不同間距比的各條曲線走向幾乎一致,而在進(jìn)入高約化速度范圍Ur>10后,錯(cuò)列兩圓柱的頻率比在數(shù)值上存在較大差距,上游圓柱明顯大于下游圓柱。

    圖4 不同間距比下頻率比隨Ur的變化關(guān)系Fig.4 Variation of frequency ratio with reduced velocity at different spacing ratios

    不同間距比下錯(cuò)列圓柱的橫流向振動(dòng)響應(yīng)隨約化速度的變化關(guān)系如圖5a和圖5b所示,順流向振動(dòng)響應(yīng)隨約化速度的變化關(guān)系如圖5c和圖5d所示。

    圖5 圓柱渦激振動(dòng)響應(yīng)隨Ur的變化關(guān)系Fig.5 Variation of vortex-induced vibration response of cylinder with reduced velocity

    從圖5可以看出,當(dāng)約化速度Ur=0~4時(shí)圓柱振動(dòng)處于上分支,Ur=4~10為鎖振區(qū)間,Ur=10~15時(shí)處于下分支。但是間距比S/D=3時(shí)下游圓柱的曲線格外突出,在兩個(gè)流向上的振幅會(huì)明顯大于其他圓柱。這是因?yàn)殚g距較近,間隙流相當(dāng)于增大了來流流速,因此下游圓柱位移幅值明顯增大。隨著間距比的增大,圓柱間的干涉作用減弱,上、下游圓柱的振幅曲線不斷向孤立圓柱的振幅曲線靠近,橫流向振幅變化呈現(xiàn)三段式變化趨勢,順流向平均振幅(x/D)隨著來流速度的增大持續(xù)增大。

    3.1.2 運(yùn)動(dòng)軌跡分析

    圓柱運(yùn)動(dòng)軌跡如圖6所示。由圖6可知,當(dāng)間距比S/D=3時(shí),兩圓柱的運(yùn)動(dòng)軌跡較為雜亂,無明顯“8”字形,下游圓柱受上游圓柱的流動(dòng)干涉影響出現(xiàn)了許多不穩(wěn)定的非周期性運(yùn)動(dòng)軌跡,表明在S/D=3時(shí)錯(cuò)列兩圓柱的流動(dòng)干涉作用較為強(qiáng)烈。

    圖6 兩圓柱的運(yùn)動(dòng)軌跡(來流角度不變)Fig.6 Kinematic trajectories of double cylinders (at constant flow angle)

    隨著間距比的不斷增大,“8”字形、“桃”形等運(yùn)動(dòng)軌跡開始占據(jù)主導(dǎo),圓柱運(yùn)動(dòng)的規(guī)律性得到提升,運(yùn)動(dòng)軌跡由繁到簡,兩圓柱軌跡間的差異越來越小,上游圓柱和下游圓柱的運(yùn)動(dòng)軌跡隨著S/D的增大慢慢趨于一致。S/D=4和5時(shí),上游圓柱對(duì)下游圓柱的干涉效應(yīng)顯著高于下游圓柱對(duì)上游圓柱的干涉反饋效應(yīng)。

    3.1.3 流場分析

    圖7反映了不同間距比下錯(cuò)列兩圓柱的尾渦模態(tài)隨約化速度的變化關(guān)系。

    圖7 不同間距比下尾渦模態(tài)隨Ur的變化Fig.7 Variation of trailing vortex mode with reduced velocity at different spacing ratios

    由圖7可知:在間距比S/D=3時(shí),隨著約化速度增大,流體與圓柱之間的耦合作用強(qiáng)烈,不能在圓柱后方形成規(guī)則脫落的漩渦;當(dāng)約化速度大于10后,兩圓柱脫離“鎖振”狀態(tài),振幅減小。S/D=3與其他間距比相比,下游圓柱受到上游圓柱尾渦脫落的影響更大,在橫流向產(chǎn)生比其他2種情況更大的振幅;下游圓柱也對(duì)上游圓柱產(chǎn)生“抑制”,使上游圓柱的尾渦變得細(xì)長。在高約化速度下,隨著間距比的增大,下游圓柱對(duì)上游圓柱的“抑制”作用減弱,上游圓柱出現(xiàn)渦旋脫落,渦街形態(tài)逐漸變得規(guī)則,直到S/D=5時(shí),上、下游圓柱的頻率比曲線逐漸接近孤立圓柱。

    3.2 錯(cuò)列圓柱來流角度研究

    3.2.1 圓柱振動(dòng)頻率與幅值

    保持間距比S/D=4不變,研究來流角度和約化速度對(duì)流動(dòng)干涉作用的影響。角度越小錯(cuò)列圓柱就會(huì)越接近串聯(lián)圓柱,反之,角度越大兩圓柱的相對(duì)位置將趨于并聯(lián)排列。因此當(dāng)來流角度從小到大變化時(shí),即可視為從串聯(lián)圓柱向并聯(lián)圓柱的逐漸過渡。圖8為不同來流角度下錯(cuò)列兩圓柱的頻率比λ隨約化速度的變化關(guān)系。觀察圖8發(fā)現(xiàn),隨著來流角度α的增大,圓柱的“鎖振”開始位置出現(xiàn)前移傾向,且45°鎖振區(qū)間相比孤立圓柱及其他雙圓柱較為狹窄;在同一來流角度α?xí)r,上游圓柱與下游圓柱的“鎖振”區(qū)間長度接近。

    圖8 不同來流角度下頻率比隨Ur的變化關(guān)系Fig.8 Variation of frequency ratio with reduced velocity at different flow angles

    圖9為不同來流角度α下錯(cuò)列圓柱的渦激振動(dòng)響應(yīng)隨約化速度的變化關(guān)系。

    由圖9可以看出,α=15°時(shí)錯(cuò)列兩圓柱的幅值顯著大于孤立圓柱和其他來流角度時(shí)的數(shù)值。此外下游圓柱的振幅最終保持在0.28D附近,遠(yuǎn)遠(yuǎn)大于上游圓柱。不同來流角度下錯(cuò)列兩圓柱的順流向振幅隨約化速度的增加持續(xù)增大,各條振幅曲線相近,且與孤立圓柱差別不大,因此判斷錯(cuò)列圓柱來流角的改變對(duì)順流向平均振幅的影響比橫流向振幅的影響小。

    圖9 圓柱渦激振動(dòng)響應(yīng)隨Ur的變化關(guān)系Fig.9 Variation of vortex-induced vibration response of cylinder with reduced velocity

    3.2.2 圓柱運(yùn)動(dòng)軌跡

    圓柱運(yùn)動(dòng)的軌跡如圖10所示。由圖10可知:當(dāng)來流角度α=15°時(shí),錯(cuò)列兩圓柱的運(yùn)動(dòng)軌跡差異顯著,上游圓柱的運(yùn)動(dòng)軌跡構(gòu)成較為簡單,下游圓柱則以非周期性運(yùn)動(dòng)為主,軌跡較為雜亂且相互重合,說明此時(shí)下游圓柱所受的干涉影響較大;當(dāng)來流角度α=45°時(shí),圓柱運(yùn)動(dòng)的周期性有所改善,“8”字形軌跡開始出現(xiàn),下游圓柱的運(yùn)動(dòng)軌跡逐漸清晰簡單,原因是來流角度增大后受到上游圓柱流動(dòng)干涉作用的影響明顯降低;當(dāng)來流角度α=75°時(shí),錯(cuò)列兩圓柱運(yùn)動(dòng)軌跡的周期性進(jìn)一步提升,“8”字形軌跡占據(jù)主導(dǎo),上游、下游圓柱的運(yùn)動(dòng)軌跡愈發(fā)接近,并顯現(xiàn)出并聯(lián)兩圓柱的軌跡對(duì)稱特性。

    圖10 兩圓柱的運(yùn)動(dòng)軌跡(間距比不變)Fig.10 Kinematic trajectories of double cylinders (at constant spacing ratio)

    3.2.3 流場分析

    圖11為不同來流角度下錯(cuò)列兩圓柱的尾渦模態(tài)隨約化速度的變化關(guān)系。由圖11可知:當(dāng)來流角度α=15°時(shí),下游圓柱處于上游圓柱的尾流區(qū)中,兩圓柱漩渦產(chǎn)生干擾,漩渦結(jié)構(gòu)顯得比較復(fù)雜;隨著來流角度變大,下游圓柱逐漸脫離上游圓柱的尾流區(qū),兩圓柱的尾渦模態(tài)也變得簡單、規(guī)則,從單一渦街模式變化為雙渦街模式;當(dāng)進(jìn)入高約化速度且來流角度α=15°時(shí),上游圓柱的兩個(gè)剪切層并沒有形成漩渦脫落,但由于上游圓柱對(duì)下游圓柱的干涉作用,推遲了下游圓柱的游渦脫落,這也與下游圓柱橫向振幅受到限制逐漸減小和上、下游圓柱頻率比逐漸接近的現(xiàn)象相對(duì)應(yīng);而隨著來流角度增大,上、下游圓柱脫落的漩渦不再相互干擾,各自后方出現(xiàn)規(guī)則的渦街,呈現(xiàn)出與并聯(lián)圓柱接近的“同向同步”模態(tài)。

    圖11 不同來流角度下尾渦模態(tài)隨Ur的變化Fig.11 Variation of trailing vortex mode with reduced velocity at different flow angles

    4 結(jié) 論

    (1)當(dāng)間距比和來流角度較小時(shí),錯(cuò)列兩圓柱的流動(dòng)干涉作用較為強(qiáng)烈,兩圓柱的運(yùn)動(dòng)軌跡較為雜亂;隨著間距比和來流角度的增大,圓柱的運(yùn)動(dòng)軌跡逐漸變得規(guī)則,兩圓柱軌跡間的差異越來越小,上游圓柱和下游圓柱的運(yùn)動(dòng)軌跡將慢慢趨于一致。

    (2)在間距比S/D=3時(shí),流體與圓柱之間的耦合作用非常強(qiáng)烈,不能在圓柱后方形成規(guī)則脫落的漩渦;在高約化速度下,與其他間距比相比,下游圓柱受到上游圓柱尾渦脫落的影響更大,在橫流向產(chǎn)生更大的振幅;下游圓柱也對(duì)上游圓柱產(chǎn)生“抑制”,使上游圓柱的尾渦變得細(xì)長。

    (3)當(dāng)來流角度α=15°時(shí),上游圓柱的2個(gè)剪切層并沒有形成漩渦脫落,由于上游圓柱對(duì)下游圓柱的干涉作用,推遲了下游圓柱的渦脫落,這也與下游圓柱橫向振幅受到限制逐漸減小和上、下游圓柱頻率比逐漸接近的現(xiàn)象相對(duì)應(yīng)。

    猜你喜歡
    渦激來流振幅
    不同間距比下串聯(lián)圓柱渦激振動(dòng)數(shù)值模擬研究
    兩種典型來流條件下風(fēng)力機(jī)尾跡特性的數(shù)值研究
    能源工程(2022年2期)2022-05-23 13:51:48
    渦激振動(dòng)發(fā)電裝置及其關(guān)鍵技術(shù)
    不同來流條件對(duì)溢洪道過流能力的影響
    盤球立管結(jié)構(gòu)抑制渦激振動(dòng)的數(shù)值分析方法研究
    電子制作(2018年14期)2018-08-21 01:38:42
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    柔性圓管在渦激振動(dòng)下的模態(tài)響應(yīng)分析
    身体一侧抽搐| 又大又爽又粗| 妹子高潮喷水视频| 国语自产精品视频在线第100页| 国产亚洲精品综合一区在线观看 | 精品国产亚洲在线| 在线十欧美十亚洲十日本专区| 看免费av毛片| 午夜福利成人在线免费观看| 午夜老司机福利片| 麻豆av在线久日| 亚洲一码二码三码区别大吗| 久久精品国产亚洲av高清一级| 亚洲精品粉嫩美女一区| 色综合亚洲欧美另类图片| 国产主播在线观看一区二区| 亚洲欧美日韩无卡精品| 日韩欧美国产一区二区入口| 色哟哟哟哟哟哟| 国产精品亚洲一级av第二区| 国产亚洲av高清不卡| 欧美绝顶高潮抽搐喷水| 精品午夜福利视频在线观看一区| 97超级碰碰碰精品色视频在线观看| 真人一进一出gif抽搐免费| 亚洲人成电影观看| 麻豆av在线久日| 中亚洲国语对白在线视频| 亚洲精品中文字幕一二三四区| 国产精品精品国产色婷婷| 在线天堂中文资源库| 亚洲成人国产一区在线观看| 国产xxxxx性猛交| 日韩成人在线观看一区二区三区| 亚洲精品在线观看二区| 一边摸一边做爽爽视频免费| 亚洲色图av天堂| 亚洲一区二区三区不卡视频| 国产精品电影一区二区三区| 精品第一国产精品| 国产精品永久免费网站| 欧美成人免费av一区二区三区| 久9热在线精品视频| 两性夫妻黄色片| 青草久久国产| 色精品久久人妻99蜜桃| 国产精品 国内视频| 俄罗斯特黄特色一大片| 国产一区二区三区在线臀色熟女| 国产精品免费一区二区三区在线| 欧美午夜高清在线| 国产高清激情床上av| 美女午夜性视频免费| 村上凉子中文字幕在线| 亚洲全国av大片| 成人av一区二区三区在线看| 亚洲 欧美 日韩 在线 免费| 国产精品免费视频内射| 国语自产精品视频在线第100页| 久久国产精品影院| 757午夜福利合集在线观看| a级毛片在线看网站| 老汉色av国产亚洲站长工具| 午夜免费观看网址| 一级a爱片免费观看的视频| 伊人久久大香线蕉亚洲五| 国产高清videossex| 正在播放国产对白刺激| 亚洲第一av免费看| 国产亚洲av高清不卡| 动漫黄色视频在线观看| 国产精品电影一区二区三区| 久久久久精品国产欧美久久久| 欧美大码av| 啦啦啦韩国在线观看视频| 国产成人av激情在线播放| 中文字幕久久专区| 久久久久久久久中文| 亚洲午夜理论影院| 热re99久久国产66热| 色综合婷婷激情| 一级毛片高清免费大全| 午夜成年电影在线免费观看| 午夜福利成人在线免费观看| 午夜福利成人在线免费观看| 午夜成年电影在线免费观看| 正在播放国产对白刺激| 亚洲情色 制服丝袜| 亚洲精品av麻豆狂野| 看黄色毛片网站| 亚洲va日本ⅴa欧美va伊人久久| 一本久久中文字幕| 1024视频免费在线观看| 两性夫妻黄色片| 一进一出抽搐动态| 国产区一区二久久| 国产激情久久老熟女| 我的亚洲天堂| 国产1区2区3区精品| 成人欧美大片| 淫妇啪啪啪对白视频| 国产亚洲精品av在线| 757午夜福利合集在线观看| 国产精品亚洲美女久久久| 美女国产高潮福利片在线看| 精品国产超薄肉色丝袜足j| 亚洲精品国产区一区二| 99国产极品粉嫩在线观看| 亚洲伊人色综图| 色在线成人网| 国产一区在线观看成人免费| 淫秽高清视频在线观看| 91成人精品电影| 亚洲国产精品999在线| 在线观看免费视频网站a站| 黄片小视频在线播放| 亚洲男人天堂网一区| 一个人免费在线观看的高清视频| 欧美日韩乱码在线| 丝袜美腿诱惑在线| 99国产综合亚洲精品| av超薄肉色丝袜交足视频| 亚洲黑人精品在线| 午夜福利在线观看吧| 后天国语完整版免费观看| 日本vs欧美在线观看视频| 久久久国产成人精品二区| 色av中文字幕| 国产成人av激情在线播放| 亚洲中文日韩欧美视频| 午夜视频精品福利| 手机成人av网站| 国产亚洲精品久久久久久毛片| 999久久久国产精品视频| 身体一侧抽搐| 国产视频一区二区在线看| 日日爽夜夜爽网站| 啦啦啦观看免费观看视频高清 | 亚洲国产毛片av蜜桃av| 1024香蕉在线观看| 国产精品日韩av在线免费观看 | 人妻丰满熟妇av一区二区三区| 久久国产亚洲av麻豆专区| 色婷婷久久久亚洲欧美| 亚洲中文字幕日韩| 国产男靠女视频免费网站| 久久狼人影院| 免费在线观看完整版高清| 淫秽高清视频在线观看| 69av精品久久久久久| 欧美乱妇无乱码| 亚洲一区高清亚洲精品| 欧美日韩亚洲国产一区二区在线观看| 国产高清videossex| 国产99久久九九免费精品| 丝袜美足系列| 不卡一级毛片| 婷婷六月久久综合丁香| 少妇裸体淫交视频免费看高清 | 亚洲成a人片在线一区二区| 老司机福利观看| 夜夜躁狠狠躁天天躁| 麻豆一二三区av精品| 一区二区三区国产精品乱码| 成年版毛片免费区| 99久久久亚洲精品蜜臀av| 国产一区在线观看成人免费| 亚洲av成人av| 性色av乱码一区二区三区2| 1024视频免费在线观看| 久久久精品国产亚洲av高清涩受| 精品欧美一区二区三区在线| 看免费av毛片| 男男h啪啪无遮挡| 欧美最黄视频在线播放免费| 一边摸一边抽搐一进一出视频| 黄色 视频免费看| 国产真人三级小视频在线观看| 国产黄a三级三级三级人| 黄频高清免费视频| 中国美女看黄片| x7x7x7水蜜桃| 中文字幕最新亚洲高清| 精品久久久久久久人妻蜜臀av | 亚洲全国av大片| 91麻豆精品激情在线观看国产| 国产精品精品国产色婷婷| 91国产中文字幕| 99久久综合精品五月天人人| 国产1区2区3区精品| 一边摸一边抽搐一进一出视频| 男人舔女人的私密视频| 国产精品久久久久久精品电影 | 亚洲性夜色夜夜综合| 午夜免费激情av| 色综合欧美亚洲国产小说| 午夜久久久在线观看| 精品第一国产精品| 欧美激情久久久久久爽电影 | 日日干狠狠操夜夜爽| 久久性视频一级片| 欧美成人免费av一区二区三区| 可以在线观看毛片的网站| 岛国在线观看网站| 国产男靠女视频免费网站| 国产欧美日韩精品亚洲av| 高清黄色对白视频在线免费看| 好男人在线观看高清免费视频 | 国产精品久久电影中文字幕| 岛国在线观看网站| 午夜福利高清视频| 男女下面进入的视频免费午夜 | 国产成人欧美| 91成人精品电影| 岛国视频午夜一区免费看| 男女下面插进去视频免费观看| 免费在线观看黄色视频的| 国产成人啪精品午夜网站| 黑人操中国人逼视频| 日本黄色视频三级网站网址| 成人三级黄色视频| 久久久国产精品麻豆| 在线观看免费日韩欧美大片| 久久久久久久精品吃奶| 国产成人系列免费观看| 午夜免费观看网址| 中文字幕另类日韩欧美亚洲嫩草| 99香蕉大伊视频| 黑人欧美特级aaaaaa片| 不卡av一区二区三区| 97碰自拍视频| www.999成人在线观看| 国产精品 国内视频| 亚洲成国产人片在线观看| 国产激情欧美一区二区| 手机成人av网站| 搡老岳熟女国产| 国产一区二区三区在线臀色熟女| 一区二区三区精品91| 正在播放国产对白刺激| 99国产精品一区二区三区| 乱人伦中国视频| 欧美激情极品国产一区二区三区| 两性夫妻黄色片| 色尼玛亚洲综合影院| 一级黄色大片毛片| 美女免费视频网站| 国产精品精品国产色婷婷| 国产国语露脸激情在线看| 国产精品一区二区在线不卡| 曰老女人黄片| 久久人人精品亚洲av| 久久久国产精品麻豆| 一级作爱视频免费观看| 淫妇啪啪啪对白视频| 亚洲熟女毛片儿| 国产主播在线观看一区二区| 91av网站免费观看| 三级毛片av免费| 高潮久久久久久久久久久不卡| 黑人巨大精品欧美一区二区mp4| 午夜a级毛片| 国产成人精品无人区| 在线av久久热| 欧美乱色亚洲激情| 欧美+亚洲+日韩+国产| 长腿黑丝高跟| 成年人黄色毛片网站| 久久影院123| 一本综合久久免费| 999久久久精品免费观看国产| 日韩欧美一区视频在线观看| 国内毛片毛片毛片毛片毛片| 精品一区二区三区四区五区乱码| 脱女人内裤的视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲第一青青草原| 欧美在线一区亚洲| 国产亚洲欧美精品永久| 两性午夜刺激爽爽歪歪视频在线观看 | 他把我摸到了高潮在线观看| 久久精品亚洲熟妇少妇任你| 伊人久久大香线蕉亚洲五| 国产欧美日韩精品亚洲av| 99在线视频只有这里精品首页| 亚洲精品国产色婷婷电影| 777久久人妻少妇嫩草av网站| 欧美在线黄色| 狠狠狠狠99中文字幕| 精品欧美一区二区三区在线| 一本大道久久a久久精品| 亚洲全国av大片| 国产三级黄色录像| 99久久国产精品久久久| 亚洲国产日韩欧美精品在线观看 | tocl精华| 国产高清videossex| 777久久人妻少妇嫩草av网站| 免费少妇av软件| 亚洲av第一区精品v没综合| 操美女的视频在线观看| 亚洲午夜理论影院| 亚洲精品在线观看二区| 国产精品自产拍在线观看55亚洲| 一区在线观看完整版| 日韩大尺度精品在线看网址 | 久久久久久久久久久久大奶| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕精品免费在线观看视频| 国产一级毛片七仙女欲春2 | 午夜福利欧美成人| 少妇裸体淫交视频免费看高清 | 韩国av一区二区三区四区| 999久久久国产精品视频| av有码第一页| 后天国语完整版免费观看| 又黄又粗又硬又大视频| 国产亚洲精品第一综合不卡| 美女免费视频网站| 国产精品爽爽va在线观看网站 | 美女 人体艺术 gogo| 亚洲一区中文字幕在线| 在线观看免费视频网站a站| 久久精品91无色码中文字幕| 99久久精品国产亚洲精品| 热re99久久国产66热| 97碰自拍视频| 久久精品影院6| 禁无遮挡网站| 怎么达到女性高潮| 黄色 视频免费看| 日本精品一区二区三区蜜桃| 视频在线观看一区二区三区| 级片在线观看| 亚洲 欧美 日韩 在线 免费| 99在线人妻在线中文字幕| 波多野结衣一区麻豆| 黄色a级毛片大全视频| 国产高清视频在线播放一区| 亚洲av日韩精品久久久久久密| 无人区码免费观看不卡| 熟妇人妻久久中文字幕3abv| 亚洲片人在线观看| 亚洲精品久久国产高清桃花| 91老司机精品| 欧美一区二区精品小视频在线| 亚洲精品粉嫩美女一区| bbb黄色大片| 精品久久久精品久久久| 色在线成人网| 亚洲七黄色美女视频| 久久久久精品国产欧美久久久| 亚洲 欧美 日韩 在线 免费| 国产亚洲欧美在线一区二区| e午夜精品久久久久久久| 成人国产综合亚洲| 两个人免费观看高清视频| 久久久久国内视频| 啦啦啦韩国在线观看视频| 乱人伦中国视频| 制服丝袜大香蕉在线| 自拍欧美九色日韩亚洲蝌蚪91| 精品不卡国产一区二区三区| 国产精品二区激情视频| 99久久国产精品久久久| 亚洲自偷自拍图片 自拍| 99re在线观看精品视频| 视频在线观看一区二区三区| 亚洲精品av麻豆狂野| 亚洲自拍偷在线| 久久天躁狠狠躁夜夜2o2o| 一级毛片女人18水好多| 波多野结衣高清无吗| 满18在线观看网站| 欧美大码av| 纯流量卡能插随身wifi吗| or卡值多少钱| 日本vs欧美在线观看视频| 国产亚洲精品一区二区www| 久久久国产欧美日韩av| av超薄肉色丝袜交足视频| 日韩大尺度精品在线看网址 | 国产成人精品无人区| 正在播放国产对白刺激| 一二三四社区在线视频社区8| 亚洲国产欧美一区二区综合| 国产在线精品亚洲第一网站| 黑人巨大精品欧美一区二区蜜桃| 18禁美女被吸乳视频| 最新美女视频免费是黄的| 韩国精品一区二区三区| 亚洲av熟女| 亚洲 欧美一区二区三区| 无遮挡黄片免费观看| 男人的好看免费观看在线视频 | 亚洲熟妇中文字幕五十中出| 三级毛片av免费| 好看av亚洲va欧美ⅴa在| av视频免费观看在线观看| 精品无人区乱码1区二区| 一边摸一边抽搐一进一出视频| 亚洲,欧美精品.| 精品一区二区三区四区五区乱码| 美国免费a级毛片| 婷婷六月久久综合丁香| 国产极品粉嫩免费观看在线| 超碰成人久久| 国产精品一区二区三区四区久久 | 丝袜美腿诱惑在线| 欧美成人免费av一区二区三区| 变态另类成人亚洲欧美熟女 | 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 天堂√8在线中文| 久久精品国产99精品国产亚洲性色 | 色精品久久人妻99蜜桃| 中文字幕av电影在线播放| 久久国产精品男人的天堂亚洲| 久久香蕉激情| 免费无遮挡裸体视频| 亚洲精品国产精品久久久不卡| 久久伊人香网站| 免费在线观看影片大全网站| 51午夜福利影视在线观看| 狂野欧美激情性xxxx| 亚洲欧美日韩另类电影网站| 久久人妻福利社区极品人妻图片| www.精华液| 男女午夜视频在线观看| 一级片免费观看大全| 一级毛片精品| 午夜福利欧美成人| 99香蕉大伊视频| 午夜激情av网站| 国产成人精品在线电影| 亚洲精品国产色婷婷电影| 国产亚洲精品第一综合不卡| 麻豆av在线久日| 亚洲av片天天在线观看| 国产精品久久久人人做人人爽| 好男人在线观看高清免费视频 | 好看av亚洲va欧美ⅴa在| 婷婷六月久久综合丁香| 99久久精品国产亚洲精品| 亚洲熟女毛片儿| 亚洲av第一区精品v没综合| 啦啦啦观看免费观看视频高清 | 最新在线观看一区二区三区| 12—13女人毛片做爰片一| 日日干狠狠操夜夜爽| 精品国内亚洲2022精品成人| 久久久久亚洲av毛片大全| 18美女黄网站色大片免费观看| 国产成人精品久久二区二区免费| 日本欧美视频一区| 亚洲五月婷婷丁香| 中文字幕精品免费在线观看视频| 黑人欧美特级aaaaaa片| 99国产精品一区二区三区| 国产一区二区在线av高清观看| 精品久久久久久成人av| 免费看a级黄色片| 真人一进一出gif抽搐免费| 亚洲国产精品合色在线| 久久久久久大精品| 成人18禁高潮啪啪吃奶动态图| 欧美av亚洲av综合av国产av| 国产av又大| 777久久人妻少妇嫩草av网站| 国产精品av久久久久免费| 一级毛片女人18水好多| 国产精品久久久久久精品电影 | 欧美日韩福利视频一区二区| 在线观看66精品国产| 嫩草影院精品99| 黄色a级毛片大全视频| 日本 av在线| 性欧美人与动物交配| 日韩免费av在线播放| 视频在线观看一区二区三区| 国产亚洲精品一区二区www| 国产成年人精品一区二区| 国产亚洲欧美精品永久| 黑人巨大精品欧美一区二区蜜桃| 国产精品久久视频播放| 欧美 亚洲 国产 日韩一| 免费高清视频大片| 久久热在线av| 黄色女人牲交| 91麻豆av在线| 97碰自拍视频| 18美女黄网站色大片免费观看| 天天一区二区日本电影三级 | 一级毛片女人18水好多| 国产一区在线观看成人免费| 精品一区二区三区四区五区乱码| 亚洲国产高清在线一区二区三 | 国产一卡二卡三卡精品| 老司机午夜福利在线观看视频| 亚洲专区中文字幕在线| tocl精华| 黑人巨大精品欧美一区二区mp4| 国内毛片毛片毛片毛片毛片| 91麻豆av在线| 免费看a级黄色片| 99在线视频只有这里精品首页| 一边摸一边抽搐一进一小说| tocl精华| 精品福利观看| 男女之事视频高清在线观看| 亚洲av成人一区二区三| 国产三级黄色录像| 一级毛片女人18水好多| 岛国视频午夜一区免费看| 国产av一区在线观看免费| 日本 av在线| 精品久久蜜臀av无| 涩涩av久久男人的天堂| 日本vs欧美在线观看视频| 999久久久国产精品视频| 亚洲人成电影免费在线| 精品国产亚洲在线| av网站免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产亚洲精品综合一区在线观看 | 91大片在线观看| 巨乳人妻的诱惑在线观看| 久久久久久亚洲精品国产蜜桃av| 嫩草影视91久久| 国产麻豆成人av免费视频| 国产一区二区三区视频了| 久久久久精品国产欧美久久久| 日韩国内少妇激情av| 欧美日韩乱码在线| 美国免费a级毛片| 欧美日韩亚洲综合一区二区三区_| 国产精品野战在线观看| 亚洲男人天堂网一区| 久久人人爽av亚洲精品天堂| 久久精品aⅴ一区二区三区四区| 色综合欧美亚洲国产小说| 一级a爱片免费观看的视频| 国产亚洲欧美98| 国产男靠女视频免费网站| 性少妇av在线| 亚洲伊人色综图| 国产亚洲欧美98| 99国产极品粉嫩在线观看| 久久久精品国产亚洲av高清涩受| 久久人人97超碰香蕉20202| 国产极品粉嫩免费观看在线| 97碰自拍视频| 久99久视频精品免费| 亚洲精品国产区一区二| 国产精华一区二区三区| 亚洲欧美激情在线| 国产伦人伦偷精品视频| 久久精品国产综合久久久| 亚洲狠狠婷婷综合久久图片| 热re99久久国产66热| 露出奶头的视频| 色综合婷婷激情| 91av网站免费观看| 又黄又爽又免费观看的视频| 免费女性裸体啪啪无遮挡网站| 亚洲最大成人中文| 免费女性裸体啪啪无遮挡网站| 日韩国内少妇激情av| 国产成人精品无人区| 成人国产综合亚洲| 亚洲精华国产精华精| 亚洲一区中文字幕在线| av超薄肉色丝袜交足视频| 中文字幕精品免费在线观看视频| 亚洲精品中文字幕一二三四区| 手机成人av网站| 一边摸一边抽搐一进一小说| 黑人巨大精品欧美一区二区mp4| 禁无遮挡网站| 男人舔女人的私密视频| 三级毛片av免费| 亚洲精品一区av在线观看| 日本vs欧美在线观看视频| 男人舔女人下体高潮全视频| 国内精品久久久久久久电影| 男女午夜视频在线观看| 女人被狂操c到高潮| 国产不卡一卡二| 亚洲欧洲精品一区二区精品久久久| 国产精品,欧美在线| 我的亚洲天堂| 女生性感内裤真人,穿戴方法视频| 久久国产亚洲av麻豆专区| 19禁男女啪啪无遮挡网站| 国产亚洲精品av在线| 制服丝袜大香蕉在线| 手机成人av网站| 丝袜人妻中文字幕| 午夜福利高清视频| 可以在线观看的亚洲视频| av中文乱码字幕在线| 色av中文字幕| 国产精品久久久久久精品电影 | 高清在线国产一区| 欧美激情 高清一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 女性生殖器流出的白浆| 国产伦一二天堂av在线观看| 国产精品一区二区在线不卡| 免费高清视频大片| 亚洲中文日韩欧美视频| 黄频高清免费视频| 侵犯人妻中文字幕一二三四区| 两性夫妻黄色片| 此物有八面人人有两片|