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

    并列雙圓柱渦激振動的經(jīng)驗性模型研究

    2013-11-22 00:59:30孫麗萍宋儒鑫
    海洋工程 2013年2期
    關(guān)鍵詞:渦激尾流振子

    秦 偉,康 莊,孫麗萍,宋儒鑫

    (哈爾濱工程大學 船舶工程學院深海工程技術(shù)研究中心,黑龍江 哈爾濱 150001)

    由于圓柱體之間的相互影響,雙圓柱體的渦激振動機理及現(xiàn)象相對于單圓柱體而言更加復(fù)雜。目前,國內(nèi)外的研究者對雙圓柱體的渦激振動機理和現(xiàn)象進行了一定的研究,發(fā)現(xiàn)了影響串列或并列排列雙圓柱渦激振動特性的主要參數(shù),比如兩圓柱體的中心間距、雷諾數(shù)以及兩個圓柱體的約束方式等。

    Zdravkovich[1]對雙圓柱的渦激振動受力特性以及尾流形式做了綜述性的介紹,同時指出:對于雙圓柱體,由于兩個圓柱周圍流場的相互影響作用,在同一雷諾數(shù)下,也會表現(xiàn)出和單個孤立圓柱體不同的振動特性和旋渦脫落現(xiàn)象,而且隨著橫向和流向間距比的不同,兩個圓柱體的尾渦特性和受力特性也會有較大的不同。姚熊亮[2]對均勻流中并列剛性和彈性圓柱的流場進行了分類,同時認為脈動升力和阻力系數(shù)與柱的振幅、約化阻尼及兩柱間距有關(guān),對柱型結(jié)構(gòu)而言,合理的響應(yīng)預(yù)報模型取決于流體載荷的系統(tǒng)分析。Williamson[3]對并列雙圓柱體的尾流特性做了深入的研究,發(fā)現(xiàn)兩個圓柱體在T/D >2.2 的特定間距比條件下,并列雙圓柱體下游的尾渦域內(nèi)存在平行渦街同步的現(xiàn)象,即兩個平行渦街表現(xiàn)為反相位(in anti-phase)或者是同相位(in phase)這兩種形式。Sumner[4]也對并列雙圓柱體尾渦進行了研究,發(fā)現(xiàn)其尾渦模態(tài)可能存在3 種形式,分別為單一渦街(Single-bluff-body behaviors)、偏向流動(Biased flow pattern)和平行雙卡門渦街(Parallel vortex streets)形式。具體每個影響區(qū)域的中心間距會隨著雷諾數(shù)、實驗條件和固定方式等的不同而發(fā)生一定的變化。Wang[5]發(fā)現(xiàn)當并列排列的兩個圓柱體之間的中心間距比為1 <T/D <1.1 ~1.2 時,兩個圓柱體可以基本看成一個圓柱體,出現(xiàn)單個卡門渦街。Sumner[4]研究發(fā)現(xiàn)在雷諾數(shù)Re=1 000 ~3 000,中心間距比T/D=1.5 時,出現(xiàn)了偏向流動的現(xiàn)象,圓柱體出現(xiàn)不對稱的尾流運動和高頻率的旋渦脫落模態(tài),其中一側(cè)的尾渦脫落頻率是另一側(cè)的脫落頻率的倍數(shù)關(guān)系,同時Sumner 還指出在T/D =2.5 和T/D =4.5 時發(fā)現(xiàn)了并排雙卡門渦街現(xiàn)象,此條件下兩個圓柱體相互影響比較小,同時兩圓柱體的旋渦脫落頻率相當。

    由于渦致雙圓柱振動的過程比單圓柱情況更加復(fù)雜,公開發(fā)表的有關(guān)預(yù)報雙圓柱渦激振動響應(yīng)的經(jīng)驗性模型十分有限。Facchinetti 等[6]采用尾流振子模型對串列排列的雙圓柱渦激振動進行描述,得到的計算結(jié)果能夠在定性和定量上與實驗和數(shù)值模擬結(jié)果保持一致。

    本文主要研究二維并列雙圓柱的渦激振動,同時建立一個能夠預(yù)測結(jié)構(gòu)渦激振動響應(yīng)的經(jīng)驗性模型:由于特定間距比條件下的旋渦脫落頻率出現(xiàn)分支,因此假設(shè)用兩個具有不同固有頻率(即斯特勞哈爾頻率)的尾流振子來描述流體對結(jié)構(gòu)的非線性動力作用,同時兩個流體振子均滿足van der Pol 方程,進而得到了結(jié)構(gòu)振子和流體振子的耦合方程組,聯(lián)立組成并列雙圓柱渦激振動的尾流振子模型。采用實例基本數(shù)據(jù),對經(jīng)驗性模型進行了數(shù)值求解,并與實驗測量值進行對比分析,結(jié)果表明模型預(yù)報的響應(yīng)能夠很好地預(yù)計并列雙圓柱渦激振動的結(jié)構(gòu)振動特性和一般規(guī)律。

    1 并列圓柱尾流振子模型

    1.1 結(jié)構(gòu)振動方程

    圖1 橫向并列圓柱布置Fig.1 A schematic layout of two side-by-side cylinders in a cross-flow

    二維并列雙圓柱系統(tǒng)的布置如圖1 所示。兩個圓柱的直徑均為D,圓心之間的間距為T,同時考慮在流速為U∞的均勻來流作用下,圓柱沿垂直流向的橫向(Y 軸)運動。Facchinetti 等[7]給出了孤立彈性支承剛性圓柱的結(jié)構(gòu)振動方程,對其進行拓展可以得到雙圓柱系統(tǒng)的結(jié)構(gòu)振動方程為:

    式中:m=(ms+ma),ms為結(jié)構(gòu)的質(zhì)量,ma為流體的附加質(zhì)量表示無粘流體的慣性作用力;ma=CaπD2ρ/4,Ca為流體的附加質(zhì)量系數(shù);r =(rs+rf),為結(jié)構(gòu)粘性阻尼;rf為流體附加阻尼,表示流體的粘性作用力;k 為結(jié)構(gòu)的剛度;Y1和Y2分別表示圓柱1 和圓柱2 的橫向位置坐標,分別為結(jié)構(gòu)橫向運動速度和加速度;FVY1和FVY2分別表示脫落的旋渦對圓柱1 和圓柱2 的脈動升力,通常將升力定義為無量綱的升力系數(shù)形式,即FVY1=ρU2DCL1/2 和FVY2=ρU2DCL2/2。

    1.2 尾流振子方程

    尾流振子方程將結(jié)構(gòu)近壁的尾渦動力特性通過非線性振子來模擬,例如采用van der Pol 方程或者Rayleigh 方程來描述其自激自限的諧振動過程,并與結(jié)構(gòu)振子聯(lián)立構(gòu)成方程組,共同預(yù)測流固耦合系統(tǒng)的響應(yīng)。Bishop 和Hassan[8]首先提出采用一個非線性自激的振子來表示作用在圓柱體上的升力。Hartlen 和Currie[9]以升力系數(shù)作為van der Pol 方程的狀態(tài)變量,給出了尾流振子方程的數(shù)學表達式。國內(nèi)外諸多學者補充并推廣該方法,主要應(yīng)用于研究孤立圓柱體,海洋立管和管道的渦激振動,如郭海燕[10]。

    根據(jù)尾流振子模型建立的基本方法,引入無量綱參數(shù)q =2CL/CL0(其中的CL0為參考脈動升力系數(shù)[6])為尾流振子方程變量,同時滿足受強迫激勵的van der Pol 方程,其中強迫項表示結(jié)構(gòu)對流體振子的耦合作用。與單圓柱情況不同,并列圓柱的間距比對斯特勞哈爾數(shù)可能存在影響,圖2 給出了斯特勞哈爾數(shù)隨間距比變化的函數(shù)關(guān)系,可以看到尾渦模態(tài)為偏向流動情況下,可能存在兩個不同的主斯特勞哈爾頻率。因此,假設(shè)在特定的T/D 范圍內(nèi),流場中存在兩個具有不同固有頻率的尾流振子,其表達為:

    式中:ε 為van der Pol 方程參數(shù);q1和q2分別表示圓柱1 和圓柱2 的尾流振子變量;ωs1和ωs2分別為圓柱1和圓柱2 相對應(yīng)的旋渦脫落頻率,同時也是流體振子的固有頻率,其與斯特勞哈爾數(shù)的關(guān)系表達為ωs1,2=分別為耦合項。

    Facchinetti 等[7]充分研究了單圓柱情況下結(jié)構(gòu)與流體振子的耦合作用,分別對位移、速度、加速度三種耦合形式進行了對比和討論。研究表明強迫項采用結(jié)構(gòu)的加速度形式最為理想,能夠反映共振特性,因此本文中采用的形式寫為:

    其中:A1和A2為實驗參數(shù)。

    因此,脈動升力系數(shù)可以表示為:

    其中:(α1,β1,α2,β2)為模型的經(jīng)驗性參數(shù),分別是間距比T/D 的函數(shù)并滿足:

    1.3 耦合振子模型

    通過對結(jié)構(gòu)振子方程式(1a)和式(1b)以及尾流振子方程式(2a)和式(2b)進行無量綱化,同時分別引入無量綱時間τ=(k/m)1/2t 和無量綱坐標y =Y/D,最終得到橫向并列雙圓柱的衰減階(Reduced-Order)耦合振子模型為:

    其中:此處的導(dǎo)數(shù)理解為對無量綱時間τ 求導(dǎo);(α1,β1,α2,β2,a1,a2)為引入的無量綱參數(shù);λ =r/(mk)為總阻尼比;ωn= (k/m)0.5為結(jié)構(gòu)的固有頻率;μ = m/ρD2為質(zhì)量比;ur= U/fnD 為來流的約化速度;ω1,2=ωs1,2/ωn為頻率比。

    模型中的部分參數(shù)可以借鑒Facchinetti 等[7]尾流振子模型參數(shù)確定,本文中取a1=6,a2=6,ε =0.3。尾流振子方程中的頻率ω1和ω2分別代表旋渦的脫落頻率,需要根據(jù)不同尾渦模態(tài)中St隨間距比的變化規(guī)律來確定。參考脈動升力系數(shù)CL01和CL02一般是在固定圓柱條件下測得的,通常與雷諾數(shù)Re 有關(guān)。模型的經(jīng)驗性參數(shù)(α1,β1,α2,β2)需要根據(jù)已有大量實驗數(shù)據(jù)通過曲線擬合或最小二乘法確定其函數(shù)變化規(guī)律,由于現(xiàn)階段可用的實驗數(shù)據(jù)有限,因此這些參數(shù)在具體實驗中給出。

    圖2 斯特勞哈爾數(shù)與間距比T/D 的關(guān)系Fig.2 Strouhal number data as a function of transverse pitch ratio T/D

    2 計算實例與結(jié)果分析

    2.1 實例基本數(shù)據(jù)

    文中對比分析的數(shù)據(jù)來源于Liu 等[11]中高質(zhì)量比并列雙圓柱的數(shù)值模擬實驗結(jié)果和低質(zhì)量比的并列排列雙圓柱的模型實驗測量結(jié)果。

    Liu 等對不同間距比條件下的并列雙圓柱進行了有限元CFD 數(shù)值模擬。分析中選取了三組不同的間距比,分別為T/D=1.8、3.0 和1.1。前兩種工況的設(shè)計參數(shù)為:無量綱質(zhì)量比m*=4ms/(ρπD2)=10,質(zhì)量阻尼比m*ξ=0.3。而在工況T/D=1.1 條件下,兩圓柱之間的距離很小,為了避免圓柱相撞需要調(diào)整結(jié)構(gòu)的質(zhì)量,其相應(yīng)的設(shè)計參數(shù)取為無量綱質(zhì)量比為m*=66,質(zhì)量阻尼比為m*ξ =2.0。約化速度均控制在ur=5.0,相應(yīng)的雷諾數(shù)為Re=200。

    低質(zhì)量比并列雙圓柱的模型試驗在船模拖曳水池開展。船模拖曳水池長108 m,池寬7 m,水深3.5 m。拖車車速0 ~6.5 m/s,拖車測量橋在垂直方向上的可調(diào)行程為0.4 m。實驗中使用的管模型為圓柱體,采用的材料為硬聚氯乙烯(UPVC),直徑10 cm,壁厚5 mm,管長1.0 m,管底部有配重,兩端用水密封頭密封。模型的輔助設(shè)備包括與拖車固定的框架結(jié)構(gòu)和與圓柱管相連的可擺動連桿架結(jié)構(gòu)組成,圖3 給出了實驗?zāi)P偷难b置圖和三維仿真模擬圖。實驗的相應(yīng)設(shè)計參數(shù)包括,無量綱質(zhì)量比m*=1.4,質(zhì)量阻尼比m*ξ =0.047 6。約化速度的范圍在ur=4.2 ~16.8,相應(yīng)的雷諾數(shù)范圍為Re=2 ×104~8 ×104。

    圖3 并列雙圓柱實驗?zāi)P秃蛯嶒炄S仿真模擬Fig.3 The experimental model for two side-by-side cylinders and 3-D simulation

    2.2 模型的數(shù)值求解

    并列雙圓柱的尾流振子模型即式(6a)~(6d)的數(shù)值求解使用的是自適應(yīng)步長的四階龍格-庫塔方法。為了得到無量綱振幅與約化速度的關(guān)系,分別對增加流速和減小流速兩種情況進行計算,其流速增量的大小取為0.01,同時需要滿足后一流速條件下的初始條件取為前一流速的最終結(jié)果。時域模擬的時間控制為至少50 個振動周期,取穩(wěn)定后的結(jié)果進行傅里葉分析從而得到振動的幅值與頻率。

    Liu 等[11]首先對不同間距比條件下的固定并列雙圓柱進行了數(shù)值模擬實驗,實驗中可以測得參考脈動升力系數(shù)(CL0)和旋渦脫落頻率,為本模型中經(jīng)驗性參數(shù)的確定提供參考。圖4(a)和4(b)分別給出了參考升力系數(shù)和模型中的經(jīng)驗性參數(shù)隨間距比的變化規(guī)律。

    在T/D=3.0 工況下,Liu 等得到了彈性并列雙圓柱的橫向位移(Y/D)響應(yīng)隨時間的變化,如圖5(a)。結(jié)果顯示,兩個圓柱的位移響應(yīng)均為一般的簡諧運動曲線,而且兩者之間的振動互為反相位。根據(jù)已確定的模型參數(shù),對本文模型的方程組(6a)~(6d)進行數(shù)值求解,圖5(b)給出了經(jīng)驗?zāi)P偷念A(yù)報結(jié)果,兩個圓柱的位移響應(yīng)均為簡諧曲線,而且相位相差半個周期,即互為反相位。對比Liu 等[11]的實驗值,可以看到經(jīng)驗?zāi)P皖A(yù)測的振動頻率與其吻合很好,但振動幅值偏小。

    圖4 參考升力系數(shù)和經(jīng)驗性參數(shù)隨間距比的變化規(guī)律Fig.4 The reference lift coefficients and the empirical parameters versus spacing ratios

    在T/D=1.8 條件下,相應(yīng)的尾渦形為偏向流模式(biased flow pattern)。Liu 等的實驗得到的升力和阻力時序曲線表現(xiàn)為強烈的不規(guī)則性,他們將產(chǎn)生的原因歸結(jié)于尾渦的相互作用,而非轉(zhuǎn)換性質(zhì)(transition behavior)。同時,實驗中監(jiān)測出升力CL 譜存在兩個主脫落頻率。圖5(c)給出了Liu 數(shù)值實驗得到的圓柱位移響應(yīng),可以看到響應(yīng)曲線表現(xiàn)為不規(guī)則性。與此不同,經(jīng)驗?zāi)P皖A(yù)報的圓柱振動中存在“拍”,即表明流體振子的一個頻率與結(jié)構(gòu)的固有頻率相接近,如圖5(d)所示。因此,實驗中表現(xiàn)的強不規(guī)則性可能還包含了流-固之間的耦合作用。

    當間距比取為T/D =1.1 時,兩個圓柱之間的距離很近,通常圓柱的位移幅值隨質(zhì)量比的增加而減小,因此選用大質(zhì)量比來保證兩圓柱不會發(fā)生碰撞。Liu 等設(shè)m*=66,并給出了圓柱的位移響應(yīng)曲線,如圖5(e)。從圓柱剪切層脫落的發(fā)展旋渦存在很強的相互作用,同時還受到間隙流的影響,導(dǎo)致位移響應(yīng)中有另一個額外峰值,位移譜也顯示出存在兩個主頻率。圖5(f)給出了經(jīng)驗?zāi)P偷念A(yù)報結(jié)果,其得到的振動頻率和幅值都與實驗結(jié)果一致。

    2.3 低質(zhì)量比響應(yīng)分析

    低質(zhì)量比并列雙圓柱渦激振動的模型實驗分別對T/D=4.7,4.2,3.9,3.5,3.0 和2.6 共六種工況進行了研究和討論。圖6(a)和6(b)給出了在間距比T/D=2.6,約化速度ur=6.3 條件下圓柱1 和圓柱2 借助加速度傳感器得到的結(jié)構(gòu)頂部橫向位移加速度結(jié)果,由于提取的信號中存在高頻噪聲,需要對其進行光順濾波處理,本文采用低通與高通濾波方式對原始數(shù)據(jù)進行處理[12]。濾波后的加速度結(jié)果通過數(shù)值積分可以得到橫向位移隨時間的變化曲線,如圖6(c)和6(d)所示。經(jīng)過譜分析從而得到了響應(yīng)的振幅和頻率,如圖6(e),以及斯特勞哈爾數(shù)隨間距比的變化關(guān)系,如圖6(f)所示。

    低質(zhì)量比渦激振動的顯著特點是圓柱的位移響應(yīng)幅值相對較大,因此有必要對圓柱體橫向無量綱最大振幅(Y/D)進行討論。數(shù)值求解經(jīng)驗?zāi)P偷姆匠探M,式(6a)~(6d),其中參考升力系數(shù)取為經(jīng)驗值CL01=CL02=0.3,模型的其它經(jīng)驗性參數(shù)的取值及其隨間距比的變化規(guī)律在圖6(g)中給出,對比中高質(zhì)量比的經(jīng)驗性參數(shù)(圖4),可以看到兩者的變化趨勢相近。

    圖6(h)給出了T/D=4.7 和單圓柱條件下,模型實驗得到的圓柱最大振幅與約化速度的變化關(guān)系??梢钥吹讲⒘须p圓柱的振幅響應(yīng)及峰值位置與單圓柱的測量結(jié)果接近,因此可以認為兩個圓柱之間的距離較遠,相互之間的影響較小。經(jīng)驗?zāi)P偷念A(yù)報結(jié)果與實驗結(jié)果的趨勢一致,振幅峰值都出現(xiàn)在約化速度ur=6.3附近,隨后最大振幅逐漸減小,但是預(yù)測的峰值數(shù)值偏小。

    圖6(i)給出了T/D=2.6 條件下,模型實驗得到的圓柱最大振幅與約化速度的關(guān)系。實驗中發(fā)現(xiàn)圓柱的振幅響應(yīng)在隨約化速度變化的過程中出現(xiàn)了兩個峰值,分別是在約化速度ur=6.3 和12.6 附近,其產(chǎn)生的原因可能是由于該間距比條件下結(jié)構(gòu)所受的流體力存在兩個主頻率,隨著流速的增加兩個主頻率分別激起結(jié)構(gòu)的共振從而出現(xiàn)兩個峰值。經(jīng)驗?zāi)P鸵差A(yù)報出兩個振動峰值,出現(xiàn)的位置和數(shù)值均小于實驗的觀測結(jié)果,但是最大振幅隨約化速度的變化趨勢與實驗結(jié)果基本一致。

    圖5 不同間距比條件下數(shù)值模擬和經(jīng)驗性模型得到的Y/D 隨時間的變化Fig.5 Y/D time series for numerical simulation and empirical model under different spacing ratios

    3 結(jié) 語

    本文主要針對并列雙圓柱的渦激振動,建立了一種新的尾流振子模型,用來預(yù)測結(jié)構(gòu)的渦激振動響應(yīng)?,F(xiàn)階段的研究表明,在特定間距比條件下,旋渦的脫落頻率(或斯特勞哈爾頻率)出現(xiàn)分支,即可能存在兩個主脫落頻率。因此,本文提出以兩個具有不同固有頻率的尾流振子來共同描述結(jié)構(gòu)的近壁尾渦動力特性,同時兩個流體振子均滿足van der Pol 方程,進一步推導(dǎo)出結(jié)構(gòu)振子和流體振子耦合的方程組,用來模擬并列雙圓柱的流固耦合現(xiàn)象。

    本文同時給出了經(jīng)驗?zāi)P偷臄?shù)值求解方法,并結(jié)合Liu 等[11]的中高質(zhì)量比數(shù)值模擬實驗和低質(zhì)量比模型實驗提供的設(shè)計數(shù)據(jù),對該經(jīng)驗?zāi)P瓦M行了數(shù)值計算,得到的預(yù)測結(jié)果與實驗結(jié)果進行了對比。分析結(jié)果表明,預(yù)測的結(jié)構(gòu)位移響應(yīng)和最大振動幅值等變化規(guī)律與相應(yīng)的實驗結(jié)果趨勢一致,能夠反映并列雙圓柱渦激振動的結(jié)構(gòu)振動特性和一般規(guī)律。

    [1]Zdravkovich M M.The effects of interference between circular cylinders in cross flow[J].Journal of Fluids and Structures.1987,1(2):239-261.

    [2]姚熊亮,陳起富.均勻流場中并列圓柱渦激振動的現(xiàn)狀與預(yù)測[J].海洋工程,1995,13(3):37-47.

    [3]Williamson C H K.Evolution of a single wake behind a pair of bluff bodies[J].Journal of Fluid Mechanics.1985,159:1-18.

    [4]Sumner D,Wong S S T,Price S J,et al.Fluid behavior of side-by-side circular cylinders in steady cross-flow[J].Journal of Fluids and Structures.1999,13:309-338.

    [5]Wang Z J,Zhou Y,Li H.Flow-visualization of a two side-by-side cylinder wake[J].Journal of Flow Visualization & Image Processing,2001,9:123-138.

    [6]Facchinetti M L,Langre E,F(xiàn)ontaine E,et al.VIV of two cylinders in tandem arrangement:analytical and numerical modeling[C]∥Proc. of the 12th (2002)International Offshore and Polar Engineering Conference.2002:524-531.

    [7]Facchinetti M L,Langre E,Biolley F.Coupling of structure and wake oscillators in vortex-induced vibrations[J].Journal of Fluids and Structures,2004,19:123-140.

    [8]Bishop R E D,Hassan A Y.The lift and drag forces on a circular cylinder in a flowing fluid[C]∥Proc.of Royal Society Series A,1963,277:32-50.

    [9]Hartlen R T,Currie I G.Lift-oscillator model of vortex-induced vibration[J].Journal of the Engineering Mechanics,1970,96:577-591.

    [10]Li Xiao-min,Guo Hai-yan,Meng Fan-shun.Nonlinear coupled in-line and cross-flow vortex-induced vibration analysis of top tensioned riser[J].China Ocean Engineering,2010,24(4):749-758.

    [11]Liu Y,So R M C,Lau Y L,et al.Numberical studies of two side-by-side elastic cylinders in a cross-flow[J].Journal of Fluids and Structures.2001,15:1009-1030.

    [12]張建僑,宋吉寧,呂林,等.柔性立管渦激振動實驗的數(shù)據(jù)分析[J].中國海洋平臺,2009,24(4):26-32.

    猜你喜歡
    渦激尾流振子
    不同間距比下串聯(lián)圓柱渦激振動數(shù)值模擬研究
    彈簧振子問題的分析與求解
    渦激振動發(fā)電裝置及其關(guān)鍵技術(shù)
    盤球立管結(jié)構(gòu)抑制渦激振動的數(shù)值分析方法研究
    電子制作(2018年14期)2018-08-21 01:38:42
    非線性Duffing擾動振子共振機制的研究
    飛機尾流的散射特性與探測技術(shù)綜述
    雷達學報(2017年6期)2017-03-26 07:53:06
    基于近似熵和混沌振子的電力諧波檢測與估計
    錐形流量計尾流流場分析
    柔性圓管在渦激振動下的模態(tài)響應(yīng)分析
    電磁彈簧振子實驗裝置的改進
    物理實驗(2015年10期)2015-02-28 17:36:55
    大片电影免费在线观看免费| 男女无遮挡免费网站观看| 免费在线观看完整版高清| 人成视频在线观看免费观看| 日韩中文字幕视频在线看片| 欧美人与性动交α欧美软件| 国产亚洲av片在线观看秒播厂| 99国产精品免费福利视频| 婷婷色麻豆天堂久久| 色吧在线观看| 国产免费视频播放在线视频| 九九爱精品视频在线观看| 国产一卡二卡三卡精品 | 久久影院123| 9热在线视频观看99| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成人av在线免费| 少妇被粗大猛烈的视频| 一二三四中文在线观看免费高清| 九九爱精品视频在线观看| 国产欧美亚洲国产| 久久久久久久精品精品| 成人漫画全彩无遮挡| 男的添女的下面高潮视频| 国产精品国产三级专区第一集| 久久久久精品久久久久真实原创| 免费久久久久久久精品成人欧美视频| 国产成人av激情在线播放| 亚洲综合精品二区| 热re99久久精品国产66热6| 新久久久久国产一级毛片| 国产精品二区激情视频| 亚洲欧美一区二区三区黑人| 久久这里只有精品19| 观看美女的网站| 国产黄色视频一区二区在线观看| tube8黄色片| 久久久国产精品麻豆| 成人毛片60女人毛片免费| 久久久久精品久久久久真实原创| 校园人妻丝袜中文字幕| 狂野欧美激情性bbbbbb| 午夜影院在线不卡| 国产av精品麻豆| av福利片在线| 精品久久久久久电影网| 制服人妻中文乱码| 久久久久网色| www.自偷自拍.com| 熟女av电影| 欧美在线一区亚洲| 日本vs欧美在线观看视频| 一边亲一边摸免费视频| 少妇猛男粗大的猛烈进出视频| 免费黄色在线免费观看| 国产野战对白在线观看| 看十八女毛片水多多多| 制服诱惑二区| 黄色视频在线播放观看不卡| 久久人人97超碰香蕉20202| 蜜桃在线观看..| 国产高清不卡午夜福利| 黄片无遮挡物在线观看| 在线观看三级黄色| av卡一久久| 建设人人有责人人尽责人人享有的| videos熟女内射| 国产午夜精品一二区理论片| 亚洲人成77777在线视频| 视频在线观看一区二区三区| 操美女的视频在线观看| 91精品伊人久久大香线蕉| 国产精品一二三区在线看| 亚洲伊人色综图| 久久天堂一区二区三区四区| 在线天堂最新版资源| 超色免费av| 一本色道久久久久久精品综合| 国产高清国产精品国产三级| 丁香六月欧美| 精品少妇一区二区三区视频日本电影 | 精品一区二区三区四区五区乱码 | 国产精品成人在线| 亚洲一级一片aⅴ在线观看| 制服人妻中文乱码| 一级毛片我不卡| 亚洲熟女毛片儿| 精品人妻一区二区三区麻豆| 免费在线观看完整版高清| 国产成人精品无人区| 国产一区二区在线观看av| 黄色怎么调成土黄色| 中文乱码字字幕精品一区二区三区| 嫩草影视91久久| 国产色婷婷99| 午夜福利免费观看在线| 大码成人一级视频| 久久人人爽人人片av| 人妻一区二区av| 成人黄色视频免费在线看| av在线播放精品| 久久久久久人妻| 最近中文字幕2019免费版| 少妇的丰满在线观看| 曰老女人黄片| 亚洲视频免费观看视频| 国产熟女午夜一区二区三区| 美国免费a级毛片| 国产又爽黄色视频| 国产极品天堂在线| 亚洲成人手机| 国产精品久久久久久精品古装| 亚洲欧美日韩另类电影网站| 免费在线观看黄色视频的| 青草久久国产| 人体艺术视频欧美日本| 亚洲激情五月婷婷啪啪| 亚洲欧美成人综合另类久久久| 又粗又硬又长又爽又黄的视频| av卡一久久| 亚洲av综合色区一区| 国产精品香港三级国产av潘金莲 | 男女无遮挡免费网站观看| 黑丝袜美女国产一区| 大香蕉久久成人网| www.精华液| 久久精品国产综合久久久| 国产精品蜜桃在线观看| 80岁老熟妇乱子伦牲交| 在线看a的网站| 又大又黄又爽视频免费| 久久久久网色| 国产欧美亚洲国产| 亚洲国产中文字幕在线视频| 亚洲精品成人av观看孕妇| 国产野战对白在线观看| 精品少妇黑人巨大在线播放| 亚洲精品aⅴ在线观看| 各种免费的搞黄视频| 黄色视频在线播放观看不卡| 成人免费观看视频高清| 亚洲国产中文字幕在线视频| 亚洲精品中文字幕在线视频| 制服丝袜香蕉在线| 黄片播放在线免费| tube8黄色片| 色精品久久人妻99蜜桃| 久久国产亚洲av麻豆专区| 男女午夜视频在线观看| 街头女战士在线观看网站| 丝袜美足系列| 亚洲视频免费观看视频| 日韩一卡2卡3卡4卡2021年| www日本在线高清视频| 午夜免费男女啪啪视频观看| 国产成人91sexporn| 久久久久精品性色| 精品酒店卫生间| 十分钟在线观看高清视频www| 美女国产高潮福利片在线看| 国产在线免费精品| 亚洲国产中文字幕在线视频| 老熟女久久久| 国产精品偷伦视频观看了| 成人18禁高潮啪啪吃奶动态图| 一边摸一边抽搐一进一出视频| 亚洲成人手机| 黄频高清免费视频| 少妇人妻 视频| 99久久综合免费| 亚洲图色成人| 老汉色av国产亚洲站长工具| 看非洲黑人一级黄片| 卡戴珊不雅视频在线播放| 在线亚洲精品国产二区图片欧美| 亚洲成色77777| 少妇人妻 视频| 亚洲婷婷狠狠爱综合网| 高清欧美精品videossex| 亚洲国产精品成人久久小说| 国产爽快片一区二区三区| 最近中文字幕2019免费版| 久久这里只有精品19| 国语对白做爰xxxⅹ性视频网站| 97在线人人人人妻| 久久女婷五月综合色啪小说| 亚洲欧美成人综合另类久久久| 嫩草影院入口| 老司机影院毛片| 免费观看人在逋| 尾随美女入室| h视频一区二区三区| 日韩人妻精品一区2区三区| 欧美国产精品一级二级三级| 一级a爱视频在线免费观看| 伊人久久国产一区二区| 男女午夜视频在线观看| 熟女av电影| 国产精品无大码| 天美传媒精品一区二区| 一级黄片播放器| 欧美黑人欧美精品刺激| 一级片免费观看大全| 最近最新中文字幕免费大全7| √禁漫天堂资源中文www| 欧美精品一区二区大全| 午夜老司机福利片| 国产av国产精品国产| 国产精品嫩草影院av在线观看| 国产一区有黄有色的免费视频| 日韩av不卡免费在线播放| 久久97久久精品| 国产精品久久久久成人av| 日日啪夜夜爽| 久久天躁狠狠躁夜夜2o2o | 大码成人一级视频| 精品少妇一区二区三区视频日本电影 | 另类亚洲欧美激情| 如日韩欧美国产精品一区二区三区| 午夜久久久在线观看| 日韩一区二区三区影片| 黄色视频不卡| 国产成人免费观看mmmm| 9热在线视频观看99| 侵犯人妻中文字幕一二三四区| 国产欧美日韩综合在线一区二区| 十分钟在线观看高清视频www| 欧美国产精品一级二级三级| 啦啦啦在线免费观看视频4| 免费黄网站久久成人精品| 国产成人精品在线电影| av片东京热男人的天堂| 丝袜脚勾引网站| 国产精品国产av在线观看| 国产精品无大码| 高清不卡的av网站| 国产免费又黄又爽又色| 国产成人精品久久二区二区91 | 国产亚洲一区二区精品| 十八禁网站网址无遮挡| 国产一区二区在线观看av| 97人妻天天添夜夜摸| 国产成人一区二区在线| 亚洲av欧美aⅴ国产| 一级毛片 在线播放| 一边摸一边抽搐一进一出视频| 黄片小视频在线播放| 日韩制服骚丝袜av| 免费日韩欧美在线观看| 深夜精品福利| 国产一区二区在线观看av| 国产免费一区二区三区四区乱码| 久久狼人影院| 深夜精品福利| 亚洲国产av新网站| 成年人午夜在线观看视频| 久久亚洲国产成人精品v| av片东京热男人的天堂| 丝袜脚勾引网站| 999久久久国产精品视频| 免费在线观看完整版高清| 亚洲欧洲国产日韩| 久久久久网色| 99热全是精品| 久久久国产精品麻豆| 免费不卡黄色视频| 王馨瑶露胸无遮挡在线观看| 自线自在国产av| 51午夜福利影视在线观看| 欧美另类一区| 亚洲国产成人一精品久久久| 美女脱内裤让男人舔精品视频| 精品一区在线观看国产| 亚洲在久久综合| 亚洲国产av新网站| 最近2019中文字幕mv第一页| 精品人妻熟女毛片av久久网站| 十分钟在线观看高清视频www| 五月天丁香电影| 一级爰片在线观看| 亚洲av中文av极速乱| 免费观看性生交大片5| 国产亚洲av片在线观看秒播厂| 永久免费av网站大全| 爱豆传媒免费全集在线观看| 我要看黄色一级片免费的| avwww免费| 最近最新中文字幕大全免费视频 | 一区二区av电影网| 91精品伊人久久大香线蕉| 国产精品久久久av美女十八| 制服诱惑二区| 亚洲国产精品成人久久小说| 性高湖久久久久久久久免费观看| 日韩中文字幕视频在线看片| 妹子高潮喷水视频| 国产男女超爽视频在线观看| 精品第一国产精品| 欧美av亚洲av综合av国产av | 亚洲国产欧美日韩在线播放| 国产女主播在线喷水免费视频网站| 亚洲 欧美一区二区三区| 男女无遮挡免费网站观看| 久久97久久精品| 少妇被粗大的猛进出69影院| 免费高清在线观看视频在线观看| 亚洲自偷自拍图片 自拍| 国产在线视频一区二区| 一级毛片 在线播放| 美女视频免费永久观看网站| 国产成人精品久久二区二区91 | 天天影视国产精品| 国产精品一国产av| 亚洲精品国产av成人精品| 久久天堂一区二区三区四区| 久久久精品国产亚洲av高清涩受| 97人妻天天添夜夜摸| 亚洲av日韩精品久久久久久密 | 国产一区二区在线观看av| 水蜜桃什么品种好| 国产精品.久久久| av网站在线播放免费| 亚洲国产欧美日韩在线播放| 大陆偷拍与自拍| 国产亚洲午夜精品一区二区久久| 男人舔女人的私密视频| 母亲3免费完整高清在线观看| 天堂8中文在线网| tube8黄色片| 国产精品无大码| 国产片特级美女逼逼视频| 啦啦啦 在线观看视频| 精品一区二区免费观看| 一级黄片播放器| 亚洲人成77777在线视频| 成人漫画全彩无遮挡| 国产日韩欧美亚洲二区| 国产黄色免费在线视频| 毛片一级片免费看久久久久| 久热这里只有精品99| 欧美人与性动交α欧美精品济南到| 亚洲欧美清纯卡通| 亚洲在久久综合| 午夜日本视频在线| 国产精品成人在线| 丝袜美足系列| 97在线人人人人妻| 波多野结衣一区麻豆| 国产片特级美女逼逼视频| 99re6热这里在线精品视频| 国产精品三级大全| 午夜福利视频在线观看免费| 日本午夜av视频| 水蜜桃什么品种好| 美女高潮到喷水免费观看| 国产精品免费大片| 嫩草影视91久久| 午夜av观看不卡| 国产精品久久久久久精品电影小说| 欧美日韩综合久久久久久| 国产精品一区二区在线不卡| 精品一品国产午夜福利视频| 精品国产一区二区久久| 丝袜脚勾引网站| av.在线天堂| 日本猛色少妇xxxxx猛交久久| 午夜福利视频在线观看免费| 国产欧美日韩综合在线一区二区| 国产精品免费视频内射| 成年美女黄网站色视频大全免费| 国产99久久九九免费精品| 久久久欧美国产精品| 狂野欧美激情性xxxx| 人人妻,人人澡人人爽秒播 | 国产又色又爽无遮挡免| 成人国产麻豆网| 免费在线观看视频国产中文字幕亚洲 | 国产欧美日韩综合在线一区二区| av不卡在线播放| 欧美黑人精品巨大| 超色免费av| 日日爽夜夜爽网站| 国产深夜福利视频在线观看| 一区二区三区精品91| av国产久精品久网站免费入址| 91成人精品电影| 久久99一区二区三区| 国产精品欧美亚洲77777| 欧美成人午夜精品| 国产欧美亚洲国产| 国产麻豆69| 婷婷色综合www| 国产日韩欧美在线精品| 国产有黄有色有爽视频| 欧美另类一区| 伊人久久国产一区二区| 免费观看a级毛片全部| 在线观看www视频免费| 不卡视频在线观看欧美| 黑人猛操日本美女一级片| 女人高潮潮喷娇喘18禁视频| 人体艺术视频欧美日本| 久久热在线av| 少妇精品久久久久久久| 国产女主播在线喷水免费视频网站| 一级毛片电影观看| 麻豆精品久久久久久蜜桃| 岛国毛片在线播放| 亚洲国产欧美在线一区| 亚洲精品一区蜜桃| 婷婷色综合大香蕉| 国产亚洲av高清不卡| 国产精品久久久久久人妻精品电影 | 精品久久蜜臀av无| av女优亚洲男人天堂| 1024香蕉在线观看| av免费观看日本| 最近2019中文字幕mv第一页| 一级毛片黄色毛片免费观看视频| 伦理电影免费视频| 欧美黑人精品巨大| av又黄又爽大尺度在线免费看| 女人高潮潮喷娇喘18禁视频| bbb黄色大片| 悠悠久久av| 国产一区二区三区av在线| 国产精品嫩草影院av在线观看| 国产免费现黄频在线看| 精品第一国产精品| 精品亚洲成国产av| 日本av免费视频播放| 国产精品一区二区精品视频观看| 不卡视频在线观看欧美| 免费av中文字幕在线| 狠狠精品人妻久久久久久综合| 亚洲欧美精品自产自拍| 男女高潮啪啪啪动态图| 精品人妻一区二区三区麻豆| 亚洲熟女精品中文字幕| 青春草国产在线视频| 欧美在线一区亚洲| 亚洲自偷自拍图片 自拍| 精品国产超薄肉色丝袜足j| 19禁男女啪啪无遮挡网站| 色吧在线观看| 麻豆精品久久久久久蜜桃| 天堂中文最新版在线下载| 国产精品女同一区二区软件| 国产在线一区二区三区精| 天天添夜夜摸| 各种免费的搞黄视频| 亚洲精品国产av蜜桃| 汤姆久久久久久久影院中文字幕| 国产99久久九九免费精品| 午夜福利乱码中文字幕| 国产一区二区在线观看av| 欧美日韩av久久| 成人亚洲欧美一区二区av| 午夜免费鲁丝| 另类亚洲欧美激情| 久久久国产一区二区| 日韩制服丝袜自拍偷拍| 尾随美女入室| 亚洲精品第二区| 老司机影院毛片| 中文字幕av电影在线播放| 亚洲人成网站在线观看播放| 狂野欧美激情性xxxx| 热re99久久精品国产66热6| 国产av码专区亚洲av| 丰满乱子伦码专区| 18禁观看日本| 国产成人精品久久二区二区91 | 青春草亚洲视频在线观看| 久久精品国产a三级三级三级| 亚洲天堂av无毛| 无限看片的www在线观看| 精品国产一区二区三区久久久樱花| 国精品久久久久久国模美| 男女无遮挡免费网站观看| 看非洲黑人一级黄片| 欧美成人午夜精品| 久久人妻熟女aⅴ| 国产成人精品在线电影| 色婷婷av一区二区三区视频| av线在线观看网站| 色网站视频免费| 亚洲精品日韩在线中文字幕| 亚洲综合精品二区| 叶爱在线成人免费视频播放| 午夜日韩欧美国产| 人人妻人人澡人人爽人人夜夜| 99久久99久久久精品蜜桃| 精品亚洲成a人片在线观看| 欧美变态另类bdsm刘玥| 少妇人妻久久综合中文| 三上悠亚av全集在线观看| 高清黄色对白视频在线免费看| 2021少妇久久久久久久久久久| 国产毛片在线视频| xxx大片免费视频| 嫩草影视91久久| 久久性视频一级片| 男男h啪啪无遮挡| 国产成人欧美| a级毛片在线看网站| 久久久国产一区二区| 日韩av不卡免费在线播放| 永久免费av网站大全| 久久热在线av| 久久久久精品性色| avwww免费| 满18在线观看网站| 精品一区在线观看国产| 韩国高清视频一区二区三区| 香蕉国产在线看| 久久精品国产综合久久久| 巨乳人妻的诱惑在线观看| 成人黄色视频免费在线看| 亚洲欧洲精品一区二区精品久久久 | av免费观看日本| av女优亚洲男人天堂| 九草在线视频观看| www.精华液| 麻豆精品久久久久久蜜桃| 免费女性裸体啪啪无遮挡网站| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美日韩另类电影网站| 女性生殖器流出的白浆| 久久久久人妻精品一区果冻| 久久久国产精品麻豆| 夜夜骑夜夜射夜夜干| 午夜福利网站1000一区二区三区| 婷婷色av中文字幕| 一级毛片 在线播放| 亚洲人成网站在线观看播放| 啦啦啦中文免费视频观看日本| 中文字幕人妻丝袜制服| 麻豆精品久久久久久蜜桃| 操美女的视频在线观看| 另类亚洲欧美激情| 日日爽夜夜爽网站| 成人毛片60女人毛片免费| 婷婷色av中文字幕| 亚洲综合色网址| 男女床上黄色一级片免费看| 亚洲综合精品二区| 亚洲人成电影观看| 高清在线视频一区二区三区| 伦理电影免费视频| 亚洲av电影在线进入| 日韩 亚洲 欧美在线| 久久精品久久久久久噜噜老黄| 美女午夜性视频免费| 久久综合国产亚洲精品| 亚洲中文av在线| 中文乱码字字幕精品一区二区三区| 香蕉国产在线看| 日本欧美视频一区| 天美传媒精品一区二区| 国产熟女欧美一区二区| 侵犯人妻中文字幕一二三四区| 少妇人妻 视频| 亚洲精品久久成人aⅴ小说| 亚洲,一卡二卡三卡| 国产成人精品无人区| 国产精品蜜桃在线观看| 午夜激情久久久久久久| 亚洲国产日韩一区二区| 男人爽女人下面视频在线观看| 日本午夜av视频| 欧美97在线视频| 看免费av毛片| 男女下面插进去视频免费观看| 如何舔出高潮| 美女扒开内裤让男人捅视频| 岛国毛片在线播放| 下体分泌物呈黄色| 成人国产麻豆网| 两个人免费观看高清视频| 91精品国产国语对白视频| 老鸭窝网址在线观看| 免费黄网站久久成人精品| svipshipincom国产片| 午夜福利免费观看在线| 久久久久国产精品人妻一区二区| 国产麻豆69| av视频免费观看在线观看| 亚洲婷婷狠狠爱综合网| 青草久久国产| 搡老岳熟女国产| 啦啦啦中文免费视频观看日本| 国产精品 国内视频| av线在线观看网站| 成人影院久久| 你懂的网址亚洲精品在线观看| 天堂俺去俺来也www色官网| 亚洲精品乱久久久久久| 大香蕉久久网| 久久久久国产精品人妻一区二区| 成人国产av品久久久| 日韩中文字幕视频在线看片| 精品少妇黑人巨大在线播放| 国产片特级美女逼逼视频| 精品一区二区三区四区五区乱码 | 女性生殖器流出的白浆| 国产亚洲av高清不卡| 18禁裸乳无遮挡动漫免费视频| 亚洲男人天堂网一区| 2021少妇久久久久久久久久久| 日本午夜av视频| 久久久欧美国产精品| 丝袜美足系列| 亚洲一区二区三区欧美精品| 精品亚洲成国产av|