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

    三角分布的海洋立管渦激振動數(shù)值分析

    2019-10-21 03:24:58趙桂峰王軍雷
    振動與沖擊 2019年19期
    關鍵詞:約化旋渦尾流

    張 猛, 劉 冉, 趙桂峰, 王軍雷

    (1. 鄭州大學 土木工程學院, 鄭州 450001; 2. 鄭州大學 化工與能源學院, 鄭州 450001)

    海洋立管是深海油氣開發(fā)系統(tǒng)中連接海面作業(yè)平臺和海底鉆采設施的關鍵設備,如鉆井隔水套管、輸液立管等,是系統(tǒng)中薄弱易損的構件之一[1]。當海流經(jīng)過立管結構時,在立管兩側會產(chǎn)生交替的旋渦脫落,柱體下游的流動分離會在其尾流中產(chǎn)生環(huán)流,根據(jù)貝努利方程,在柱體上下側會出現(xiàn)壓強差,由此產(chǎn)生橫向力,隨著旋渦脫落的交替出現(xiàn),上下橫向力也交替產(chǎn)生,從而在結構表面形成周期性的脈動作用力,很容易誘發(fā)海洋立管渦激振動[2]。若流體流速增大,旋渦脫落頻率接近于結構的固有頻率fn時,兩者會發(fā)生共振,尾渦的非定常壓力就會作用在結構上,引起結構加速運動,大大增加其動能和振幅。長時間的振動極易引起海洋立管結構疲勞破壞,會對結構安全運行構成很大威脅[3]。

    實際工程中,海洋立管往往不是單獨存在的,常以管群的方式出現(xiàn)[4],如圖1所示。

    圖1 管群結構形式

    對于典型的單圓柱鈍體結構,已引起學者們的廣泛關注,現(xiàn)階段已對單圓柱繞流及渦激振動響應進行了大量試驗與數(shù)值模擬研究[5]。而當流體流過多圓柱時其流動特性要比單圓柱復雜的多,當立管彼此相互靠近時,會發(fā)生流場干涉效應,這種情況下的激勵誘發(fā)特性與單根立管是大不相同的。目前對管群干擾效應的研究大多數(shù)是針對雙圓柱體系[6]。對于各種布置的雙圓柱繞流,Zdravkovich等[7-9]通過實驗方法研究發(fā)現(xiàn)根據(jù)圓柱間距將流體干擾可分為三類,當兩個圓柱彼此足夠靠近時會發(fā)生接近干擾,當下游圓柱部分或全部浸沒在上游圓柱尾流時會發(fā)生尾流干擾,在第三類中,接近和尾流干擾被合并。黃鈺期等[10]使用分塊耦合求解方法來數(shù)值模擬單圓柱以及不同間距下的串列雙圓柱繞流情況,研究分析了改變雙圓柱中心間距對上下游圓柱的升阻力系數(shù)和脈動頻率所產(chǎn)生的影響

    目前國內外研究學者對于管群中的三圓柱繞流問題還少有研究,三圓柱是一種具有代表性的結構排列組合單位,由于成三角形分布的三圓柱可以看成是對多對雙圓柱的組合,所以對于三圓柱繞流現(xiàn)象的研究具有重要的學術意義和工程價值[11]。基于此,本文擬對成三角形布置的等直徑三圓柱體系在均勻來流中的渦激振動特性進行研究,模擬分析來流約化速度、圓柱間距等參數(shù)對圓柱升阻力系數(shù)、振幅及尾流旋渦脫落模式的影響,以期為海洋立管管群布置提供參考依據(jù)。

    1 結構模型

    XFlow是基于格子Boltzmann粒子網(wǎng)格技術,用介觀模型來模擬流體宏觀行為,具有易于使用、高效并行、邊界條件處理簡單、模擬精確等特點[12]。本文采用XFlow進行海洋立管渦激振動計算,計算域如圖2所示。

    該流場域有5個邊界條件,進出口邊界,上下邊界和三圓柱外表面邊界。把三圓柱看成一個整體,形心到進口邊界和上下邊界的距離相同均設置為10D,形心到出口邊界的距離設置為25D。對于三個圓柱,下游兩并聯(lián)圓柱固定不動,解除上游圓柱的橫向約束,使其可以發(fā)生振動。上游圓柱到下游圓柱的水平凈距離設為L1,下游兩圓柱的垂直凈距離設為L2,本文中L1與L2始終保持相等,令L1=L2=L如圖2所示。

    圖2 計算域

    由于上游圓柱是非固定的,在流場中該圓柱在橫向力作用下會產(chǎn)生橫向單自由度渦激振動,我們把它簡化為質量-彈簧-阻尼系統(tǒng),如圖3所示。

    (a) 三維模型

    (b) 二維模型

    為了使圓柱繞流模擬結果更精確,本文采用自適應網(wǎng)格劃分方式,分別對圓柱周圍與尾流旋渦進行捕捉加密,如圖4所示。

    2 數(shù)學模型

    本文主要研究來流約化速度及圓柱間距比對圓柱繞流特性的影響,以海水為流場介質,假設密度為常數(shù),則該介質為定常不可壓縮黏性流體。

    2.1 流動控制方程

    由于流體流動要受到物理守恒定律的支配,所以需要滿足質量守恒定律(連續(xù)性方程)、動量守恒定律。

    (a) 圓柱周圍粒子網(wǎng)格加密

    (b) 圓柱及尾流區(qū)粒子網(wǎng)格捕捉加密

    (1)

    (2)

    寫成矢量形式

    (3)

    式中:μ為流體黏性系數(shù);fx,fy,fz為外力分量;ux,uy,uz為來流速度分量;ρf為流體密度;P為微元體上的壓強。

    2.2 振動方程

    單自由度質量-彈簧-阻尼系統(tǒng)振動方程為

    (4)

    (5)

    (6)

    (7)

    式中:Fl(t)、Fd(t)為圓柱所受到的升力、阻力;Cl(t)、Cd(t)為與升力所對應的升力系數(shù)、阻力系數(shù);ξ為阻尼比。各流體動力參數(shù)如表1所示。

    表1 流體動力參數(shù)

    式中:fs為漩渦脫落頻率, 一般是通過對升力的傅里葉變換求得;fn為圓柱固有頻率;D為圓柱的直徑;v為運動黏度;U∞為來流速度;m為海洋立管質量;k為系統(tǒng)剛度。立管和流體參數(shù)見表2。

    表2 立管和海水的相關參數(shù)[13]

    由于海洋立管的特點是低阻尼,低質量比,在模擬計算時根據(jù)Carmo等[14]的結果取阻尼比ξ= 0.007,將立管和海水參數(shù)取值代入質量比表達式可以求出m*= 2.5。Lee等考慮立管張力的自振頻率計算公式為[15]

    (8)

    將參數(shù)代入可求得fn= 0.76 Hz。

    3 流體求解器的驗證

    本文分別對單圓柱和三圓柱進行了模擬分析,為驗證單圓柱流固耦合計算的正確性,其各參數(shù)取值與前文所定義的海洋立管和海水參數(shù)相同, 模擬了單圓柱在雷諾數(shù)等于200、3 900、10 000時的阻力系數(shù)與斯特勞哈爾數(shù),并將結果與其他文獻中數(shù)據(jù)進行對比。如表3所示,其結果吻合較好,說明單圓柱數(shù)值模擬的精度較高。

    表3 單圓柱繞流模擬結果對比

    對于三圓柱,Gu等[25]對成三角形布置的三圓柱進行了風洞試驗研究,為驗證本文中三圓柱流固耦合計算的正確性,其模擬參數(shù)取值與試驗保持一致,圓柱直徑D=48 mm,圓柱間距比L/D=1.7,來流風速為18 m/s,空氣密度ρf=1.205 kg/m3,動力黏度系數(shù)u=1.81×10-5Pa·s,模擬分析了三個圓柱的平均阻力系數(shù)與平均升力系數(shù),并與試驗結果進行對比,如表4所示,其結果吻合較好。

    表4 數(shù)值模擬結果與風洞試驗結果對比

    4 計算結果及分析

    為了對比分析在流場中三圓柱間的相互干擾強度,首先模擬了相同參數(shù)的單圓柱渦激振動。在三圓柱渦激振動數(shù)值模擬中,圓柱間距分別根據(jù)小間距比,中間距比和大間距比取L/D= 0.5、2.0、4.0??紤]到來流約化速度對圓柱繞流的影響,在每一種圓柱間距工況下又分別選取了Ur= 1、1.5、2、2.5、3、3.5、4、4.5、5、5.5、6、6.5、7、7.5、8、8.5、9等工況進行模擬計算。本文中用CD1(CL1),CD2(CL2),CD3(CL3)分別代表上游圓柱,下游上圓柱,下游下圓柱的阻(升)力系數(shù)。

    4.1 渦激振動特性

    以小間距比L/D= 0.5 mm,Ur= 4為例,圖5~圖7中分別給出了上游圓柱的渦激振動位移時程曲線,三個圓柱的升阻力時程曲線以及各圓柱升力頻譜圖。

    從圖5可知,在上游圓柱的前幾個振動周期內,振幅逐漸增大,隨后趨于穩(wěn)定,體現(xiàn)了流體與固體的耦合作用,隨著上游圓柱振幅的增大,其對周圍流體的擾動也越來越大,流體的變化反過來又會作用在圓柱上,抑制圓柱的振動,所以圓柱的振幅不會一直增加,當流體與固體耦合穩(wěn)定時,圓柱的振幅也會趨于穩(wěn)定[26]。

    圖5 上游圓柱位移時程曲線

    圖6列舉了三個圓柱的升阻力系數(shù)時程曲線。從圖6可知,升力系數(shù)周期是阻力系數(shù)的兩倍,對于上游圓柱升力系數(shù)脈動幅值明顯大于阻力系數(shù),由此可見升力對上游圓柱的影響更大。

    將圖6各圓柱升力系數(shù)經(jīng)過傅里葉變換可得到圖7的頻譜圖,從圖7可知,三圓柱的旋渦脫落頻率均等于0.596 Hz。對比上下游圓柱的頻譜圖可以發(fā)現(xiàn),上游圓柱只有一個對應于fs= 0.596 Hz的頻率幅值突出點,但下游兩并列圓柱則有多個頻率幅值突出點,主要是因為小間距比時,上下游圓柱距離太近,使得上游圓柱后方來不及形成渦脫就直接作用在下游圓柱上,干擾下游圓柱的旋渦脫落。

    (a) 上游圓柱

    (b) 下游上圓柱

    (c) 下游下圓柱

    圖6 升阻力系數(shù)時程圖

    Fig.6 Lift-drag coefficient time history curve

    如圖8所示,分別列出了三圓柱不同間距時,上游圓柱振幅均方根和圓截面特征長度的比值Yrms/D與來流約化速度的關系曲線,并與單圓柱進行對比。從圖8可知,當來流約化速度較小時(Ur< 3)上游圓柱振幅均趨于0,激起圓柱振動不顯著,當Ur達到3.5時上游圓柱振幅均有跳躍性增加,然后隨來流約化速度的增加,上游圓柱均保持在一個較高幅值振動,當Ur進一步增大時,振幅又逐漸減小。由此可以看出L/D=0.5、2.0、4.0時其上游圓柱渦激振動的鎖振區(qū)間分別為3.5~7.0、3.5~7.5、3.5~6.0,最大振幅分別為Yrms/D= 0.372, 0.546, 0.470,單圓柱大約在3.5~5范圍內發(fā)生鎖振,最大振幅為Yrms/D= 0.472。通過對比可以發(fā)現(xiàn)圓柱間距比的大小對上游圓柱的鎖振區(qū)間和最大振幅的影響很明顯,L/D= 4.0時上游圓柱的鎖振區(qū)間、最大振幅與單圓柱很接近,證明對于大間距比,上游圓柱受下游圓柱的干擾很小。L/D= 2.0時上游圓柱的鎖振區(qū)間和最大幅值均明顯大于單圓柱,證明中間距比,下游兩并列圓柱會促進上游圓柱保持較高振幅振動。而L/D= 0.5時上游圓柱鎖振區(qū)間和最大幅值均明顯小于單圓柱,證明小間距比,下游兩并列圓柱會抑制上游圓柱的高振幅振動。所以中間距比時海洋立管振幅最大,最不安全。

    (a) 上游圓柱

    (b) 下游上圓柱

    (c) 下游下圓柱

    圖7 升力系數(shù)頻譜圖

    Fig.7 Lift coefficient spectrum diagram

    圖8 上游圓柱振幅隨來流約化速度的變化曲線

    圖9列出了不同間距比時上游圓柱渦脫頻率比fs/fn與來流約化速度Ur的關系曲線,由圖可以看出,L/D= 0.5和L/D= 2.0關系曲線變化規(guī)律相同,L/D= 4.0與單圓柱關系曲線變化規(guī)律相同。以L/D= 2.0和L/D= 4.0為例來分別介紹:當L/D= 2.0時,在Ur= 4~7.5范圍內,頻率比fs/fn≈ 1,這表明上游圓柱的渦脫頻率被鎖定在了結構的自振頻率附近,上游圓柱與流體發(fā)生了共振。當約化速度超過共振區(qū)間時,頻率比隨約化流速的增大而減小,其變化規(guī)律較單圓柱有較大差異,這是由于上游圓柱無明顯振動,上下游圓柱間距始終較小,圓柱間的接近效應和尾流效應的影響抑制了上游圓柱的旋渦脫落。當L/D= 4.0時,在Ur= 4.5~6范圍內,頻率比fs/fn≈ 1,當約化速度超過共振區(qū)間時,頻率比變化規(guī)律與單圓柱相同,均隨約化流速的增加而增大,因為對于大間距比,上下游圓柱始終可以產(chǎn)生充分發(fā)展的旋渦脫落。

    通過對上游圓柱鎖振區(qū)與共振區(qū)對比可以發(fā)現(xiàn),某些約化速度雖處于上游圓柱的鎖振區(qū)間,但其渦脫頻率并未鎖定在柱體自振頻率附近,這與本文數(shù)值模擬所采用的質量比有關,由于柱體結構與周圍流場有很強的非線性關系,可能會有附加質量對質量彈簧系統(tǒng)固有頻率和發(fā)生渦激共振時的振動頻率產(chǎn)生影響[27]。

    圖9 頻率比與來流約化速度的關系曲線

    4.2 流動特性

    如表5所示,分別給出了L/D= 0.5、2.0、4.0,來流約化速度Ur= 1.5、7、8時在同一時刻的旋渦脫落云圖,所選約化速度對應著上游圓柱渦激振動鎖振與非鎖振階段。

    從表5云圖可知,在小間距比時(L/D= 0.5),圓柱之間的干擾較為嚴重,上游圓柱后方均沒有出現(xiàn)回流區(qū)和旋渦,是由于下游兩并列圓柱限制了上游圓柱自由剪切層的卷起,使其不能產(chǎn)生旋渦脫落,而是與下游圓柱的分離層相互作用,只在下游圓柱后方出現(xiàn)規(guī)則的旋渦脫落現(xiàn)象。其流體流動特性主要受接近效應的影響。

    在中間距比時(L/D= 2.0),當約化速度在非鎖振區(qū)范圍內,在上游圓柱后方產(chǎn)生了旋渦,但下游圓柱抑制了其旋渦的脫落,上游圓柱后方的尾流被下游兩并列圓柱的尾流夾在中間,在其間隙中變窄,可見其流體流動特性受接近效應和尾流效應的共同作用。當約化速度處于鎖振區(qū)時,由于上游圓柱與流體發(fā)生了共振使其橫向振動加劇,振幅增大,這時上游圓柱與下游兩并列圓柱的間距比成周期性變化,當間距比達到一定值時,上游圓柱后方可產(chǎn)生旋渦脫落,其流動特性僅受尾流效應的影響。

    在大間距比時(L/D= 4.0),上、下游圓柱后方始終會產(chǎn)生充分發(fā)展的旋渦脫落的現(xiàn)象。在三個圓柱后方形成一個接近同步的尾流模式,來自上游圓柱的自由剪切層在下游兩并列圓柱的間隙中波動并與下游圓柱的尾流發(fā)生干涉,影響了下游圓柱的旋渦脫落,因而使得下游圓柱后方的渦脫形式變得不規(guī)則。而這種作用隨來流約化速度的增大而減小,當Ur= 8時,上下游圓柱的尾流互不干擾,周圍的流體流動模式相對于自由流方向是不對稱的,且呈現(xiàn)出明顯的偏流現(xiàn)象。

    表5 不同圓柱間距比和約化流速下圓柱渦量圖

    4.3 升阻力特性

    圓柱間流體的流動特性決定著圓柱的升阻力特性。為研究來流約化速度、圓柱間距比對各圓柱升阻力系數(shù)的影響,如圖10~圖11分別列出了L/D=0.5、2.0、4.0時的數(shù)值模擬結果。

    圖10給出了不同間距比時三圓柱阻力均值與來流約化速度的關系曲線,從圖中可以看出:上下游圓柱的阻力均值在不同間距比時隨流體約化速度的變化趨勢相同。對于上游圓柱,當Ur< 3.5時,其阻力均值隨流體約化速度的增加無明顯變化,當Ur= 3.5時發(fā)生突變,阻力均值突然大幅增加致最大值(CD,mean=0.640, 0.593, 0.566),隨后又隨流體約化速度的增加而減小。對應圖8,我們可以發(fā)現(xiàn)在Ur= 3.5時上游圓柱阻力均值的突變與振幅的突變有關。對于下游兩并列圓柱的阻力均值除了在Ur< 3.5區(qū)間內,其余均小于上游圓柱,這是由于當上游圓柱進入鎖振區(qū)間后振幅增大,渦脫頻率增加,其尾流更大程度上包裹了下游圓柱。

    圖11給出了不同間距比時三圓柱升力均方根與來流約化速度的關系曲線,從圖11可知:① 上下游圓柱的升力均方根隨來流約化速度的變化大體一致,Ur< 3.5范圍內各圓柱的升力均方根均保持較小值,在Ur= 3.5時同阻力均值一樣由于上游圓柱進入鎖振區(qū)而跳躍增加致最大值(CL,rms= 0.412, 0.430, 0.382),下游圓柱在Ur=4.0時達到最大值,隨著來流約化速度的繼續(xù)增加,各圓柱升力均方根逐漸減小;② 下游兩并列圓柱的升力均方根普遍大于上游圓柱,是由于上游圓柱后方的尾流或旋渦作用在下游圓柱表面,會對其產(chǎn)生一個附加力;③ 圓柱間距比足夠大時,上下游圓柱升力均方根僅在鎖振區(qū)間存在較小的差值,因為對于非鎖振區(qū)間,上游圓柱無明顯振動,下游兩并列圓柱間距較大,由表5可知,上游圓柱后方產(chǎn)生的旋渦直接從下游圓柱間隙中流過,不能充分作用在下游圓柱表面。對于鎖振區(qū)間,由于上游圓柱發(fā)生大幅振動,其后方產(chǎn)生的旋渦會周期性的作用在下游圓柱表面,且由于上下游圓柱間距比較大,上游圓柱產(chǎn)生的旋渦需經(jīng)過一段長距離才能作用于下游圓柱,其能量有較大的損耗。由此可以得出,當圓柱達到大間距比時,隨圓柱間距比的增大,上下游圓柱的升力系數(shù)差值會逐漸減小;④ 在小間距比的鎖振區(qū)間及中間距比的超鎖振區(qū)間,下游兩并列圓柱的升力均方根出現(xiàn)明顯差異,這是由于上游圓柱后方的尾流是不對稱的。

    (a) L/D = 0.5

    (b) L/D = 2.0

    (c) L/D = 4.0

    5 結 論

    對阻尼比ξ= 0.007,質量比m*= 2.5,三角形分布的海洋立管在間距比為L/D= 0.5、2.0、4.0,來流約化速度為Ur= 1~9進行了渦激振動數(shù)值模擬,并與單圓柱渦激振動進行了對比,將結果歸納如下:

    (1) 海洋立管間距比的大小對上游立管的鎖振區(qū)間和最大振幅有很大影響,中間距比時其鎖振區(qū)間和最大振幅最大,最不安全。

    (a) L/D = 0.5

    (b) L/D = 2.0

    (c) L/D = 4.0

    (2) 由三圓柱渦量云圖對比可知,立管間距比較大或較小時,立管間流體的流動模式分別受尾流效應和接近效應的影響,中間距比時受兩種效應的共同作用。

    (3) 由于上游立管的尾流或旋渦會對下游立管產(chǎn)生附加力,使得上游立管的升力系數(shù)普遍低于下游立管,且當立管達到大間距比時,隨立管間距比的增大,上下游立管的升力系數(shù)差值會逐漸減小。在小間距比的鎖振區(qū)間和中間距比的超鎖振區(qū)間,下游兩并列立管表現(xiàn)出明顯的不對稱性,這是由于上游立管后方流體流動的不對稱性,尾流區(qū)出現(xiàn)了偏流現(xiàn)象。

    猜你喜歡
    約化旋渦尾流
    約化的(3+1)維Hirota方程的呼吸波解、lump解和半有理解
    小心,旋渦來啦
    大班科學活動:神秘的旋渦
    旋渦笑臉
    山間湖
    揚子江(2019年1期)2019-03-08 02:52:34
    飛機尾流的散射特性與探測技術綜述
    雷達學報(2017年6期)2017-03-26 07:53:06
    錐形流量計尾流流場分析
    水面艦船風尾流效應減弱的模擬研究
    M-強對稱環(huán)
    (3+1)-維廣義Kadomtsev-Petviashvili方程的對稱約化與精確解
    日本撒尿小便嘘嘘汇集6| 久久久欧美国产精品| 赤兔流量卡办理| 国产成人一区二区在线| 亚洲欧美成人精品一区二区| 天天一区二区日本电影三级| 国产成人精品久久久久久| 国产亚洲精品久久久久久毛片| 人妻夜夜爽99麻豆av| 国产视频一区二区在线看| 高清午夜精品一区二区三区 | 欧美人与善性xxx| 国产伦在线观看视频一区| 天天一区二区日本电影三级| 精品少妇黑人巨大在线播放 | 一级黄片播放器| 午夜福利成人在线免费观看| 久久午夜亚洲精品久久| 插阴视频在线观看视频| 3wmmmm亚洲av在线观看| 国产午夜精品久久久久久一区二区三区 | 午夜日韩欧美国产| 91在线观看av| 少妇熟女aⅴ在线视频| 成年女人看的毛片在线观看| 日产精品乱码卡一卡2卡三| 久久精品国产自在天天线| 两个人的视频大全免费| 变态另类成人亚洲欧美熟女| 在线观看66精品国产| 国产极品精品免费视频能看的| 在线观看av片永久免费下载| 一级a爱片免费观看的视频| 亚洲乱码一区二区免费版| av天堂在线播放| 国产一区二区在线观看日韩| 国产精品,欧美在线| 国产精品日韩av在线免费观看| 欧美性猛交╳xxx乱大交人| 一个人免费在线观看电影| 插阴视频在线观看视频| 熟女人妻精品中文字幕| 午夜激情欧美在线| 午夜久久久久精精品| 大又大粗又爽又黄少妇毛片口| 在线国产一区二区在线| 免费看光身美女| 黄色日韩在线| 日韩人妻高清精品专区| 国产男靠女视频免费网站| 国产精品av视频在线免费观看| 国产伦精品一区二区三区视频9| 一a级毛片在线观看| 亚洲精品日韩av片在线观看| 国产 一区精品| 我的女老师完整版在线观看| 舔av片在线| 亚洲婷婷狠狠爱综合网| 在线看三级毛片| 精品99又大又爽又粗少妇毛片| 成年女人毛片免费观看观看9| 女的被弄到高潮叫床怎么办| 成人美女网站在线观看视频| 亚洲av中文字字幕乱码综合| 国产白丝娇喘喷水9色精品| av天堂中文字幕网| 你懂的网址亚洲精品在线观看 | 成人av在线播放网站| 午夜激情欧美在线| 午夜福利成人在线免费观看| 国产色爽女视频免费观看| 中国美白少妇内射xxxbb| 小说图片视频综合网站| 国产免费一级a男人的天堂| 少妇被粗大猛烈的视频| 国模一区二区三区四区视频| 青春草视频在线免费观看| 国产精品一及| 精品午夜福利在线看| 亚洲av免费高清在线观看| 天堂网av新在线| 国产女主播在线喷水免费视频网站 | 久久亚洲精品不卡| 亚洲五月天丁香| 精品人妻一区二区三区麻豆 | 一个人观看的视频www高清免费观看| 国产91av在线免费观看| 久久精品91蜜桃| 亚洲熟妇中文字幕五十中出| 亚洲成人久久爱视频| 啦啦啦韩国在线观看视频| 国产精品综合久久久久久久免费| 丰满乱子伦码专区| 国产精品久久久久久久久免| 九九在线视频观看精品| 色哟哟·www| 99热6这里只有精品| 国产成人一区二区在线| 午夜福利18| 男女之事视频高清在线观看| 观看免费一级毛片| 亚洲av中文av极速乱| 国国产精品蜜臀av免费| 欧美成人精品欧美一级黄| 婷婷亚洲欧美| 国产私拍福利视频在线观看| 久久精品国产亚洲网站| 天天一区二区日本电影三级| 亚洲av美国av| 乱系列少妇在线播放| 波多野结衣巨乳人妻| 2021天堂中文幕一二区在线观| 国产黄色小视频在线观看| 校园春色视频在线观看| 国产成人精品久久久久久| 国产一区二区激情短视频| 麻豆乱淫一区二区| 国语自产精品视频在线第100页| 免费电影在线观看免费观看| 简卡轻食公司| 久久久久国产精品人妻aⅴ院| 搡女人真爽免费视频火全软件 | 亚洲欧美日韩东京热| 色噜噜av男人的天堂激情| 一级毛片我不卡| 黄色欧美视频在线观看| 日产精品乱码卡一卡2卡三| 日日干狠狠操夜夜爽| 又爽又黄a免费视频| 国产一级毛片七仙女欲春2| a级毛片免费高清观看在线播放| a级毛色黄片| 久久久久国产网址| 免费av毛片视频| 最近在线观看免费完整版| 人人妻人人看人人澡| av国产免费在线观看| videossex国产| 亚洲欧美日韩高清在线视频| 亚洲成人中文字幕在线播放| 禁无遮挡网站| 国产精品人妻久久久久久| 18禁黄网站禁片免费观看直播| 午夜福利在线观看吧| 欧美一级a爱片免费观看看| 嫩草影视91久久| 精品免费久久久久久久清纯| 亚洲欧美日韩高清在线视频| 村上凉子中文字幕在线| 嫩草影院入口| 中国美白少妇内射xxxbb| 一区福利在线观看| 熟女人妻精品中文字幕| 亚洲av二区三区四区| 日韩三级伦理在线观看| 久久久国产成人精品二区| 精品福利观看| 国国产精品蜜臀av免费| 直男gayav资源| 国产免费一级a男人的天堂| 听说在线观看完整版免费高清| 人妻久久中文字幕网| 国产午夜福利久久久久久| 91久久精品国产一区二区三区| 亚洲精品456在线播放app| 变态另类成人亚洲欧美熟女| 综合色av麻豆| 精品久久久久久久久久久久久| 国产免费一级a男人的天堂| 美女xxoo啪啪120秒动态图| 蜜桃久久精品国产亚洲av| 身体一侧抽搐| 99国产精品一区二区蜜桃av| 看片在线看免费视频| 美女大奶头视频| 一级毛片我不卡| 少妇的逼水好多| 黄片wwwwww| 亚洲最大成人手机在线| 少妇熟女aⅴ在线视频| 亚洲精品色激情综合| 久久精品国产亚洲网站| 欧美区成人在线视频| 久久久久久久午夜电影| 欧美国产日韩亚洲一区| 久久热精品热| 亚洲国产精品成人综合色| 亚洲在线观看片| 午夜影院日韩av| 精品免费久久久久久久清纯| 国产91av在线免费观看| 国内少妇人妻偷人精品xxx网站| 国国产精品蜜臀av免费| 在线看三级毛片| 亚洲精品一卡2卡三卡4卡5卡| 国产精品三级大全| 亚洲天堂国产精品一区在线| 波多野结衣高清无吗| 性插视频无遮挡在线免费观看| 欧美潮喷喷水| 免费av毛片视频| 日本撒尿小便嘘嘘汇集6| 小说图片视频综合网站| av在线播放精品| 国产探花在线观看一区二区| 草草在线视频免费看| 国产真实乱freesex| 久久久久国产精品人妻aⅴ院| 久久人人爽人人爽人人片va| 极品教师在线视频| 欧美+亚洲+日韩+国产| 六月丁香七月| 亚洲欧美成人精品一区二区| 亚洲熟妇熟女久久| 别揉我奶头 嗯啊视频| 日韩一本色道免费dvd| 一个人观看的视频www高清免费观看| av专区在线播放| 久久精品久久久久久噜噜老黄 | 免费观看的影片在线观看| 国产精品不卡视频一区二区| 天堂√8在线中文| 老司机福利观看| av在线蜜桃| 丝袜美腿在线中文| 搡老熟女国产l中国老女人| 久久久久久久久久黄片| 久久精品国产亚洲av天美| 日日啪夜夜撸| 毛片女人毛片| 免费看光身美女| 麻豆国产97在线/欧美| 联通29元200g的流量卡| 草草在线视频免费看| 韩国av在线不卡| 在线天堂最新版资源| 国产精品女同一区二区软件| 亚洲av成人精品一区久久| 黄片wwwwww| 亚洲人与动物交配视频| 成人亚洲欧美一区二区av| 久久久精品欧美日韩精品| 精品久久久久久成人av| 国产午夜精品久久久久久一区二区三区 | 老女人水多毛片| 国产三级在线视频| 精品欧美国产一区二区三| 久久久久国产网址| 免费观看人在逋| 最近最新中文字幕大全电影3| 人人妻人人澡欧美一区二区| 简卡轻食公司| 身体一侧抽搐| 国产精华一区二区三区| .国产精品久久| 成人漫画全彩无遮挡| 国产成人a区在线观看| 五月伊人婷婷丁香| 日韩 亚洲 欧美在线| 免费av不卡在线播放| 久久久精品欧美日韩精品| 精品一区二区三区av网在线观看| 免费av毛片视频| 老司机午夜福利在线观看视频| 亚洲欧美精品自产自拍| 天天躁日日操中文字幕| 尤物成人国产欧美一区二区三区| 22中文网久久字幕| 亚洲av免费高清在线观看| 国产精品电影一区二区三区| 精品人妻视频免费看| 国产麻豆成人av免费视频| 国产高清不卡午夜福利| 天天躁日日操中文字幕| 亚洲av五月六月丁香网| 综合色丁香网| av在线亚洲专区| 天堂√8在线中文| 精品人妻一区二区三区麻豆 | 久久精品国产亚洲网站| 又粗又爽又猛毛片免费看| 俺也久久电影网| 人妻制服诱惑在线中文字幕| 99九九线精品视频在线观看视频| 免费av观看视频| 国产精品一区二区三区四区久久| 国产精品一区二区免费欧美| 国产精品国产三级国产av玫瑰| 精品午夜福利视频在线观看一区| 成人三级黄色视频| 成人二区视频| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 国产精品一区二区性色av| 97热精品久久久久久| 久久鲁丝午夜福利片| 露出奶头的视频| 天堂网av新在线| 香蕉av资源在线| 成人亚洲精品av一区二区| 99国产极品粉嫩在线观看| 老司机影院成人| 欧美日韩综合久久久久久| 久久6这里有精品| 久久亚洲国产成人精品v| 变态另类丝袜制服| 国产精品电影一区二区三区| 九色成人免费人妻av| 久久人妻av系列| 在线观看一区二区三区| 99热这里只有精品一区| 男女下面进入的视频免费午夜| 亚洲一区高清亚洲精品| 日韩欧美精品v在线| 午夜激情欧美在线| 久久欧美精品欧美久久欧美| 久久人人爽人人片av| 悠悠久久av| 国产亚洲精品久久久com| h日本视频在线播放| 免费观看人在逋| 在线播放无遮挡| 久久久精品94久久精品| 午夜免费男女啪啪视频观看 | 成人高潮视频无遮挡免费网站| 乱人视频在线观看| 91狼人影院| 激情 狠狠 欧美| 精品久久久久久久人妻蜜臀av| 久久人人爽人人爽人人片va| 免费人成在线观看视频色| 日韩三级伦理在线观看| 国产色婷婷99| 女生性感内裤真人,穿戴方法视频| 亚洲av成人精品一区久久| 高清日韩中文字幕在线| 国产综合懂色| 在线免费观看的www视频| 久久精品国产亚洲av天美| 久久精品久久久久久噜噜老黄 | 久久久久性生活片| 国产三级在线视频| 久久精品国产99精品国产亚洲性色| 久久欧美精品欧美久久欧美| 国产视频内射| 亚洲成人精品中文字幕电影| 欧美中文日本在线观看视频| 亚洲国产精品sss在线观看| 欧美丝袜亚洲另类| 嫩草影视91久久| 熟女人妻精品中文字幕| 亚洲av中文av极速乱| 日韩 亚洲 欧美在线| 午夜福利高清视频| 一进一出抽搐动态| 此物有八面人人有两片| 成年免费大片在线观看| 国产真实伦视频高清在线观看| 久久久久性生活片| 国产精品无大码| 日本黄色片子视频| 97人妻精品一区二区三区麻豆| 美女高潮的动态| 久久人人精品亚洲av| 最近的中文字幕免费完整| 高清午夜精品一区二区三区 | 性插视频无遮挡在线免费观看| 18+在线观看网站| 成人特级av手机在线观看| 可以在线观看毛片的网站| 长腿黑丝高跟| 精品久久久久久久人妻蜜臀av| 一级毛片电影观看 | 非洲黑人性xxxx精品又粗又长| 22中文网久久字幕| 麻豆精品久久久久久蜜桃| 人人妻人人澡欧美一区二区| 99久国产av精品| 精品久久久久久久久亚洲| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产自在天天线| 亚洲久久久久久中文字幕| 亚洲美女视频黄频| 观看免费一级毛片| 99在线人妻在线中文字幕| 国产av麻豆久久久久久久| 色综合站精品国产| 久久久久免费精品人妻一区二区| 超碰av人人做人人爽久久| 久久久久免费精品人妻一区二区| 亚洲av第一区精品v没综合| 欧美激情久久久久久爽电影| 国产精品久久视频播放| 亚洲第一区二区三区不卡| 欧美色视频一区免费| 国产高清视频在线播放一区| 在线免费观看不下载黄p国产| 蜜臀久久99精品久久宅男| 国产日本99.免费观看| 亚洲美女黄片视频| 成年av动漫网址| 久久精品国产亚洲av香蕉五月| av黄色大香蕉| 五月玫瑰六月丁香| 久久久久久久久中文| 黄色视频,在线免费观看| 欧美绝顶高潮抽搐喷水| 卡戴珊不雅视频在线播放| 网址你懂的国产日韩在线| 成人无遮挡网站| 日本黄色视频三级网站网址| 能在线免费观看的黄片| 午夜福利18| 人妻久久中文字幕网| 国产伦精品一区二区三区视频9| 在线国产一区二区在线| 成人漫画全彩无遮挡| 欧美+日韩+精品| 狠狠狠狠99中文字幕| 网址你懂的国产日韩在线| 18禁裸乳无遮挡免费网站照片| 欧美日韩一区二区视频在线观看视频在线 | 一区二区三区四区激情视频 | 日韩中字成人| 中文在线观看免费www的网站| 国产精品国产三级国产av玫瑰| 日韩大尺度精品在线看网址| 成人特级av手机在线观看| 国产午夜精品论理片| 两个人视频免费观看高清| 看十八女毛片水多多多| 日本欧美国产在线视频| 此物有八面人人有两片| 色视频www国产| 国产一区二区在线观看日韩| 最近最新中文字幕大全电影3| 国产国拍精品亚洲av在线观看| 可以在线观看毛片的网站| 我要看日韩黄色一级片| 国内精品一区二区在线观看| 日韩欧美精品免费久久| 午夜福利高清视频| 久久草成人影院| 久久亚洲精品不卡| 午夜精品一区二区三区免费看| 色综合站精品国产| 久久鲁丝午夜福利片| 97热精品久久久久久| 国产免费一级a男人的天堂| 欧美又色又爽又黄视频| 久久久久性生活片| 精品国产三级普通话版| 亚洲人成网站高清观看| 乱人视频在线观看| 国语自产精品视频在线第100页| 久久国内精品自在自线图片| 午夜福利在线在线| 色av中文字幕| 亚州av有码| 麻豆一二三区av精品| 在线观看美女被高潮喷水网站| 成年版毛片免费区| 麻豆国产av国片精品| 久久精品91蜜桃| 国产亚洲av嫩草精品影院| 亚洲精品亚洲一区二区| 特大巨黑吊av在线直播| 国产精品,欧美在线| 日韩一区二区视频免费看| 最近视频中文字幕2019在线8| 日本熟妇午夜| 毛片女人毛片| 国产精品福利在线免费观看| 国产久久久一区二区三区| 国产精品人妻久久久影院| 国产精品人妻久久久久久| 日日摸夜夜添夜夜爱| 国产熟女欧美一区二区| 秋霞在线观看毛片| 成人国产麻豆网| 日韩大尺度精品在线看网址| 少妇裸体淫交视频免费看高清| .国产精品久久| 国产片特级美女逼逼视频| 日本与韩国留学比较| 日本在线视频免费播放| 国产一区二区亚洲精品在线观看| 少妇熟女aⅴ在线视频| 午夜福利视频1000在线观看| 波多野结衣高清无吗| 成年av动漫网址| 天堂√8在线中文| 欧美+日韩+精品| 午夜老司机福利剧场| 麻豆久久精品国产亚洲av| 日本撒尿小便嘘嘘汇集6| 草草在线视频免费看| 一级毛片aaaaaa免费看小| 日本免费a在线| 精品不卡国产一区二区三区| 欧美xxxx性猛交bbbb| 少妇丰满av| 大又大粗又爽又黄少妇毛片口| 在线观看美女被高潮喷水网站| 国产老妇女一区| 一级黄片播放器| 男女做爰动态图高潮gif福利片| 一级毛片aaaaaa免费看小| 日本爱情动作片www.在线观看 | aaaaa片日本免费| 亚洲欧美成人综合另类久久久 | 久久久a久久爽久久v久久| 嫩草影院精品99| 国产黄片美女视频| 国产麻豆成人av免费视频| 日韩精品有码人妻一区| 99精品在免费线老司机午夜| 久久6这里有精品| 亚洲成人久久性| 欧美成人免费av一区二区三区| 国产精品一区二区三区四区久久| 久久精品国产亚洲av涩爱 | 黄色欧美视频在线观看| 69人妻影院| 三级国产精品欧美在线观看| 久久久国产成人精品二区| 久久久久久久久久黄片| 国产伦在线观看视频一区| 最近手机中文字幕大全| 亚洲自拍偷在线| 久久久久久久亚洲中文字幕| 少妇裸体淫交视频免费看高清| 日本黄色视频三级网站网址| 欧美中文日本在线观看视频| 天堂动漫精品| 免费看光身美女| 国产欧美日韩一区二区精品| 少妇被粗大猛烈的视频| 美女被艹到高潮喷水动态| 少妇猛男粗大的猛烈进出视频 | 久久久精品欧美日韩精品| 欧美日韩综合久久久久久| 在线看三级毛片| 99久久中文字幕三级久久日本| 久久久久免费精品人妻一区二区| 美女cb高潮喷水在线观看| 亚洲自拍偷在线| 成人美女网站在线观看视频| 日日干狠狠操夜夜爽| 国产高清视频在线播放一区| 欧美高清性xxxxhd video| 久久久久久国产a免费观看| 国产精品久久久久久亚洲av鲁大| 国产成人精品久久久久久| 丰满的人妻完整版| 精品久久久久久久久久免费视频| 亚洲一级一片aⅴ在线观看| 国产高清三级在线| 99九九线精品视频在线观看视频| 色综合色国产| 亚洲精品国产av成人精品 | 午夜精品一区二区三区免费看| 国产精品精品国产色婷婷| 成人亚洲欧美一区二区av| 亚洲国产色片| 精品国产三级普通话版| 国产av麻豆久久久久久久| 成人漫画全彩无遮挡| 一级黄色大片毛片| 麻豆乱淫一区二区| 日韩欧美一区二区三区在线观看| 久久久欧美国产精品| 久久久久国内视频| 欧美日韩在线观看h| 黄色欧美视频在线观看| av视频在线观看入口| 成人毛片a级毛片在线播放| 极品教师在线视频| 亚洲内射少妇av| 久久精品综合一区二区三区| 亚洲国产精品合色在线| 亚洲欧美中文字幕日韩二区| 搡女人真爽免费视频火全软件 | 久久精品国产清高在天天线| 中文字幕av在线有码专区| 99久久无色码亚洲精品果冻| 日韩一本色道免费dvd| 99久久成人亚洲精品观看| 女生性感内裤真人,穿戴方法视频| 免费无遮挡裸体视频| 欧美xxxx性猛交bbbb| 久久婷婷人人爽人人干人人爱| 亚洲人成网站在线观看播放| 成人无遮挡网站| 亚洲美女黄片视频| 亚洲在线自拍视频| 亚洲性夜色夜夜综合| 成人特级黄色片久久久久久久| 欧美中文日本在线观看视频| 五月玫瑰六月丁香| 黄色配什么色好看| 最近2019中文字幕mv第一页| 丰满的人妻完整版| 国产免费男女视频| 国产真实乱freesex| 久久这里只有精品中国| 人人妻人人看人人澡| 亚洲欧美清纯卡通| 97人妻精品一区二区三区麻豆| av专区在线播放| 久久天躁狠狠躁夜夜2o2o| 午夜影院日韩av| 国内精品久久久久精免费| 免费搜索国产男女视频| 午夜福利在线在线|