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

    潛艇阻力及流場(chǎng)數(shù)值仿真策略優(yōu)化分析

    2022-05-07 06:16:50孫華偉郭春雨
    關(guān)鍵詞:模型

    李 鵬, 王 超, 孫華偉, 郭春雨

    (哈爾濱工程大學(xué) 船舶工程學(xué)院, 哈爾濱 150001)

    隨著計(jì)算機(jī)的發(fā)展,計(jì)算流體力學(xué)(Computational Fluid Dynamics, CFD)被越來(lái)越多地應(yīng)用到海洋結(jié)構(gòu)物的性能研究中.潛艇是重要的海洋結(jié)構(gòu)物,準(zhǔn)確了解其不同航速下的水動(dòng)力性能、精確掌握其流場(chǎng)特征,對(duì)提高潛艇自身性能具有重要意義[1].

    國(guó)內(nèi)外學(xué)者針對(duì)潛艇性能開(kāi)展過(guò)大量的數(shù)值研究.Pan等[2]采用RANSk-ω(k為湍動(dòng)能,ω為單位湍動(dòng)能耗散率)和雷諾應(yīng)力湍流模型(Reynolds Stress Model, RSM)對(duì)潛艇的阻力進(jìn)行預(yù)報(bào),取得了較好的結(jié)果.Posa等[3-5]基于大渦模擬湍流模型(Large Eddy Simulation, LES)對(duì)潛艇的尾流場(chǎng)進(jìn)行仿真計(jì)算,結(jié)果表明LES能很好地捕捉流動(dòng)特征,同時(shí)為流噪聲的預(yù)報(bào)提供精確輸入.Fu等[6-8]就分離渦湍流模型(Detach Eddy Simulation, DES)下的潛艇水動(dòng)力性能及流場(chǎng)特征進(jìn)行仿真,結(jié)果證實(shí)前者對(duì)垂直壁面方向的無(wú)量綱化距離y+的敏感性微弱,后者對(duì)空間離散精度的高需求性.此外,李士強(qiáng)等[8]還就k-ω、k-ε(ε為湍流耗散率)模型對(duì)潛艇附體變動(dòng)的影響進(jìn)行仿真,并證實(shí)關(guān)于舵失速的問(wèn)題,在失速前k-ε獲取的渦流脫落頻率較k-ω更高.

    綜上,國(guó)內(nèi)外學(xué)者就直航潛艇的性能進(jìn)行較多分析,但未闡明不同湍流模型在潛艇水動(dòng)力性能及流場(chǎng)特征計(jì)算中的優(yōu)勢(shì)或短板.本文針對(duì)這一問(wèn)題,首先展多套網(wǎng)格及時(shí)間步長(zhǎng)的收斂性分析.接著,針對(duì)不同的湍流模型開(kāi)展特定工況的數(shù)值仿真計(jì)算,比較預(yù)報(bào)結(jié)果的差異.最后,基于LES-Smagorinsky模型的預(yù)報(bào)結(jié)果對(duì)潛艇渦系特征進(jìn)行分析.

    1 理論基礎(chǔ)

    1.1 控制方程

    時(shí)均化后不可壓縮流體的連續(xù)性方程與動(dòng)量方程為

    (1)

    (2)

    1.2 湍流模型

    由于在時(shí)間和空間上對(duì)計(jì)算成本的需求較高,目前主要的CFD求解器均采用求解平均量或?yàn)V波量的方法來(lái)模擬小尺度流動(dòng)以節(jié)省計(jì)算成本.本文就多個(gè)湍流模型開(kāi)展?jié)撏У乃畡?dòng)力性能及流場(chǎng)特征計(jì)算,對(duì)比各種湍流模型的計(jì)算結(jié)果,進(jìn)而篩選適合潛艇水動(dòng)力性能及流場(chǎng)特征仿真計(jì)算的湍流模型.

    雷諾平均納維爾-斯托克斯(Reynolds-Averaged Navier-Stokes, RANS)是一類湍流模型集合的總稱,涉及多個(gè)子模型,這些模型的特征都是在N-S(Navier-Stokes)方程的基礎(chǔ)上添加了一個(gè)額外的黏度項(xiàng)且都包含一個(gè)通用項(xiàng)湍動(dòng)能k,即單位質(zhì)量湍流脈動(dòng)的動(dòng)能,而這些模型的區(qū)別在于對(duì)額外的黏度項(xiàng)處理方式不同:SA(Spalart-Allmaras)模型是一方程湍流模型,無(wú)壁面函數(shù)要求但增加一個(gè)新的變量SA黏度,其優(yōu)點(diǎn)是對(duì)計(jì)算內(nèi)存要求較低,有良好的收斂性,缺點(diǎn)是不適合求解剪切流與分離流;k-ε模型基于k和ε建立,是具有良好收斂性的兩方程模型,自提出后逐漸演化出realizablek-ε,RNGk-ε等,局限性體現(xiàn)在對(duì)無(wú)滑移壁面、逆壓梯度等的流動(dòng)模擬不精確,耗散不易計(jì)算;k-ω模型計(jì)算結(jié)果對(duì)初始條件較為敏感且收斂相對(duì)困難,但相對(duì)于ε,ω更容易求解.LES是基于自相似理論,通過(guò)過(guò)濾函數(shù)將大渦和小渦分離開(kāi),大渦直接解析,小渦是各項(xiàng)同性,因此可對(duì)其進(jìn)行模化計(jì)算,在解析大渦和?;u的過(guò)程中,產(chǎn)生的亞格子尺度應(yīng)力項(xiàng)通過(guò)亞格子(Subgrid Scale, SGS) 模型計(jì)算.DES模型結(jié)合RANS/LES模型各自的優(yōu)點(diǎn),在近壁面采用RANS進(jìn)行計(jì)算而利用LES模型的算法對(duì)湍流核心區(qū)域計(jì)算,經(jīng)過(guò)多年發(fā)展逐漸演變出延遲分離渦湍流 (Delayed Detach Eddy Simulation, DDES)、改進(jìn)的延遲分離渦湍流 (Improved Delayed Detached Eddy Simulation, IDDES) 等模型變體,對(duì)非定常湍流流動(dòng)的模擬有較大的優(yōu)勢(shì).

    本文基于植入到STAR-CCM+中的湍流模型:RANSk-ω、k-ε、SA、LES-Samgorinsky(為行文簡(jiǎn)潔,后續(xù)以LES-S替代)、WALE(LES-W)、Dynamic Smagorinsky(LES-DS)以及DES-SA、DES-EBKE、DDES、IDDES模型開(kāi)展Sub-off水動(dòng)力性能及流場(chǎng)特征的計(jì)算,進(jìn)而甄選合適的湍流模型開(kāi)展?jié)撏畡?dòng)力性能及流場(chǎng)特征的計(jì)算.

    2 仿真基礎(chǔ)

    2.1 坐標(biāo)系和計(jì)算域

    本文采用的幾何模型是Sub-off模型,詳細(xì)的模型資料可參考文獻(xiàn)[9-11],后續(xù)的研究中,在原始幾何模型尺度的基礎(chǔ)上進(jìn)行縮放,縮放比λ=0.688 7,模型總長(zhǎng)L=3 m,最大直徑D=0.35 m,若無(wú)特殊說(shuō)明本文的計(jì)算結(jié)果均是基于這一參數(shù)下的幾何模型.

    本文設(shè)定笛卡爾坐標(biāo)系Oxyz,其原點(diǎn)位于艇首,如圖1(a)所示,圖中Uref為自由流場(chǎng)的自由來(lái)流速度.本文設(shè)定計(jì)算域?yàn)閳A柱,具體尺寸及邊界條件的設(shè)置見(jiàn)圖1(b),參考文獻(xiàn)[12]的研究,當(dāng)目標(biāo)幾何模型在計(jì)算域流向投影面積小于1%時(shí),數(shù)值計(jì)算結(jié)果可忽略阻塞效應(yīng)的影響,通過(guò)簡(jiǎn)單的數(shù)據(jù)計(jì)算可知,Sub-off模型在計(jì)算域上下底面的投影面積占比為0.7%(小于1%),故可忽略阻塞效應(yīng)的影響.

    圖1 幾何模型及計(jì)算域

    2.2 網(wǎng)格劃分和工況設(shè)定

    圖2 網(wǎng)格方案3簡(jiǎn)圖

    為統(tǒng)一和方便作者進(jìn)行數(shù)值和試驗(yàn)研究,結(jié)合計(jì)算和實(shí)際條件需求,對(duì)航速為2.53 m/s的潛艇模型的水動(dòng)力性能及流場(chǎng)特征進(jìn)行收斂性計(jì)算,以對(duì)比不同湍流模型獲得的仿真結(jié)果.

    3 收斂性分析

    表2 4套網(wǎng)格方案的計(jì)算結(jié)果

    網(wǎng)格收斂率K=εi, i+1/εi, i-1,εi, i+1是較密兩方案結(jié)果差,εi, i-1是較稀兩方案結(jié)果差.K的計(jì)算結(jié)果可能是K>1,0

    表3 網(wǎng)格收斂性計(jì)算結(jié)果

    P=ln (εi, i+1/εi, i-1)lnδ

    評(píng)估3套網(wǎng)格的收斂指數(shù):

    Ii, i+1=DS|ei, i+1|/(δP-1)

    式中:DS=1.25為網(wǎng)格收斂安全系數(shù)[18];ei, i+1=εi, i+1.理論上,當(dāng)P=2時(shí),網(wǎng)格正交性或湍流模型非線性對(duì)計(jì)算結(jié)果準(zhǔn)確性的影響程度最小[14,18],選用中等密度網(wǎng)格進(jìn)行數(shù)值計(jì)算具有較高可靠性.可知方案2~4具有更好的收斂率且離散精度更接近理論值,故后續(xù)的計(jì)算中采用網(wǎng)格方案3進(jìn)行計(jì)算.

    表4是不同時(shí)間步長(zhǎng)計(jì)算得到的判定系數(shù)結(jié)果,以0.001 s的計(jì)算結(jié)果為基準(zhǔn)進(jìn)行對(duì)比.可以看出, 不同時(shí)間步長(zhǎng)對(duì)計(jì)算結(jié)果的影響較小,為提高計(jì)算效率,后續(xù)的計(jì)算中均采用0.001 s為時(shí)間步長(zhǎng).

    表4 不同時(shí)間步長(zhǎng)計(jì)算結(jié)果

    4 結(jié)果與討論

    4.1 水動(dòng)力性能

    不同湍流模型對(duì)近壁面流動(dòng)、黏性應(yīng)力等的解析不同,造成仿真得到的阻力有所差異,表5為不同湍流模型計(jì)算得到的潛艇阻力.

    表5 不同湍流模型的艇體阻力仿真結(jié)果

    表6 不同湍流模型計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比

    4.2 流場(chǎng)特征

    采用Q標(biāo)準(zhǔn)來(lái)判定流體為渦旋或應(yīng)變主導(dǎo)[19],圖3所示為不同湍流模型仿真計(jì)算得到的Sub-off周圍的渦場(chǎng).圖3清晰地表明LES-W、LES-DS對(duì)潛艇重要渦特征——馬蹄渦的捕捉不足,不同的是除此外的8個(gè)湍流模型均很好地捕捉到了圍殼附近馬蹄渦的流動(dòng)特征,另外這兩種湍流模型預(yù)報(bào)的潛艇尾肩部的渦結(jié)構(gòu)亦沒(méi)有其他模型強(qiáng)烈.選用的10種湍流模型中,RANS類湍流模型對(duì)潛艇尾部渦強(qiáng)度的預(yù)報(bào)結(jié)果顯然低于其他模型,這可能是RANS速度平均后的結(jié)果,是RANS模型的算法特點(diǎn)所決定.值得注意的是,10個(gè)模型均沒(méi)能很好地捕捉或預(yù)報(bào)圍殼稍渦的演化,但對(duì)于尾舵稍渦的演化均能較好地捕捉.推測(cè)原因:一是圍殼稍渦的強(qiáng)度沒(méi)有尾舵稍渦強(qiáng)度大,二是空間分辨率或時(shí)間分辨率不足以捕捉其演化過(guò)程.限于計(jì)算資源以及本文研究重點(diǎn),在此不做深入探討.

    圖3 不同湍流模型計(jì)算艇體周圍Q采樣體積(Q=5 s-2)

    伴流是考核潛艇綜合性能的重要指標(biāo),對(duì)伴流的準(zhǔn)確預(yù)報(bào)十分重要.圖4所示為不同湍流模型計(jì)算得到的Sub-off尾流伴流分?jǐn)?shù)δ隨距離變化的曲線,圖中,0為潛艇中軸線的最尾端,D、2D、3D、4D為沿著x正方向與潛艇中軸線最尾端的距離.可見(jiàn),LES-DS、LES-W的預(yù)報(bào)結(jié)果與其他模型有明顯偏差,但同屬大渦模擬類中LES-S的預(yù)報(bào)結(jié)果有良好的“隨群”性,結(jié)合對(duì)于宏觀力的預(yù)報(bào)結(jié)果可知,LES-S模型在預(yù)報(bào)潛艇宏觀變量方面有著較好的適用性.

    圖4 不同湍流模型預(yù)報(bào)尾流伴流(0位于艇體尾端)

    為了深入探究不同湍流模型對(duì)Sub-off流場(chǎng)的預(yù)報(bào)結(jié)果, 本文提取了如圖5虛線所示位置的流場(chǎng)時(shí)均變量并加以分析,時(shí)均流場(chǎng)的采樣總長(zhǎng)是流經(jīng)一倍艇長(zhǎng)時(shí)間. 圖6所示為時(shí)均軸向速度的周向變化曲線,圖中θ為圖5中提取截面沿著白色虛線環(huán)向位置的角度,u為流場(chǎng)局部軸向速度,通過(guò)u/Uref對(duì)軸向速度進(jìn)行無(wú)因次化對(duì)比, 分析其變化規(guī)律,為了有效說(shuō)明流場(chǎng),將Sub-off周圍流場(chǎng)分為內(nèi)流場(chǎng)(r/D<0.501)和外流場(chǎng)(r/D>0.501),其中r為對(duì)應(yīng)曲線的半徑值.可知,10種湍流模型對(duì)Sub-off的外流場(chǎng)及其平行中體內(nèi)流場(chǎng)的預(yù)報(bào)結(jié)果吻合程度非常高,推測(cè)由于遠(yuǎn)離壁面或壁面流動(dòng)簡(jiǎn)單,各湍流模型對(duì)其處理無(wú)明顯差異,但即使對(duì)平行中體內(nèi)流場(chǎng)的流動(dòng)預(yù)報(bào)差異微弱,Sub-off尾內(nèi)流場(chǎng)的預(yù)報(bào)差異卻十分明顯.這表明,湍流模型對(duì)因?yàn)閹缀翁卣髯兓l(fā)流動(dòng)變化更為敏感,對(duì)流場(chǎng)“傳遞過(guò)程”引發(fā)的變化敏感性較弱.圖7所示為流場(chǎng)時(shí)均壓力的周向分布曲線,可以看出,Sub-off周圍流場(chǎng)的時(shí)均壓力系數(shù)Cp(Cp=p/(0.5ρUref))分布曲線對(duì)湍流模型地變化不敏感,各個(gè)模型均有良好的“隨群”性,由于尾部附體以及尾部幾何收縮造成的負(fù)壓力梯度,在x/L=0.9處,流場(chǎng)壓力波動(dòng)仍較為劇烈.圖6和7的結(jié)果表明:時(shí)均速度較時(shí)均壓力對(duì)湍流模型的改變更敏感,尤其是近壁面流動(dòng),如圖6(a)、7(a)所示,這是不同湍流模型調(diào)用不同壁面解析函數(shù)造成的.圖6(b)、7(b)中不同湍流模型得到的流場(chǎng)流速和壓力的周向分布基本相同,說(shuō)明在遠(yuǎn)流場(chǎng)各個(gè)湍流模型預(yù)報(bào)的數(shù)值仿真結(jié)果基本相同,但圖7(b)中x/L=0.90處的流速和壓力分布仍有相對(duì)明顯的波動(dòng),這是由于艇體尾部逆壓力梯度造成尾部流速降低,邊界層增厚,從而改變流場(chǎng)流速和壓力基本對(duì)應(yīng)的變化趨勢(shì).

    圖5 流場(chǎng)提取截面

    圖6 流場(chǎng)時(shí)均軸向速度的周向分布

    圖7 Cp的周向分布

    據(jù)上推斷,10種湍流模型在Sub-off流場(chǎng)時(shí)均變量的捕捉上存在差異,且主要應(yīng)該體現(xiàn)在流場(chǎng)黏性上,進(jìn)而導(dǎo)致其他流場(chǎng)變量的變化,在本文采用的時(shí)間及空間分辨率的情況下,LES-DE、LES-W同比效果更差,原因可能是空間或時(shí)間分辨率不足導(dǎo)致的湍流模型對(duì)流場(chǎng)黏性特征的預(yù)報(bào)存在差異.

    不同湍流模型對(duì)潛艇周圍流場(chǎng)參數(shù)的平均量預(yù)測(cè)結(jié)果存在差異,同時(shí)對(duì)流場(chǎng)的二次量如雷諾應(yīng)力的預(yù)測(cè)結(jié)果亦有不同.圖8所示為不同湍流模型仿真得到的潛艇周圍流場(chǎng)的雷諾正應(yīng)力RNS,其余湍流模型的結(jié)果由于數(shù)值過(guò)小,出于文章簡(jiǎn)潔的考慮未在文中列出.如同文獻(xiàn)[17],LES模型對(duì)流場(chǎng)二次特征的捕捉較好,DES模型次之,RANS模型最差,故研究流場(chǎng)二次量時(shí)模型的選擇會(huì)引起數(shù)值仿真結(jié)果的巨大差異.

    綜上可知,就模型尺度下的潛艇而言,RANS類湍流模型不能有效地捕捉潛艇流場(chǎng)的流動(dòng)細(xì)節(jié),對(duì)湍流應(yīng)力有明顯的抑制;DES類湍流模型普遍會(huì)抑制潛艇尾流的流動(dòng)渦結(jié)構(gòu),同時(shí)預(yù)報(bào)的宏觀力亦存在明顯的偏差;LES類湍流模型在捕捉潛艇流動(dòng)渦結(jié)構(gòu)、湍流應(yīng)力方面表現(xiàn)優(yōu)異,但是除LES-S模型外,針對(duì)宏觀力的預(yù)測(cè)均出現(xiàn)明顯的誤差.因此在開(kāi)展?jié)撏畡?dòng)力及流場(chǎng)特征分析的綜合研究中,LES-S模型是一個(gè)較為優(yōu)異的選擇.

    4.3 潛艇周圍渦場(chǎng)特征

    基于上述分析得到的數(shù)值仿真策略對(duì)Sub-off的流場(chǎng)特征進(jìn)行計(jì)算.由于艇體表面附體的存在,使得其周圍流場(chǎng)較為復(fù)雜.圖9所示為潛艇周圍的渦系簡(jiǎn)圖,渦系通過(guò)Q標(biāo)準(zhǔn)的重新采樣體積進(jìn)行表達(dá),Q標(biāo)準(zhǔn)是Hunt等[19]提出的一項(xiàng)渦的判斷方法.其中附著在圍殼或尾舵附近的稍渦、馬蹄渦(項(xiàng)鏈渦)會(huì)改變螺旋槳的來(lái)流均勻性而增加螺旋槳的載荷波動(dòng)、提高空化風(fēng)險(xiǎn)及增加水動(dòng)力噪聲,湍流度巨大的渦結(jié)構(gòu)會(huì)提高艇體表面激勵(lì)力,引發(fā)艇身結(jié)構(gòu)振動(dòng),削弱潛艇的隱身性能.艇體附著渦是由邊界層流動(dòng)引起的速度梯度變化形成,在邊界層較為穩(wěn)定的平行中體段強(qiáng)度微弱,在幾何特征變化明顯的首尾段有較強(qiáng)的艇體附著渦.圖10所示為圍殼附近多截面流場(chǎng)變量,圖中:wi,wj,wk為渦量在3個(gè)方向的分量,圖11對(duì)比圍殼和尾舵表面的壓力變化.從圖10可以看出,圍殼在軸向首先出現(xiàn)正壓力梯度,在尾段出現(xiàn)負(fù)壓力梯度,而其附近的渦量表明稍渦的脫落和壓力梯度的變化存在密切關(guān)系,同時(shí)值得注意的是wk可能是稍渦變化的重要變量,存在明顯脫落和分離的wi、wj是圍殼隨邊渦變化產(chǎn)生的重要變量.由于尾舵和圍殼相似的幾何特征,本文認(rèn)為尾舵稍渦和圍殼稍渦的形成及影響因素相似,但尾舵稍渦的強(qiáng)度明顯大于圍殼稍渦.從圖11可以看出,圍殼和尾舵頂面和側(cè)面的壓力梯度差異可能是導(dǎo)致強(qiáng)度不同的關(guān)鍵因素, 需要注意的是圖10和11中均對(duì)變量進(jìn)行無(wú)因次化處理以方便對(duì)比分析.文獻(xiàn)[17]的研究表明馬蹄渦的存在與圍殼(尾舵)、艇體邊界層的三維流動(dòng)有關(guān),本文不對(duì)其進(jìn)行深入研究.

    圖9 Sub-off周圍渦系

    圖10 圍殼附近多截面流場(chǎng)變量

    圖11 圍殼和尾舵表面Cp分布

    5 結(jié)語(yǔ)

    本文以Sub-off為幾何模型,針對(duì)特定工況,就10種湍流模型計(jì)算獲得的結(jié)果進(jìn)行分析,得出如下結(jié)論:

    通過(guò)空間及時(shí)間收斂性分析,基于理查森外推法判定本文選用的網(wǎng)格方案具有良好的網(wǎng)格收斂性,最終采用了合適的網(wǎng)格開(kāi)展后續(xù)的數(shù)值計(jì)算.本文選用的湍流模型中,LES-S模型在潛艇的宏觀量如阻力、時(shí)均流場(chǎng)、伴流方面具有良好的性能或“隨群性”;同時(shí),在流場(chǎng)的二次量預(yù)報(bào)方面,LES類模型獲得的結(jié)果較其他模型更加豐富.

    潛艇渦系因其表面附體的存在較為復(fù)雜,但由于圍殼和尾舵的幾何相似性,其渦系具有內(nèi)部形成機(jī)理的相似性,但因?yàn)楦襟w裝配位置及參數(shù)的差異,尾舵的渦系耦合和演化有較大差異.

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    好男人视频免费观看在线| 国产一区二区亚洲精品在线观看| 午夜福利视频精品| 男人狂女人下面高潮的视频| 80岁老熟妇乱子伦牲交| 日韩欧美 国产精品| 亚洲国产av新网站| 亚洲国产精品成人综合色| a级毛色黄片| 亚洲丝袜综合中文字幕| 亚洲成色77777| 亚洲精品成人久久久久久| 嫩草影院精品99| 美女内射精品一级片tv| 最近2019中文字幕mv第一页| 欧美丝袜亚洲另类| a级毛片免费高清观看在线播放| 蜜臀久久99精品久久宅男| 黄片wwwwww| 亚洲欧美日韩无卡精品| 你懂的网址亚洲精品在线观看| 三级国产精品片| 欧美丝袜亚洲另类| 香蕉精品网在线| 国产成人精品一,二区| 日韩成人伦理影院| 你懂的网址亚洲精品在线观看| 欧美+日韩+精品| 日日摸夜夜添夜夜添av毛片| 日韩电影二区| 在线亚洲精品国产二区图片欧美 | 久久韩国三级中文字幕| 亚洲va在线va天堂va国产| 卡戴珊不雅视频在线播放| 99久国产av精品国产电影| 欧美丝袜亚洲另类| av黄色大香蕉| 一级爰片在线观看| 亚洲av中文字字幕乱码综合| 熟女av电影| 日韩在线高清观看一区二区三区| 成年版毛片免费区| 又大又黄又爽视频免费| 国产高清有码在线观看视频| 欧美成人午夜免费资源| 精品国产一区二区三区久久久樱花 | 亚洲精品自拍成人| 爱豆传媒免费全集在线观看| 国内精品宾馆在线| a级一级毛片免费在线观看| 国产成人精品久久久久久| 亚洲一级一片aⅴ在线观看| 亚洲婷婷狠狠爱综合网| 天堂网av新在线| 亚洲成色77777| 99久久精品热视频| 一本一本综合久久| 亚洲丝袜综合中文字幕| av播播在线观看一区| 中国国产av一级| 韩国av在线不卡| 亚洲精华国产精华液的使用体验| 精品人妻视频免费看| 青青草视频在线视频观看| 欧美3d第一页| 成人亚洲欧美一区二区av| 亚洲国产av新网站| 国产老妇伦熟女老妇高清| 在线观看一区二区三区激情| 身体一侧抽搐| 91久久精品国产一区二区三区| 久久久久国产精品人妻一区二区| 尾随美女入室| 精品国产三级普通话版| 欧美日韩国产mv在线观看视频 | 99精国产麻豆久久婷婷| 高清av免费在线| av女优亚洲男人天堂| 亚洲自拍偷在线| 国产精品av视频在线免费观看| 在线观看免费高清a一片| 亚洲美女搞黄在线观看| 亚洲第一区二区三区不卡| 各种免费的搞黄视频| 最近手机中文字幕大全| 午夜视频国产福利| 免费观看a级毛片全部| 最近中文字幕高清免费大全6| 如何舔出高潮| 日本免费在线观看一区| 各种免费的搞黄视频| 国产淫片久久久久久久久| 亚洲丝袜综合中文字幕| 亚洲天堂av无毛| 亚洲欧美日韩卡通动漫| 中文字幕人妻熟人妻熟丝袜美| 人妻制服诱惑在线中文字幕| 22中文网久久字幕| 51国产日韩欧美| 久久久久精品性色| 最后的刺客免费高清国语| 欧美zozozo另类| 王馨瑶露胸无遮挡在线观看| 禁无遮挡网站| 亚洲成人一二三区av| 国产淫语在线视频| 久久久久久久大尺度免费视频| 人人妻人人看人人澡| 国产永久视频网站| 男女边吃奶边做爰视频| 免费观看的影片在线观看| 国产成人aa在线观看| 国产欧美日韩精品一区二区| 人妻夜夜爽99麻豆av| 国产一区亚洲一区在线观看| 性色avwww在线观看| 狠狠精品人妻久久久久久综合| 人妻 亚洲 视频| 亚洲欧美一区二区三区黑人 | 精品午夜福利在线看| 日本猛色少妇xxxxx猛交久久| 一二三四中文在线观看免费高清| 乱系列少妇在线播放| 少妇熟女欧美另类| 性色av一级| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人精品一区久久| 麻豆成人av视频| 久久久久国产网址| 国产精品一区二区在线观看99| 亚洲av国产av综合av卡| 高清日韩中文字幕在线| 亚洲,一卡二卡三卡| 性色avwww在线观看| 麻豆乱淫一区二区| 久久精品国产a三级三级三级| 国产探花在线观看一区二区| 一级毛片黄色毛片免费观看视频| 国产午夜精品久久久久久一区二区三区| 久久精品综合一区二区三区| 热99国产精品久久久久久7| 欧美 日韩 精品 国产| 在线观看一区二区三区激情| 久久久欧美国产精品| 亚洲性久久影院| 亚洲av欧美aⅴ国产| 日韩三级伦理在线观看| 在线观看免费高清a一片| 日韩一区二区三区影片| 涩涩av久久男人的天堂| 久久久久久久久久成人| 成人国产av品久久久| 女人久久www免费人成看片| 国产精品久久久久久精品电影| 一本一本综合久久| 国产毛片a区久久久久| 国产精品一及| 成年免费大片在线观看| 国产有黄有色有爽视频| 免费大片18禁| 日韩视频在线欧美| 午夜福利在线在线| 国产亚洲精品久久久com| 99久久精品一区二区三区| 亚洲成人精品中文字幕电影| 国产精品熟女久久久久浪| kizo精华| 亚洲精品色激情综合| 在线免费观看不下载黄p国产| 日韩国内少妇激情av| 精品久久久久久电影网| 青春草国产在线视频| 欧美最新免费一区二区三区| 亚洲精品国产成人久久av| 狂野欧美激情性bbbbbb| 国产乱来视频区| 久久久久久久久大av| 国产成年人精品一区二区| 国产黄片视频在线免费观看| 最近的中文字幕免费完整| 一级黄片播放器| 精品99又大又爽又粗少妇毛片| 搞女人的毛片| 亚洲成人中文字幕在线播放| 激情 狠狠 欧美| 国产亚洲精品久久久com| 国产一区有黄有色的免费视频| 伦理电影大哥的女人| 爱豆传媒免费全集在线观看| 女的被弄到高潮叫床怎么办| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 国产在线男女| 少妇人妻久久综合中文| 欧美xxxx性猛交bbbb| 精品酒店卫生间| 九色成人免费人妻av| 极品教师在线视频| 麻豆成人午夜福利视频| 一个人看的www免费观看视频| 少妇的逼好多水| 亚洲精品中文字幕在线视频 | 能在线免费看毛片的网站| 免费人成在线观看视频色| 精品人妻视频免费看| 日韩免费高清中文字幕av| 免费av观看视频| 亚洲美女视频黄频| 热re99久久精品国产66热6| 欧美最新免费一区二区三区| 午夜亚洲福利在线播放| 99久久精品一区二区三区| 国内揄拍国产精品人妻在线| 建设人人有责人人尽责人人享有的 | 校园人妻丝袜中文字幕| 高清在线视频一区二区三区| 日韩av不卡免费在线播放| 精品少妇久久久久久888优播| 国产精品国产三级国产专区5o| www.av在线官网国产| 亚洲怡红院男人天堂| 大香蕉97超碰在线| 亚洲精品456在线播放app| 国产探花极品一区二区| 欧美日韩视频精品一区| 亚洲精品亚洲一区二区| 精品一区二区免费观看| 国国产精品蜜臀av免费| av在线播放精品| 国产 精品1| 国产淫语在线视频| 国产精品99久久久久久久久| 日韩,欧美,国产一区二区三区| 亚洲自拍偷在线| 久久久精品欧美日韩精品| 毛片一级片免费看久久久久| 免费观看的影片在线观看| 建设人人有责人人尽责人人享有的 | 日本三级黄在线观看| 日本熟妇午夜| 国产伦精品一区二区三区视频9| 禁无遮挡网站| 亚洲丝袜综合中文字幕| 成人一区二区视频在线观看| 亚洲精品国产成人久久av| 亚洲欧洲日产国产| 成人国产麻豆网| 亚洲av电影在线观看一区二区三区 | 国产精品国产三级专区第一集| 97超视频在线观看视频| 欧美xxxx黑人xx丫x性爽| 国产v大片淫在线免费观看| 中国美白少妇内射xxxbb| 亚洲av中文av极速乱| 91久久精品国产一区二区成人| 午夜免费观看性视频| 日韩制服骚丝袜av| 国产爽快片一区二区三区| 午夜日本视频在线| 国产精品久久久久久精品电影小说 | 国产伦在线观看视频一区| 一级毛片 在线播放| 亚洲精品中文字幕在线视频 | 成人午夜精彩视频在线观看| 欧美精品国产亚洲| 国产精品麻豆人妻色哟哟久久| 日本色播在线视频| 干丝袜人妻中文字幕| 制服丝袜香蕉在线| 美女cb高潮喷水在线观看| 亚洲国产日韩一区二区| 99九九线精品视频在线观看视频| 小蜜桃在线观看免费完整版高清| 国产成人精品婷婷| 亚洲在线观看片| 亚洲天堂av无毛| 在线 av 中文字幕| 中国美白少妇内射xxxbb| 亚洲av成人精品一二三区| 中文资源天堂在线| av天堂中文字幕网| 哪个播放器可以免费观看大片| 激情 狠狠 欧美| 极品教师在线视频| 三级经典国产精品| 国产成人免费观看mmmm| 日本爱情动作片www.在线观看| 超碰97精品在线观看| 亚洲国产日韩一区二区| 精品一区二区免费观看| 国产精品精品国产色婷婷| 性色avwww在线观看| 狠狠精品人妻久久久久久综合| 久久精品国产a三级三级三级| 亚洲最大成人中文| 久久久久久九九精品二区国产| 最新中文字幕久久久久| 卡戴珊不雅视频在线播放| 亚洲美女搞黄在线观看| 肉色欧美久久久久久久蜜桃 | 免费大片18禁| 精品99又大又爽又粗少妇毛片| 人人妻人人澡人人爽人人夜夜| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲一区二区精品| 亚洲精品自拍成人| 久久精品夜色国产| 男女啪啪激烈高潮av片| 亚洲最大成人手机在线| 成年女人看的毛片在线观看| 亚洲综合精品二区| 老司机影院毛片| 天美传媒精品一区二区| 大又大粗又爽又黄少妇毛片口| 国产黄色免费在线视频| 夜夜爽夜夜爽视频| 免费观看av网站的网址| 一边亲一边摸免费视频| 如何舔出高潮| 国产高清不卡午夜福利| 国产成年人精品一区二区| 午夜老司机福利剧场| 国产爽快片一区二区三区| 亚洲成人av在线免费| 亚洲综合色惰| 亚洲高清免费不卡视频| 天堂俺去俺来也www色官网| www.av在线官网国产| av免费在线看不卡| 在线 av 中文字幕| 熟女av电影| www.色视频.com| 亚洲欧美中文字幕日韩二区| 99久久九九国产精品国产免费| 一级二级三级毛片免费看| 亚洲av中文字字幕乱码综合| 有码 亚洲区| 精品久久久噜噜| 国产精品久久久久久精品电影| 中文字幕免费在线视频6| 国产精品国产三级专区第一集| av女优亚洲男人天堂| 欧美xxⅹ黑人| 亚洲欧美清纯卡通| 成年女人看的毛片在线观看| 国产亚洲av嫩草精品影院| av天堂中文字幕网| 欧美成人一区二区免费高清观看| 大香蕉久久网| 69av精品久久久久久| 日韩国内少妇激情av| 色哟哟·www| 丝袜喷水一区| 精品午夜福利在线看| av在线app专区| 亚洲国产av新网站| 亚洲精品一区蜜桃| 最近最新中文字幕免费大全7| 丝瓜视频免费看黄片| 少妇被粗大猛烈的视频| 免费看av在线观看网站| 免费观看a级毛片全部| 久久精品国产亚洲网站| 亚洲自拍偷在线| 麻豆乱淫一区二区| 国产精品偷伦视频观看了| 日韩一区二区三区影片| 自拍偷自拍亚洲精品老妇| 日韩一区二区三区影片| freevideosex欧美| 免费看av在线观看网站| 成人黄色视频免费在线看| 性色av一级| 国产亚洲5aaaaa淫片| 一级黄片播放器| 精品人妻偷拍中文字幕| 少妇 在线观看| av网站免费在线观看视频| av女优亚洲男人天堂| 身体一侧抽搐| 日韩 亚洲 欧美在线| 成人毛片a级毛片在线播放| 日韩亚洲欧美综合| 日韩一本色道免费dvd| 97热精品久久久久久| av天堂中文字幕网| 久久久久九九精品影院| 国产伦在线观看视频一区| 日本猛色少妇xxxxx猛交久久| 欧美成人精品欧美一级黄| 一个人看的www免费观看视频| 亚洲内射少妇av| 国产精品久久久久久久久免| 肉色欧美久久久久久久蜜桃 | 在现免费观看毛片| 午夜视频国产福利| eeuss影院久久| 亚洲四区av| 五月开心婷婷网| 国产成人福利小说| 欧美成人午夜免费资源| 国产精品国产三级专区第一集| 婷婷色综合www| 高清午夜精品一区二区三区| 99久久精品国产国产毛片| 日本-黄色视频高清免费观看| 欧美日韩国产mv在线观看视频 | 99热网站在线观看| 国产高清三级在线| 精品久久久久久久久av| 少妇的逼水好多| av在线天堂中文字幕| 91狼人影院| 免费观看性生交大片5| 伊人久久精品亚洲午夜| 欧美成人a在线观看| 伦理电影大哥的女人| 亚洲高清免费不卡视频| 肉色欧美久久久久久久蜜桃 | 国产真实伦视频高清在线观看| 美女高潮的动态| 亚洲欧美日韩无卡精品| 午夜精品国产一区二区电影 | a级毛色黄片| 街头女战士在线观看网站| 如何舔出高潮| 精品久久久久久电影网| 中文字幕久久专区| 国产亚洲午夜精品一区二区久久 | 丝袜美腿在线中文| 看十八女毛片水多多多| 亚洲va在线va天堂va国产| 在线看a的网站| 国产精品久久久久久精品电影小说 | 国产精品.久久久| 一个人观看的视频www高清免费观看| 亚洲性久久影院| 综合色av麻豆| 午夜福利在线观看免费完整高清在| 久久精品国产亚洲av涩爱| 久久久国产一区二区| 观看免费一级毛片| 国产精品福利在线免费观看| 最后的刺客免费高清国语| 黄色配什么色好看| 亚洲精品乱久久久久久| 美女高潮的动态| 三级国产精品欧美在线观看| 69av精品久久久久久| 国产久久久一区二区三区| av在线蜜桃| 亚洲精品国产av成人精品| 亚洲经典国产精华液单| 久久人人爽av亚洲精品天堂 | av在线app专区| 精华霜和精华液先用哪个| 建设人人有责人人尽责人人享有的 | 一区二区三区四区激情视频| 久久久久国产网址| 一本久久精品| 男人舔奶头视频| 久久精品久久久久久久性| 日本色播在线视频| 超碰97精品在线观看| 亚洲欧美日韩无卡精品| av国产久精品久网站免费入址| 美女被艹到高潮喷水动态| 中国国产av一级| 亚洲国产日韩一区二区| 99久久精品国产国产毛片| 午夜福利高清视频| 婷婷色av中文字幕| 精品一区二区免费观看| 亚洲精品日本国产第一区| 寂寞人妻少妇视频99o| 精品人妻熟女av久视频| av国产免费在线观看| 熟妇人妻不卡中文字幕| 国产av不卡久久| 免费看不卡的av| 一二三四中文在线观看免费高清| 免费观看av网站的网址| 国产精品一二三区在线看| 黄片wwwwww| 国产一区有黄有色的免费视频| 欧美国产精品一级二级三级 | av专区在线播放| 成人漫画全彩无遮挡| 看免费成人av毛片| 国产亚洲精品久久久com| 精品人妻偷拍中文字幕| 久久久久久久国产电影| 大片电影免费在线观看免费| 在线观看一区二区三区激情| 欧美国产精品一级二级三级 | 国产精品不卡视频一区二区| 欧美zozozo另类| 国产免费福利视频在线观看| 99热6这里只有精品| 99热国产这里只有精品6| 精品亚洲乱码少妇综合久久| 久久精品综合一区二区三区| 美女被艹到高潮喷水动态| 久久99热这里只频精品6学生| 九九在线视频观看精品| 精品人妻熟女av久视频| 久久久久精品久久久久真实原创| 欧美一级a爱片免费观看看| 国产精品精品国产色婷婷| 黄色配什么色好看| 成年人午夜在线观看视频| 国产极品天堂在线| 又爽又黄a免费视频| 一级片'在线观看视频| 欧美激情国产日韩精品一区| 日本一二三区视频观看| 蜜臀久久99精品久久宅男| 在线亚洲精品国产二区图片欧美 | 国产精品国产av在线观看| 日本爱情动作片www.在线观看| 欧美日韩一区二区视频在线观看视频在线 | 如何舔出高潮| 一本色道久久久久久精品综合| 精品久久久久久久人妻蜜臀av| 一二三四中文在线观看免费高清| 三级国产精品欧美在线观看| 91精品国产九色| 亚洲自偷自拍三级| 人妻一区二区av| 男插女下体视频免费在线播放| 久久精品人妻少妇| 亚洲美女搞黄在线观看| 麻豆国产97在线/欧美| 伊人久久精品亚洲午夜| 在线 av 中文字幕| 免费观看a级毛片全部| 夜夜爽夜夜爽视频| 久久6这里有精品| 国产免费一级a男人的天堂| 亚洲成人一二三区av| 国产欧美另类精品又又久久亚洲欧美| 人妻一区二区av| 高清欧美精品videossex| 黄片无遮挡物在线观看| 国产日韩欧美在线精品| 国产高清有码在线观看视频| 国产成人精品久久久久久| 免费看日本二区| 欧美精品国产亚洲| 国产精品不卡视频一区二区| 日韩免费高清中文字幕av| 街头女战士在线观看网站| 国产大屁股一区二区在线视频| 国产高清国产精品国产三级 | 一级毛片 在线播放| 欧美高清成人免费视频www| 日韩欧美精品免费久久| 精品人妻偷拍中文字幕| 在线亚洲精品国产二区图片欧美 | 夜夜爽夜夜爽视频| 日韩制服骚丝袜av| 直男gayav资源| 国产精品av视频在线免费观看| 国产在线一区二区三区精| 久久久久久伊人网av| 亚洲在线观看片| 九色成人免费人妻av| 亚洲精品一二三| 一级毛片久久久久久久久女| 久久久久国产网址| 亚洲高清免费不卡视频| 欧美三级亚洲精品| av女优亚洲男人天堂| 久久久午夜欧美精品| 韩国高清视频一区二区三区| 日本午夜av视频| 狠狠精品人妻久久久久久综合| 成人国产av品久久久| 乱系列少妇在线播放| 国产高潮美女av| 免费观看的影片在线观看| 亚洲精品日韩在线中文字幕| 精品人妻熟女av久视频| 精品人妻视频免费看| 国产爽快片一区二区三区| 中文精品一卡2卡3卡4更新| 男女无遮挡免费网站观看| 综合色av麻豆| 色综合色国产| 老司机影院成人| 99精国产麻豆久久婷婷| 五月天丁香电影| 免费高清在线观看视频在线观看| 亚洲精品国产色婷婷电影| 久久久久久久亚洲中文字幕| 亚洲真实伦在线观看| 免费看av在线观看网站| 五月伊人婷婷丁香| tube8黄色片| 欧美少妇被猛烈插入视频| 欧美激情久久久久久爽电影| 国产精品熟女久久久久浪| 亚洲av国产av综合av卡| 亚洲人成网站在线播| 99久久中文字幕三级久久日本| 国产精品人妻久久久久久| 91久久精品国产一区二区成人| 亚洲av日韩在线播放| 99热6这里只有精品| 少妇被粗大猛烈的视频| 日日啪夜夜爽| 欧美国产精品一级二级三级 | 舔av片在线| 亚洲最大成人中文| 午夜福利视频1000在线观看|