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

    基于時(shí)間步長(zhǎng)自適應(yīng)技術(shù)的海上浮基風(fēng)電平臺(tái)繞流數(shù)值模擬分析

    2013-10-30 08:12:16王開(kāi)加程建生段金輝王景全
    船舶與海洋工程 2013年2期
    關(guān)鍵詞:升力步長(zhǎng)圓柱

    王開(kāi)加,程建生,段金輝,王景全,王 亮

    (解放軍理工大學(xué)野戰(zhàn)工程學(xué)院,南京 210007)

    0 引 言

    隨著能源問(wèn)題的日趨嚴(yán)峻,開(kāi)發(fā)海上風(fēng)能并向遠(yuǎn)海、深海發(fā)展是必然趨勢(shì),漂浮式結(jié)構(gòu)是深?;A(chǔ)的主要方案。而解決圓柱繞流問(wèn)題是分析深遠(yuǎn)海漂浮式風(fēng)電基礎(chǔ)平臺(tái)在波浪和海流的作用下的動(dòng)力特性的關(guān)鍵。對(duì)于多圓柱的研究一直是流體力學(xué)領(lǐng)域中一個(gè)重要研究課題。由于繞流場(chǎng)的位置和形狀的復(fù)雜性,數(shù)值求解成為十分重要的手段[1]。

    目前,對(duì)圓柱繞流采用了不同的方法。對(duì)不同雷諾數(shù)情況下的圓柱繞流問(wèn)題包括二維數(shù)值模擬作了研究。而在二維數(shù)值模擬研究中,易獲得較好的阻力系數(shù)DC 和斯特魯哈數(shù)St,但升力系數(shù)LC的數(shù)值計(jì)算結(jié)果與真實(shí)值相差較大,這是由于流場(chǎng)網(wǎng)格(特別是邊界層附近區(qū)域的流場(chǎng)網(wǎng)格)不夠細(xì)密及時(shí)間步長(zhǎng)的選取不合理造成的。在流固耦合的數(shù)值模擬中,網(wǎng)格設(shè)計(jì)及時(shí)間步長(zhǎng)是影響數(shù)值模擬精度的重要因素。本文在應(yīng)用自適應(yīng)時(shí)間步長(zhǎng)理論及小雷諾數(shù)(Re=100)情況下,采用有限體積法對(duì)圓柱繞流進(jìn)行了數(shù)值模擬,并對(duì)流場(chǎng)特性作了分析。

    1 數(shù)學(xué)模型

    1.1 概念化設(shè)計(jì)

    以REpower公司研制開(kāi)發(fā)的5MW的水平軸風(fēng)力發(fā)電機(jī)為設(shè)計(jì)依據(jù),自主設(shè)計(jì)了一種風(fēng)電浮式基礎(chǔ)平臺(tái)見(jiàn)圖1,結(jié)構(gòu)主體采用半潛式浮式平臺(tái),為了降低整個(gè)系統(tǒng)的重心,底部采用懸掛重物的辦法。概念設(shè)計(jì)的風(fēng)電浮基系統(tǒng)由3個(gè)部分組成:風(fēng)電機(jī)組和塔筒組成的風(fēng)機(jī)系統(tǒng),甲板、立柱和雙浮體組成的半潛式平臺(tái)和包括6根連接桿件在內(nèi)的壓載物部分。確定了初步設(shè)計(jì)尺寸參數(shù),并建立了幾何模型,見(jiàn)圖2。

    圖1 海上浮基風(fēng)電平臺(tái)概念化設(shè)計(jì)

    圖2 海上浮基風(fēng)電平臺(tái)模型

    1.2 參數(shù)的定義

    在計(jì)算流體力學(xué)中,一般用無(wú)量綱的升力系數(shù)LC和阻力系數(shù)DC 來(lái)衡量某物體的受力情況。升力系數(shù)和阻力系數(shù)定義如下:

    式中:DF和LF——分別為作用于單位長(zhǎng)度圓柱上的阻力和升力;U——入流速度。

    斯特魯哈數(shù)St是描述圓柱繞流的一個(gè)非常重要的無(wú)量綱量,表征旋渦脫離情況。其定義為:

    式中:T——圓柱單側(cè)旋渦脫落周期;D——圓柱直徑;U——未受干擾的自由來(lái)流速度。

    1.3 控制方程

    對(duì)于不可壓縮黏性流體,在笛卡爾坐標(biāo)系下,其運(yùn)動(dòng)的數(shù)值解受控于N-S方程,平面流動(dòng)的連續(xù)性方程和動(dòng)量方程為:

    式中:u、v——分別為x,y方向的速度分量,p——壓力,υ——流體的運(yùn)動(dòng)黏性系數(shù),ρ——流體密度。

    2 數(shù)值計(jì)算

    2.1 邊界條件

    流動(dòng)方向與垂直流動(dòng)方向上的均勻速度分別定義為u和v。

    進(jìn)口處:給定速度進(jìn)口,u=U,v=0;

    出口處:在計(jì)算域的出口處,設(shè)置壓力出口邊界條件;

    上下邊界:采用對(duì)稱邊界條件;

    圓柱表面:給定無(wú)滑移、無(wú)穿透邊界條件,u=U,v=0。

    2.2 離散格式及求解

    計(jì)算雷諾數(shù)(Re=100)為小雷諾數(shù),屬層流范圍,故采用Laminar模型。采用基于壓力基的分離式求解器進(jìn)行求解,計(jì)算中采用具有二階隱式時(shí)間格式的非定常流動(dòng)進(jìn)行計(jì)算。壓力項(xiàng)與速度項(xiàng)的耦合項(xiàng)計(jì)算采用PISO算法實(shí)現(xiàn),壓力項(xiàng)離散采用具有二階精度的格式離散,動(dòng)量方程采用二階迎風(fēng)格式離散。計(jì)算中壓力、密度、體積力和動(dòng)量項(xiàng)的欠松弛因子分別為0.3、1、1和0.7。

    2.3 幾何模型

    結(jié)合本文研究目標(biāo),取圓柱直徑 2mD= ,計(jì)算區(qū)域?yàn)?030DD×的矩形區(qū)域,如圖3所示。上游尺寸20D,下游尺寸30D,D為物體垂直于來(lái)流方向平面上的特征尺寸,對(duì)圓柱一般取直徑。

    圖3 計(jì)算模型

    圖4 網(wǎng)格劃分

    3 網(wǎng)格和計(jì)算時(shí)間步長(zhǎng)

    3.1 網(wǎng)格劃分

    3.2 自適應(yīng)時(shí)間步長(zhǎng)

    在流固耦合的數(shù)值模擬中,時(shí)間步長(zhǎng)是影響數(shù)值模擬精度的重要因素。如果時(shí)間步長(zhǎng)過(guò)大,則圓柱和流場(chǎng)的耦合計(jì)算不夠緊密,會(huì)導(dǎo)致計(jì)算結(jié)果精度不高,甚至?xí)霈F(xiàn)錯(cuò)誤的結(jié)果。如果時(shí)間步長(zhǎng)太小,則無(wú)謂地增加計(jì)算量。目前有關(guān)文獻(xiàn)資料中對(duì)時(shí)間步長(zhǎng)的選擇方案有兩種:一種是取圓柱自振周期的0.2%[2],另外一種是時(shí)間步長(zhǎng)取遠(yuǎn)處來(lái)流速度經(jīng)過(guò)圓柱直徑長(zhǎng)度所耗費(fèi)時(shí)間(對(duì)流時(shí)間單位)的0.2%[3]。

    在Fluent中,當(dāng)求解器采用PISO時(shí),即使時(shí)間步長(zhǎng)取得較大,計(jì)算仍然會(huì)收斂,但結(jié)果的準(zhǔn)確性無(wú)法得到保證。關(guān)于時(shí)間步長(zhǎng)的選取,引入CFL穩(wěn)定性限制條件,即

    式中:cV——網(wǎng)格單元的體積(二維時(shí)取為面積);d——維數(shù)(二維計(jì)算取為2,三維計(jì)算取為3);——相應(yīng)網(wǎng)格單元內(nèi)流速。式(6)表示計(jì)算的時(shí)間步長(zhǎng)tΔ近似等于流速?gòu)囊粋€(gè)網(wǎng)格點(diǎn)傳播到另一個(gè)網(wǎng)格點(diǎn)所需的時(shí)間。由于采用的是非均勻網(wǎng)格以及每個(gè)網(wǎng)格點(diǎn)處當(dāng)?shù)氐牧魉俨灰粯?,所以每個(gè)網(wǎng)格點(diǎn)處的tΔ是不一樣的。因此,在給定的時(shí)刻t和給定的網(wǎng)格點(diǎn)i上,式(6)應(yīng)寫成:

    為了保證流場(chǎng)在推進(jìn)求解中不發(fā)生“扭曲”現(xiàn)象,選擇在所有網(wǎng)格點(diǎn)(總數(shù)設(shè)為M)上計(jì)算(i=1,2,…,M),從中選擇最小的一個(gè),再附一調(diào)節(jié)系數(shù)α,即?。?/p>

    式中:α——自適應(yīng)時(shí)間步長(zhǎng)適應(yīng)系數(shù)。

    采用上述的自適應(yīng)時(shí)間步長(zhǎng)方案,通過(guò)用戶自定義功能UDF的二次開(kāi)發(fā)接口對(duì)Fluent軟件二次開(kāi)發(fā)來(lái)實(shí)現(xiàn)自適應(yīng)時(shí)間步長(zhǎng)方案在圓柱周圍流場(chǎng)的數(shù)值模擬中應(yīng)用。

    4 計(jì)算結(jié)果及分析

    4.1 模擬工況

    采用同一條件進(jìn)行繞流計(jì)算:入流速度為1m/s,圓柱直徑D=2m,Re=100。對(duì)自適應(yīng)時(shí)間步長(zhǎng)適應(yīng)系數(shù)α為0.5、1.0、1.5和固定時(shí)間步長(zhǎng)為0.025、0.050、0.100共6種工況進(jìn)行了數(shù)值模擬。數(shù)值模擬結(jié)果見(jiàn)圖5~7,同其他數(shù)值模擬結(jié)果和實(shí)驗(yàn)結(jié)果的比較見(jiàn)表1、2。

    4.2 結(jié)果分析

    4.2.1 時(shí)間步長(zhǎng)的影響

    針對(duì)上述6種工況進(jìn)行數(shù)值模擬,給出了采用自適應(yīng)時(shí)間步長(zhǎng)時(shí)α對(duì)計(jì)算結(jié)果的影響,如表1所示,通過(guò)與前人研究結(jié)果比較知,該計(jì)算結(jié)果是可靠的并提高了計(jì)算精度。從表1中可以看到,隨著α的增大,St數(shù)、升、阻力系數(shù)都會(huì)逐漸減小。當(dāng)α=1.0時(shí),St數(shù)、升、阻力系數(shù)相對(duì)于α=1.5時(shí)的數(shù)據(jù)變化不大,但實(shí)際運(yùn)算中計(jì)算量卻要大50%,因此全面考慮計(jì)算時(shí)間、效率以及經(jīng)驗(yàn),取α=1.5時(shí)自適應(yīng)時(shí)間步長(zhǎng)算法為宜。采用固定時(shí)間步長(zhǎng)時(shí),St數(shù)、升、阻力系數(shù)隨時(shí)間步長(zhǎng)dt增大而減小的趨勢(shì)比較明顯。如當(dāng)dt= 0 .100,計(jì)算所得的 St數(shù)和阻力系數(shù)相比于 dt= 0 .025時(shí)減小得并不太多,但升力系數(shù)的振幅卻下降了約三分之一。

    圖5是Re=100,1.5α=情況下的時(shí)間步長(zhǎng)時(shí)程曲線。從圖中可以明顯地看出,1.5α=情況下采用自適應(yīng)時(shí)間步長(zhǎng)算法的時(shí)間步長(zhǎng)均大于0.025,并且計(jì)算結(jié)果優(yōu)于后者。

    表1 升阻力系數(shù)及斯特魯哈數(shù)對(duì)比表

    圖5 Re=100,α=1.5情況下的時(shí)間步長(zhǎng)時(shí)程曲線

    圖6給出了采用自適應(yīng)時(shí)間步長(zhǎng)(1.5α=)時(shí)圓柱升力系數(shù)和阻力系數(shù)的時(shí)程曲線。由于圓柱尾流形成了卡門渦街,當(dāng)周期性的渦脫處于圓柱上方時(shí)圓柱的升力最大,而后渦脫會(huì)逐漸向下游運(yùn)動(dòng),經(jīng)過(guò)圓柱中心時(shí)升力系數(shù)為0,當(dāng)渦脫運(yùn)動(dòng)到圓柱下游時(shí),升力系數(shù)達(dá)到最小負(fù)值,如此往復(fù),故升力系數(shù)的均值為0。阻力系數(shù)在渦脫達(dá)到穩(wěn)定后的均值為1.342,這與文獻(xiàn)[4]中Re為100的計(jì)算結(jié)果非常吻合。

    由圖6阻力系數(shù)與升力系數(shù)時(shí)程曲線可以看出,升力系數(shù)發(fā)生1個(gè)周期變化的同時(shí),阻力系數(shù)就會(huì)發(fā)生2個(gè)周期的變化,即阻力系數(shù)的頻率約為升力系數(shù)頻率的2倍。這是由于圓柱發(fā)生周期性渦脫時(shí),從上下表面脫離的渦會(huì)引起阻力改變1次,而這2個(gè)渦共同影響升力變化1次。

    圖6 Re=100,α=1.5情況下的圓柱阻力系數(shù)與升力系數(shù)時(shí)程曲線

    圖7給出了圓柱邊界層附近的矢量場(chǎng),從中可以看到在邊界層內(nèi),矢量沿圓柱外法線方向逐漸增大,準(zhǔn)確地模擬出速度矢量在邊界層內(nèi)的變化規(guī)律。

    圖7 圓柱邊界層附近矢量場(chǎng)

    4.2.2 網(wǎng)格設(shè)計(jì)的影響

    在靜止圓柱繞流的數(shù)值模擬中,比較容易獲得較好的阻力系數(shù)DC和斯特哈爾數(shù)tS,升力系數(shù)LC的數(shù)值計(jì)算結(jié)果有時(shí)同真實(shí)值相比則偏低。原因是因流場(chǎng)網(wǎng)格劃分不夠精細(xì),特別是邊界層網(wǎng)格劃分不夠精細(xì)造成的。在邊界層網(wǎng)格的設(shè)計(jì)中,除了沿圓周有足夠的網(wǎng)格數(shù),而且第一層網(wǎng)格到圓柱壁面的距離也要足夠小。表2給出了一個(gè)比較:

    表2 粗糙網(wǎng)格和精細(xì)網(wǎng)格對(duì)數(shù)值模擬結(jié)果的影響

    從表2的比較可以看出,采用精細(xì)網(wǎng)格比采用粗糙網(wǎng)格能獲得更好的數(shù)值模擬結(jié)果,特別是升力系數(shù)的幅值有了提高,更加接近精確值。因此要想獲得較好的數(shù)值模擬結(jié)果,流場(chǎng)計(jì)算區(qū)域要足夠大,流場(chǎng)網(wǎng)格劃分要足夠精細(xì),邊界層附近區(qū)域的流場(chǎng)網(wǎng)格除了沿圓周要有足夠多的網(wǎng)格數(shù),而且第一層網(wǎng)格到圓柱壁面的距離要足夠小。

    5 結(jié) 語(yǔ)

    深海風(fēng)電具有風(fēng)資源豐富、無(wú)噪音污染、風(fēng)電設(shè)備利用率高、不受水深和海底地質(zhì)條件限制等優(yōu)點(diǎn),是未來(lái)海上風(fēng)電發(fā)展的主要方向。風(fēng)電浮式基礎(chǔ),是一種專門為深海區(qū)域風(fēng)電開(kāi)發(fā)的新型海洋基礎(chǔ),但由于深海區(qū)域的各種復(fù)雜環(huán)境條件引起的作用力對(duì)基礎(chǔ)結(jié)構(gòu)的安全造成很大的影響,目前國(guó)內(nèi)外對(duì)于深海風(fēng)電浮基的研究尚處于起步階段。采用時(shí)間步長(zhǎng)自適應(yīng)技術(shù)和精細(xì)的邊界層網(wǎng)格處理技術(shù)能夠獲得較為精確的數(shù)值結(jié)果。

    在雷諾數(shù)較低、圓柱周圍的流動(dòng)主要呈現(xiàn)二維流動(dòng)的情況下,采用基于有限體積法的自適應(yīng)時(shí)間步長(zhǎng)算法,流場(chǎng)網(wǎng)格劃分較好、特別是邊界層網(wǎng)格劃分較精細(xì),完全可以用Fluent軟件準(zhǔn)確地模擬靜止圓柱繞流問(wèn)題。本文應(yīng)用Fluent軟件對(duì)單圓柱繞流流場(chǎng)進(jìn)行數(shù)值模擬,所采用的自適應(yīng)時(shí)間步長(zhǎng)算法在現(xiàn)有的計(jì)算設(shè)備下可以獲得較好的數(shù)值模擬結(jié)果,與固定時(shí)間步長(zhǎng)算法相比提高了數(shù)值計(jì)算的精度。為后續(xù)風(fēng)電浮式基礎(chǔ)結(jié)構(gòu)的優(yōu)化設(shè)計(jì)和將來(lái)深海風(fēng)電浮基安裝、現(xiàn)場(chǎng)作業(yè)等具有指導(dǎo)性的意義。

    [1] 王 健,李海濤. 計(jì)算流體力學(xué)方法在船舶領(lǐng)域的實(shí)用性研究[J]. 船舶與海洋工程,2012, (4): 6-11.

    [2] 曹豐產(chǎn),項(xiàng)海帆. 圓柱非定常繞流及渦致振動(dòng)的數(shù)值計(jì)算[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯),2001, 16(1): 111-118.

    [3] Newman D J, Karniadakis G E. A direct numerical simulation study of flow past a freely vibrating cable[J]. Journal of Fluid Mechanics, 1997, 344: 95-136.

    [4] Kim J, Kim D, Choi H. An immersed-boundary finite-volume method for simulations of flow in complex geometries [J]. J Comput Phys, 2001,171:132-150.

    [5] Braza M, Chassaing P, Ha Mind H. Numerical study and physical analysis of the pressure and velocity fields in the near wake of a circular cylinder [J]. J Fluid Mech, 1986, 165: 79-130.

    [6] Lu XY, Dalton C, Yhang JF. Application of large eddy simulation to an oscillating flow past a circular cylinder [J]. Journal of Fluids Engineering, Transactions of the ASME. 1997, 119(3): 519-525.

    [7] 蘇銘德,康欽軍. 亞臨界雷諾數(shù)下圓柱繞流的大渦模擬[J]. 力學(xué)學(xué)報(bào),1999, 31(1): 100-105.

    猜你喜歡
    升力步長(zhǎng)圓柱
    工程學(xué)和圓柱
    高速列車車頂–升力翼組合體氣動(dòng)特性
    圓柱的體積計(jì)算
    基于Armijo搜索步長(zhǎng)的BFGS與DFP擬牛頓法的比較研究
    無(wú)人機(jī)升力測(cè)試裝置設(shè)計(jì)及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    削法不同 體積有異
    升力式再入飛行器體襟翼姿態(tài)控制方法
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥搜索算法
    一種新型光伏系統(tǒng)MPPT變步長(zhǎng)滯環(huán)比較P&O法
    在线免费观看不下载黄p国产| 我要看黄色一级片免费的| 性高湖久久久久久久久免费观看| 天堂俺去俺来也www色官网| 久久精品国产综合久久久| 日本爱情动作片www.在线观看| 9色porny在线观看| 天天添夜夜摸| 晚上一个人看的免费电影| 日韩免费高清中文字幕av| 久久精品亚洲熟妇少妇任你| 狂野欧美激情性xxxx| 男女边吃奶边做爰视频| 黄色毛片三级朝国网站| 色婷婷久久久亚洲欧美| 婷婷色av中文字幕| 国产日韩欧美亚洲二区| kizo精华| 叶爱在线成人免费视频播放| 天天躁夜夜躁狠狠久久av| 日韩精品有码人妻一区| svipshipincom国产片| 亚洲精品av麻豆狂野| 高清欧美精品videossex| netflix在线观看网站| 精品少妇内射三级| 少妇的丰满在线观看| 少妇人妻 视频| 99热国产这里只有精品6| 又大又黄又爽视频免费| 日本av免费视频播放| 国产黄色免费在线视频| 中文字幕人妻熟女乱码| 可以免费在线观看a视频的电影网站 | 国产视频首页在线观看| 一区在线观看完整版| 日韩一本色道免费dvd| 成人亚洲精品一区在线观看| 亚洲av日韩精品久久久久久密 | 婷婷色麻豆天堂久久| 亚洲,欧美精品.| 日本91视频免费播放| 飞空精品影院首页| 亚洲熟女毛片儿| 狂野欧美激情性bbbbbb| 亚洲欧美激情在线| 制服人妻中文乱码| 国产97色在线日韩免费| 亚洲美女黄色视频免费看| 人成视频在线观看免费观看| 亚洲,欧美,日韩| 精品一品国产午夜福利视频| 少妇精品久久久久久久| 国产亚洲精品第一综合不卡| 伊人久久大香线蕉亚洲五| 丝袜在线中文字幕| 国产成人系列免费观看| 97人妻天天添夜夜摸| 亚洲综合精品二区| 国产在视频线精品| 老司机在亚洲福利影院| 九色亚洲精品在线播放| 韩国av在线不卡| 午夜福利影视在线免费观看| 91精品国产国语对白视频| 日韩中文字幕视频在线看片| 久久精品亚洲av国产电影网| 日韩,欧美,国产一区二区三区| 1024香蕉在线观看| 久久天堂一区二区三区四区| 久热这里只有精品99| 亚洲国产精品999| 巨乳人妻的诱惑在线观看| 国产精品久久久久久久久免| 丝袜脚勾引网站| 一本大道久久a久久精品| 久久 成人 亚洲| 欧美久久黑人一区二区| 高清av免费在线| 在线观看三级黄色| 亚洲av成人精品一二三区| 国产精品一国产av| 日本色播在线视频| av在线观看视频网站免费| 国产精品一国产av| 亚洲美女黄色视频免费看| 亚洲精品自拍成人| 在线看a的网站| 97在线人人人人妻| 亚洲成人一二三区av| 老司机靠b影院| 啦啦啦啦在线视频资源| 性色av一级| 国产精品蜜桃在线观看| 久久人人爽av亚洲精品天堂| av天堂久久9| 国产精品亚洲av一区麻豆 | 国产亚洲欧美精品永久| 婷婷色综合大香蕉| 免费观看a级毛片全部| 久久精品亚洲av国产电影网| 国产极品粉嫩免费观看在线| 国产成人欧美在线观看 | 国产精品三级大全| 亚洲国产av影院在线观看| 久久精品久久久久久噜噜老黄| 免费观看a级毛片全部| 久久久久国产一级毛片高清牌| 人人妻人人添人人爽欧美一区卜| 色婷婷久久久亚洲欧美| 黄网站色视频无遮挡免费观看| 亚洲精品中文字幕在线视频| 久久久久久人人人人人| 人人妻,人人澡人人爽秒播 | av国产精品久久久久影院| 夜夜骑夜夜射夜夜干| 国产免费视频播放在线视频| 一本大道久久a久久精品| 久久精品久久久久久久性| tube8黄色片| 亚洲第一青青草原| 亚洲精华国产精华液的使用体验| 欧美黑人欧美精品刺激| 19禁男女啪啪无遮挡网站| 免费观看人在逋| 纵有疾风起免费观看全集完整版| 久久99精品国语久久久| 久久这里只有精品19| 国产精品 欧美亚洲| 国产免费视频播放在线视频| 咕卡用的链子| 免费在线观看视频国产中文字幕亚洲 | 美女扒开内裤让男人捅视频| 日本vs欧美在线观看视频| 国产国语露脸激情在线看| 亚洲美女视频黄频| 精品国产乱码久久久久久小说| 波多野结衣一区麻豆| 香蕉国产在线看| 欧美日韩视频精品一区| 黄片无遮挡物在线观看| 黄片无遮挡物在线观看| 又大又爽又粗| 91国产中文字幕| 在线观看免费日韩欧美大片| 两性夫妻黄色片| 街头女战士在线观看网站| 性少妇av在线| 日日撸夜夜添| 一本色道久久久久久精品综合| 国产片内射在线| 高清不卡的av网站| 亚洲 欧美一区二区三区| 热99国产精品久久久久久7| 在线免费观看不下载黄p国产| 国产精品久久久久久精品古装| 国产精品久久久久久精品古装| 免费少妇av软件| 午夜免费男女啪啪视频观看| 男女床上黄色一级片免费看| 你懂的网址亚洲精品在线观看| 伦理电影大哥的女人| 最近中文字幕2019免费版| 国产有黄有色有爽视频| 亚洲一码二码三码区别大吗| 美女主播在线视频| 亚洲欧美成人综合另类久久久| 国产激情久久老熟女| av在线播放精品| 热re99久久国产66热| 校园人妻丝袜中文字幕| 亚洲国产欧美网| 丰满少妇做爰视频| 亚洲熟女毛片儿| av在线app专区| 青春草视频在线免费观看| 成年动漫av网址| 另类亚洲欧美激情| 老汉色av国产亚洲站长工具| 欧美日韩亚洲综合一区二区三区_| 交换朋友夫妻互换小说| 黄片无遮挡物在线观看| 精品午夜福利在线看| 国产男人的电影天堂91| 亚洲国产欧美网| 在线观看www视频免费| 国产av在哪里看| 97碰自拍视频| 久热这里只有精品99| 欧美午夜高清在线| 99在线人妻在线中文字幕| 婷婷六月久久综合丁香| 国产三级黄色录像| 国产伦人伦偷精品视频| 91国产中文字幕| 久久久国产成人精品二区| 最近最新中文字幕大全电影3 | av有码第一页| 大陆偷拍与自拍| 免费人成视频x8x8入口观看| 999精品在线视频| 欧美绝顶高潮抽搐喷水| 在线观看66精品国产| 亚洲 国产 在线| 久久久国产欧美日韩av| 欧美老熟妇乱子伦牲交| 天天添夜夜摸| 久久久久国内视频| 黄色片一级片一级黄色片| 日韩欧美在线二视频| 在线国产一区二区在线| 久久久久国产精品人妻aⅴ院| 国产一级毛片七仙女欲春2 | 成人av一区二区三区在线看| 亚洲国产精品sss在线观看| 久久亚洲精品不卡| 久久人妻熟女aⅴ| 免费不卡黄色视频| 国产成人av激情在线播放| 黑人操中国人逼视频| 精品一区二区三区视频在线观看免费| 国产又爽黄色视频| 日韩欧美一区视频在线观看| 一区二区日韩欧美中文字幕| 国产一区二区在线av高清观看| 叶爱在线成人免费视频播放| 国产精品 国内视频| 免费在线观看黄色视频的| 精品国内亚洲2022精品成人| 国产成+人综合+亚洲专区| 在线观看66精品国产| 免费看a级黄色片| 女人高潮潮喷娇喘18禁视频| 午夜福利高清视频| 午夜影院日韩av| 国产精品久久久av美女十八| 1024香蕉在线观看| 久久久水蜜桃国产精品网| 午夜精品久久久久久毛片777| 两个人看的免费小视频| 亚洲欧美精品综合一区二区三区| 麻豆久久精品国产亚洲av| 日本 av在线| 久久香蕉国产精品| 久久久久久久精品吃奶| 亚洲欧美激情综合另类| 一卡2卡三卡四卡精品乱码亚洲| 日日夜夜操网爽| 久久婷婷成人综合色麻豆| 亚洲视频免费观看视频| 日韩三级视频一区二区三区| 久久亚洲真实| 夜夜夜夜夜久久久久| 亚洲色图综合在线观看| 搡老熟女国产l中国老女人| 亚洲专区国产一区二区| 国产野战对白在线观看| 国产主播在线观看一区二区| av在线天堂中文字幕| 大香蕉久久成人网| 国内久久婷婷六月综合欲色啪| 欧美精品亚洲一区二区| 麻豆久久精品国产亚洲av| 日本 av在线| tocl精华| 亚洲成国产人片在线观看| 性欧美人与动物交配| 久久人妻av系列| 熟妇人妻久久中文字幕3abv| 国产蜜桃级精品一区二区三区| 精品国产美女av久久久久小说| 免费高清在线观看日韩| 久久精品成人免费网站| 可以在线观看的亚洲视频| 国产精品 欧美亚洲| 国产精品98久久久久久宅男小说| 美女免费视频网站| 精品免费久久久久久久清纯| 国产97色在线日韩免费| 国产成人啪精品午夜网站| 黄片播放在线免费| 亚洲专区国产一区二区| 天堂√8在线中文| 99久久久亚洲精品蜜臀av| 侵犯人妻中文字幕一二三四区| 一区二区三区高清视频在线| 日本五十路高清| 久久久久久大精品| 十分钟在线观看高清视频www| 怎么达到女性高潮| 在线观看免费视频日本深夜| 亚洲欧美激情在线| avwww免费| 精品免费久久久久久久清纯| 无人区码免费观看不卡| 最新美女视频免费是黄的| 搡老熟女国产l中国老女人| 村上凉子中文字幕在线| videosex国产| 亚洲欧美激情综合另类| 看片在线看免费视频| 亚洲一区二区三区色噜噜| 精品人妻在线不人妻| 亚洲av第一区精品v没综合| 日本撒尿小便嘘嘘汇集6| 热re99久久国产66热| 亚洲三区欧美一区| 亚洲成av片中文字幕在线观看| 色精品久久人妻99蜜桃| 日本三级黄在线观看| 国产亚洲精品一区二区www| 老司机在亚洲福利影院| 日韩免费av在线播放| 精品乱码久久久久久99久播| av网站免费在线观看视频| 欧美色欧美亚洲另类二区 | 日韩大码丰满熟妇| 国产成人欧美在线观看| 老熟妇乱子伦视频在线观看| 村上凉子中文字幕在线| 黄色视频不卡| 国产成人av教育| 黄片播放在线免费| 中文字幕精品免费在线观看视频| 国产三级在线视频| 女性生殖器流出的白浆| 最新在线观看一区二区三区| 亚洲一区二区三区不卡视频| 免费少妇av软件| 久久人妻熟女aⅴ| 伦理电影免费视频| 国产一区二区三区视频了| 国产一区二区在线av高清观看| 久久性视频一级片| 夜夜躁狠狠躁天天躁| 亚洲 欧美 日韩 在线 免费| 99国产精品免费福利视频| 99久久国产精品久久久| 亚洲av五月六月丁香网| 国产精品二区激情视频| 中文字幕色久视频| 欧美精品亚洲一区二区| 亚洲少妇的诱惑av| 国产三级黄色录像| 亚洲一区高清亚洲精品| 欧美精品啪啪一区二区三区| 精品人妻1区二区| 国产精品日韩av在线免费观看 | 午夜免费鲁丝| 亚洲欧美精品综合一区二区三区| 国产一卡二卡三卡精品| x7x7x7水蜜桃| 欧美 亚洲 国产 日韩一| 18禁观看日本| 一二三四社区在线视频社区8| 18禁裸乳无遮挡免费网站照片 | 91国产中文字幕| 黄色丝袜av网址大全| 91精品国产国语对白视频| av电影中文网址| av中文乱码字幕在线| 欧美在线一区亚洲| 欧美最黄视频在线播放免费| 无限看片的www在线观看| 精品人妻在线不人妻| 欧美日韩一级在线毛片| 精品久久久久久久人妻蜜臀av | 男女午夜视频在线观看| 99国产精品一区二区蜜桃av| 这个男人来自地球电影免费观看| 色综合婷婷激情| 18禁观看日本| 亚洲国产精品成人综合色| 啦啦啦观看免费观看视频高清 | 免费在线观看黄色视频的| 国产精品1区2区在线观看.| 亚洲无线在线观看| 精品午夜福利视频在线观看一区| 国产欧美日韩一区二区三区在线| 亚洲国产欧美一区二区综合| 美国免费a级毛片| 国产又色又爽无遮挡免费看| 18禁裸乳无遮挡免费网站照片 | 日本 欧美在线| www.熟女人妻精品国产| 欧美日韩一级在线毛片| 看黄色毛片网站| 亚洲精品在线观看二区| 国产成年人精品一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 欧美一区二区三区| 老熟妇仑乱视频hdxx| 亚洲精品粉嫩美女一区| 亚洲一码二码三码区别大吗| 亚洲男人的天堂狠狠| 午夜两性在线视频| 亚洲人成77777在线视频| 精品国产一区二区久久| 亚洲 欧美一区二区三区| 国产精品影院久久| 日本 欧美在线| 两个人看的免费小视频| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 日韩大码丰满熟妇| 黑人操中国人逼视频| 日本欧美视频一区| 在线观看午夜福利视频| 脱女人内裤的视频| 日本五十路高清| 日本 欧美在线| 伦理电影免费视频| 每晚都被弄得嗷嗷叫到高潮| 久久久水蜜桃国产精品网| 美女 人体艺术 gogo| 韩国av一区二区三区四区| 久久狼人影院| 亚洲五月婷婷丁香| 在线观看舔阴道视频| 日韩精品免费视频一区二区三区| 久久精品成人免费网站| 国产日韩一区二区三区精品不卡| 妹子高潮喷水视频| 久久久久久国产a免费观看| 国产亚洲欧美98| 日韩精品中文字幕看吧| 欧美老熟妇乱子伦牲交| 国内精品久久久久久久电影| 精品国产国语对白av| av在线播放免费不卡| 国产色视频综合| 搡老熟女国产l中国老女人| 亚洲视频免费观看视频| 69精品国产乱码久久久| 美国免费a级毛片| 国产精品,欧美在线| 成年女人毛片免费观看观看9| av电影中文网址| 国产精品香港三级国产av潘金莲| 午夜免费激情av| 日日摸夜夜添夜夜添小说| 日韩欧美一区二区三区在线观看| 老熟妇仑乱视频hdxx| 亚洲精品美女久久av网站| 性少妇av在线| 97人妻天天添夜夜摸| 欧美日本中文国产一区发布| 可以免费在线观看a视频的电影网站| 亚洲av五月六月丁香网| 免费在线观看完整版高清| 三级毛片av免费| 9热在线视频观看99| www.999成人在线观看| bbb黄色大片| 久久精品国产综合久久久| 午夜福利一区二区在线看| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 88av欧美| 午夜成年电影在线免费观看| 韩国精品一区二区三区| 午夜福利视频1000在线观看 | 91字幕亚洲| 国产成人精品无人区| 两个人看的免费小视频| 久久欧美精品欧美久久欧美| 国产成人影院久久av| 麻豆一二三区av精品| 日韩欧美一区二区三区在线观看| 日韩精品免费视频一区二区三区| 亚洲人成77777在线视频| 成熟少妇高潮喷水视频| 身体一侧抽搐| 搡老妇女老女人老熟妇| 久久国产精品人妻蜜桃| 韩国精品一区二区三区| 一级,二级,三级黄色视频| 欧美大码av| 精品人妻1区二区| 亚洲精品久久国产高清桃花| 成人精品一区二区免费| 国产免费av片在线观看野外av| 在线观看日韩欧美| 色老头精品视频在线观看| 亚洲av美国av| 久久草成人影院| 久久婷婷成人综合色麻豆| 搡老妇女老女人老熟妇| 50天的宝宝边吃奶边哭怎么回事| 美女免费视频网站| 99精品久久久久人妻精品| 曰老女人黄片| 国产黄a三级三级三级人| 精品久久久久久成人av| 午夜福利高清视频| 电影成人av| 看黄色毛片网站| 久久中文看片网| 国产单亲对白刺激| 老汉色∧v一级毛片| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久久午夜电影| 中文字幕av电影在线播放| 亚洲激情在线av| 在线播放国产精品三级| avwww免费| 在线av久久热| 久久久久亚洲av毛片大全| 亚洲五月天丁香| 国产亚洲精品综合一区在线观看 | 午夜久久久在线观看| 午夜福利高清视频| 中国美女看黄片| 少妇熟女aⅴ在线视频| 国产一区二区三区综合在线观看| 级片在线观看| 久久中文看片网| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久亚洲av鲁大| 欧美日韩福利视频一区二区| 在线观看66精品国产| 国产高清有码在线观看视频 | 最好的美女福利视频网| 99在线视频只有这里精品首页| 老熟妇乱子伦视频在线观看| 叶爱在线成人免费视频播放| 真人做人爱边吃奶动态| 色综合亚洲欧美另类图片| 可以免费在线观看a视频的电影网站| 一进一出抽搐gif免费好疼| 日日干狠狠操夜夜爽| 真人做人爱边吃奶动态| 国产黄a三级三级三级人| 在线国产一区二区在线| 欧美黑人欧美精品刺激| 午夜久久久在线观看| 午夜福利18| 可以免费在线观看a视频的电影网站| 好男人电影高清在线观看| 桃红色精品国产亚洲av| 手机成人av网站| 久久精品国产亚洲av高清一级| 国产精品乱码一区二三区的特点 | 精品卡一卡二卡四卡免费| 色综合站精品国产| 99国产精品99久久久久| 99在线人妻在线中文字幕| 午夜久久久久精精品| av在线播放免费不卡| 中出人妻视频一区二区| 免费看a级黄色片| 在线观看免费午夜福利视频| 身体一侧抽搐| a在线观看视频网站| 精品福利观看| 久久人妻熟女aⅴ| 精品久久久精品久久久| 老汉色av国产亚洲站长工具| 国产片内射在线| 国产主播在线观看一区二区| 欧美黄色淫秽网站| 亚洲av成人不卡在线观看播放网| 男人操女人黄网站| 身体一侧抽搐| 亚洲五月婷婷丁香| 又紧又爽又黄一区二区| 日本 欧美在线| 欧美不卡视频在线免费观看 | 两性午夜刺激爽爽歪歪视频在线观看 | 激情视频va一区二区三区| 久久精品91无色码中文字幕| 欧美性长视频在线观看| 亚洲七黄色美女视频| 亚洲一区二区三区不卡视频| 成人亚洲精品一区在线观看| 亚洲中文日韩欧美视频| 欧美日本中文国产一区发布| 久久人人爽av亚洲精品天堂| 99热只有精品国产| 久久影院123| 国产91精品成人一区二区三区| 日韩欧美在线二视频| 99re在线观看精品视频| 久久草成人影院| 色在线成人网| 精品久久蜜臀av无| 亚洲成人精品中文字幕电影| 亚洲成a人片在线一区二区| 欧美日本亚洲视频在线播放| 日本撒尿小便嘘嘘汇集6| www.自偷自拍.com| 美女 人体艺术 gogo| 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 此物有八面人人有两片| e午夜精品久久久久久久| 亚洲精品一区av在线观看| 欧美中文日本在线观看视频| 韩国精品一区二区三区| 三级毛片av免费| 欧美日本亚洲视频在线播放| 国产区一区二久久| 亚洲国产精品sss在线观看| 亚洲av电影不卡..在线观看| 免费少妇av软件| 校园春色视频在线观看| 精品人妻在线不人妻| 欧美大码av| 91成年电影在线观看| 欧美激情久久久久久爽电影 | 男男h啪啪无遮挡| 淫秽高清视频在线观看| 黑人操中国人逼视频| 日本免费a在线| www日本在线高清视频|