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

    基于RANS的suboff尾跡特征數(shù)值預(yù)報(bào)

    2018-11-05 01:30:52孟慶杰鄧海華
    艦船科學(xué)技術(shù) 2018年10期
    關(guān)鍵詞:內(nèi)波尾跡計(jì)算結(jié)果

    孟慶杰,彭 亮,鄧海華

    (武漢第二船舶設(shè)計(jì)研究所,湖北 武漢 430064)

    0 引 言

    水下航行器航行時(shí)引起的水動(dòng)力尾跡具有傳播范圍廣、幅度大、持續(xù)時(shí)間長、傳播距離在幾十千米外仍會(huì)引起可識別的特征等特點(diǎn),是造成其暴露的重要物理特征。利用合成孔徑雷達(dá)(SAR)對水下航行器水動(dòng)力尾跡進(jìn)行探測,并通過一定的方法可反演出水下航行器的位置、速度及潛深等重要信息,進(jìn)而可實(shí)現(xiàn)對水下航行器的準(zhǔn)確定位與監(jiān)視。因此,水下航行器水動(dòng)力尾跡探測技術(shù)是近年來各軍事強(qiáng)國優(yōu)先、重點(diǎn)發(fā)展的非聲探潛手段之一。

    水下航行器水動(dòng)力尾跡的準(zhǔn)確預(yù)報(bào)是開展SAR探潛技術(shù)研究的前提。國外學(xué)者圍繞水下航行器水動(dòng)力學(xué)尾跡產(chǎn)生機(jī)理、傳播特性及其水面特征等方面進(jìn)行了大量的研究[1–10]。Schooley[11]基于試驗(yàn)對均勻流體及存在垂向密度梯度的流體中,自驅(qū)動(dòng)物體的湍動(dòng)尾跡進(jìn)行研究,認(rèn)為粘性是引起內(nèi)波的主要原因。由于國內(nèi)對水下航行器尾跡特征研究滯后,目前主要采用細(xì)長體簡化模型水池試驗(yàn)進(jìn)行研究。魏崗等[12]對半球體模型運(yùn)動(dòng)生成的內(nèi)波問題進(jìn)行實(shí)驗(yàn)研究,試驗(yàn)結(jié)果表明傅汝德數(shù)對內(nèi)波形成影響顯著。姚志崇等[13]對密度分層流體中拖曳球體激發(fā)內(nèi)波特性進(jìn)行試驗(yàn),對體效應(yīng)穩(wěn)態(tài)內(nèi)波和尾流效應(yīng)非穩(wěn)態(tài)內(nèi)波的特征給出了直觀表達(dá),并獲得了內(nèi)波流場速度大小隨內(nèi)傅氏數(shù)的變化規(guī)律。

    盡管簡化模型試驗(yàn)[14–18]可針對研究問題給出直接的結(jié)果,但無法獲得流場所有信息,不利于深入研究水下航行器水動(dòng)力尾跡產(chǎn)生機(jī)理及其傳播特性。因此,近年來,CFD方法以其可與試驗(yàn)媲美的精度,并可獲得所有的流場細(xì)節(jié)等特點(diǎn),正逐漸成為研究水下航行器尾跡特征的主要方法之一[19–25]。Wessel[26],Young[27]及David[28]利用數(shù)值模擬方法,研究了均勻流體及分層流體中水下航行器尾流特征,認(rèn)為尾流結(jié)構(gòu)及其發(fā)展均會(huì)受到流動(dòng)方式及環(huán)境因素的影響。Sarakinos[29]利用代碼程序Galatea-I對不同來流方向情況下,DARPA suboff模型水動(dòng)力性能進(jìn)行數(shù)值仿真。Abdilghanie等[30]對拖曳物體尾流誘發(fā)的隨機(jī)內(nèi)波進(jìn)行了細(xì)致的數(shù)值模擬,基于計(jì)算結(jié)果,作者認(rèn)為尾流誘發(fā)的隨機(jī)內(nèi)波是一個(gè)持續(xù)且復(fù)雜的過程。

    本文將基于RANS方法,結(jié)合k-ε湍流模型及VOF方法,對suboff模型在均勻流體及密度強(qiáng)分層流體中運(yùn)動(dòng)時(shí)尾跡特征進(jìn)行預(yù)報(bào)。通過對比分析2種工況下suboff模型水動(dòng)力尾跡特征,以期實(shí)現(xiàn)suboff內(nèi)波尾跡與開爾文尾跡不同的直觀展現(xiàn),并實(shí)現(xiàn)流體分層對suboff尾跡影響方式的清晰展示,進(jìn)而為利用RANS方法預(yù)報(bào)suboff內(nèi)波尾跡特征提供可行性支持及參考數(shù)據(jù)。

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

    1.1 控制方程

    本文控制方程為三維不可壓RANS方程,慣性坐標(biāo)系下其方程表達(dá)式可表述為:

    1.2 湍流模型

    本文計(jì)算采用標(biāo)準(zhǔn)k-ε湍流模型實(shí)現(xiàn)控制方程的封閉。其湍動(dòng)能k方程為

    耗散率ε方程為

    1.3 VOF自由面捕捉

    VOF(Volume of Fluid)方法以流體占據(jù)網(wǎng)格單元的體積分?jǐn)?shù)來實(shí)現(xiàn)對自由面演化的跟蹤。在該方法中,所有流體滿足同一動(dòng)量方程,通過在整個(gè)計(jì)算區(qū)域上跟蹤每一種流體在計(jì)算單元中的體積分?jǐn)?shù)φq來確定交界面的位置。假設(shè)第q種流體在單元中的體積分?jǐn)?shù)為 φq,則

    φq=0,單元中不存在第q種流體;

    0<φq<1,單元中存在多種流體,即存在流體界面;

    φq=1,單元中全部為第q種流體。

    由于φq沿流體界面法向變化最快,因此確定界面法向及φq值后,就可確定一條分割線來表達(dá)流體界面位置。

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

    2.1 研究模型

    本文選用DARPA suboff標(biāo)模作為研究對象。該模型是美國國防先進(jìn)技術(shù)研究署(Defence Advanced Research Project Agency, DARPA)為建立suboff CFD分析軟件驗(yàn)證數(shù)據(jù)庫而專門設(shè)計(jì)的模型。如圖1所示,模型由水滴形主體部分、圍殼及4個(gè)對稱尾翼組成。表1給出了該模型的主要艇型參數(shù)。

    圖 1 DARPA suboff模型Fig. 1 Geometry of the DARPA suboff model

    表 1 DARPA suboff模型艇型參數(shù)Tab. 1 Principal dimensions of the suboff model

    2.2 計(jì)算設(shè)置與邊界條件

    本文分別對均勻流體(case A)及密度強(qiáng)分層流體(case B)中運(yùn)動(dòng)suboff尾跡特征進(jìn)行預(yù)報(bào)??紤]到流場的對稱性,為節(jié)約計(jì)算資源,提高計(jì)算效率,本文采用計(jì)算一半流場的方式進(jìn)行計(jì)算。計(jì)算域設(shè)置如圖2所示。考慮到suboff尾跡傳播距離遠(yuǎn)且持續(xù)時(shí)間久等特點(diǎn),為實(shí)現(xiàn)對suboff遠(yuǎn)場尾跡特征的預(yù)報(bào),計(jì)算域大小設(shè)置為:–100<x<20,0<y<40,–8<z<10。計(jì)算采用右手直角坐標(biāo)系,坐標(biāo)原點(diǎn)位于suboff艇首在靜止水面的投影點(diǎn),x軸沿艇長方向指向艇首,z軸垂直向上。本文計(jì)算中共采用7種邊界條件:速度入口邊界條件位于船首前方x=20m處,來流速度U=2.785 6 m/s;壓力出口邊界條件位于艇后方x=–100m處;y=0平面設(shè)置為對稱面邊界條件;y=40m平面設(shè)置為零梯度邊界條件;z=–8.0m平面設(shè)置為可滑移邊壁面界條件;z=10.0m平面設(shè)置為無窮遠(yuǎn)邊界條件;艇體表面設(shè)置為壁面邊界條件。

    圖 2 計(jì)算域設(shè)置Fig. 2 Computational domain and boundary conditions

    液體密度分層設(shè)置如圖3所示。其中,對于case A而言,ρ1=ρ2=997.561,而case B中液體密度設(shè)置為:ρ1:ρ2=997.561:1 030.0。對于 2 種工況,suboff均固定于上層液體中間位置。計(jì)算中,通過自定義函數(shù)設(shè)定不同密度流體范圍及流動(dòng)速度,實(shí)現(xiàn)對suboff勻速前行時(shí)流場的仿真。2種工況下具體計(jì)算設(shè)置見表2。

    圖 3 密度分層區(qū)域設(shè)置Fig. 3 Design of stratified fluid

    表 2 計(jì)算條件設(shè)置Tab. 2 Computational conditions

    2.3 網(wǎng)格

    為節(jié)省網(wǎng)格數(shù)量,提高計(jì)算效率,本計(jì)算在近suboff模型周圍采用非結(jié)構(gòu)化網(wǎng)格,在遠(yuǎn)離suboff位置采用結(jié)構(gòu)化網(wǎng)格。此外,為保證流場求解精度并更好地捕捉suboff尾跡特征,在艇體首尾、靜水面及流體密度分層處進(jìn)行了網(wǎng)格加密。case A與case B網(wǎng)格量分別為368萬和494萬。圖4和圖5分別給出了suboff表面網(wǎng)格及流體密度強(qiáng)分層工況計(jì)算域網(wǎng)格分布示意圖。

    圖 4 suboff表面網(wǎng)格Fig. 4 Surface mesh of the suboff model

    圖 5 強(qiáng)密度分層工況計(jì)算網(wǎng)格示意圖Fig. 5 Grid distributions for stratified fluid case

    2.4 結(jié)果分析

    1)網(wǎng)格收斂性驗(yàn)證

    本文基于均勻流體中suboff運(yùn)動(dòng)粘性流場數(shù)值試驗(yàn),根據(jù)網(wǎng)格在3個(gè)方向網(wǎng)格點(diǎn)布置依次根據(jù)倍變化的方式選取粗糙、中間及細(xì)密3種網(wǎng)格計(jì)算了計(jì)算收斂性驗(yàn)證,網(wǎng)格量分別為:130萬、368萬和1 040萬。圖6給出了3種網(wǎng)格計(jì)算得到的總阻力計(jì)算結(jié)果及其與實(shí)驗(yàn)結(jié)果的對比情況。圖7為3種網(wǎng)格工況下,艇體壓力沿船長方向分布計(jì)算結(jié)果及其與試驗(yàn)結(jié)果的對比情況。圖6計(jì)算結(jié)果表明,3種網(wǎng)格計(jì)算得到的總阻力誤差都在可接受范圍之內(nèi)。細(xì)密網(wǎng)格計(jì)算得到的總阻力系數(shù)更接近實(shí)驗(yàn)值,可能是由于細(xì)密的網(wǎng)格捕捉到船體周圍的不易捕捉的小渦,使得阻力計(jì)算更加準(zhǔn)確。圖7表明3種網(wǎng)格艇體表面壓力測量點(diǎn)計(jì)算結(jié)果都與試驗(yàn)值壓力較為吻合。此外,圖7結(jié)果還表明網(wǎng)格加密能夠更好地捕捉艇體壓力的突躍??紤]到計(jì)算資源限制,為兼顧計(jì)算精度與計(jì)算效率,本文將采用中間網(wǎng)格展開suboff尾跡的數(shù)值預(yù)報(bào)。

    圖 6 3種網(wǎng)格的總阻力計(jì)算結(jié)果與試驗(yàn)值對比圖Fig. 6 Comparison of the computed resistance and experimental data

    圖 7 3種網(wǎng)格艇體壓力計(jì)算結(jié)果與試驗(yàn)值對比值Fig. 7 Comparison of the computed pressure and experimental data

    2)阻力

    為研究流體分層對suboff水動(dòng)力的影響,本文對均勻流體(case A)及強(qiáng)密度分層流體(case B)工況下,suboff阻力性能進(jìn)行仿真計(jì)算,計(jì)算結(jié)果見表3。表3中結(jié)果表明B工況下suboff總阻力較A工況略高約0.77%??梢姡黧w分層對suboff所受阻力影響較小。這可能是由于B工況中,盡管流體存在強(qiáng)密度分層,但suboff仍位于密度較低液體層,流體分層引起的剪切作用對suboff周圍流場影響有限所致。

    表 3 不同工況下,suboff阻力性能計(jì)算結(jié)果Tab. 3 Computational resistance of two cases

    3)艇體壓力分布

    為進(jìn)一步研究A,B工況下,suboff阻力性能相近的原因,本文對2種工況下,suboff艇體中剖面壓力沿艇長方向的分布及suboff艇體表面壓力進(jìn)行計(jì)算與對比分析。由圖8和圖9中計(jì)算結(jié)果可以看出,2種工況下,suboff中剖面壓力存在3個(gè)峰值,分別為艇首、指揮臺前端及尾舵前端,且以前兩者最大。一方面,壓力突躍對相應(yīng)位置強(qiáng)度提出了要求,另一方面壓力突躍強(qiáng)行改變此處水流狀態(tài),不僅加大了suboff阻力,更不利于suboff隱身性。為進(jìn)一步提高suboff快速性與隱身性,西方軍事強(qiáng)國已深入開展低矮流線型指揮臺圍殼研制工作并已在最新一代suboff上取得成功應(yīng)用。同時(shí),計(jì)算結(jié)果還表明,2種工況下suboff艇體表面壓力差別微弱,進(jìn)一步說明2種工況下,suboff周圍流場相近,也從另一角度解釋了2種工況下suboff阻力性能相近的原因。

    圖 8 2種工況下,艇中剖面壓力分布Fig. 8 Computed pressure for two cases

    圖 9 2種工況下,壓力分布對比圖Fig. 9 Computed surface pressure for two cases

    4)速度分布

    圖10給出了2種工況下,x/L=0.978槳盤面處速度分布計(jì)算結(jié)果。圖中結(jié)果表明2種工況下suboff槳盤面位置水流速度顯著低于水流速度,且存在顯著的流動(dòng)分離。此外,圖10結(jié)果還表明2種工況下,suboff槳盤面伴流場位置水流速度差別微弱,進(jìn)一步說明B工況中流體分層對suboff周圍流場影響較弱。圖11為2種工況下,suboff中心位置所在水平面(z=–0.875 m)流場速度分布計(jì)算結(jié)果。從圖中可以看出:1)由于suboff的存在改變了其后方水流的分布,其水流流速呈現(xiàn)出顯著的開爾文特征。2)2種工況下,suboff尾跡橫波波長相同。3)2種工況下,suboff周圍水流速度相近,進(jìn)一步說明B工況下,流體分層對suboff周圍流場的改變較弱。4)2種工況下,遠(yuǎn)離suboff的位置suboff縱波尾跡流速相當(dāng)。但suboff尾跡流場橫波特征呈現(xiàn)出顯著不同。具體表現(xiàn)為B工況下,suboff尾跡流場橫波特征較A工況顯著。說明流體分層引起的內(nèi)波對suboff尾跡的影響主要表現(xiàn)為引起尾跡橫波特征加劇。

    圖 10 2種工況下,槳盤面伴流場速度分布對比圖Fig. 10 Wake at x/L=0.978 for two cases

    5)尾跡特征

    圖 11 2種工況下,z=–0.875截面速度對比圖Fig. 11 velocity at z=–0.875 m for two cases

    本文對2種工況下,suboff尾跡自由面特征進(jìn)行計(jì)算分析與對比。圖12為2種工況下suboff尾跡特征計(jì)算結(jié)果對比圖。結(jié)果表明:1)2種工況下,suboff尾跡特征均呈現(xiàn)顯著的開爾文波系特征;2)2種工況下suboff近場尾跡特征相近;3)2種工況下,suboff尾跡自艇后3個(gè)波長后尾跡特征開始呈現(xiàn)不同;4)B工況下,suboff遠(yuǎn)場尾跡中橫波作用較A工況顯著,且衰減速度緩慢。圖13給出了B工況下,液體分界層處的內(nèi)波尾跡特征計(jì)算結(jié)果。計(jì)算結(jié)果表明:1)內(nèi)波尾跡通常呈現(xiàn)開爾文波系特征,但其艇后縱波的傳播被約束在約5°范圍內(nèi);2)內(nèi)波尾跡中橫波波長與A工況時(shí)橫波波長一致;3)內(nèi)波尾跡縱波沿艇后方向持續(xù)傳播,波長可達(dá)幾十倍船長,因此內(nèi)波尾跡將有望和聲波一同成為能在水中遠(yuǎn)距離傳播信息的物理場,進(jìn)而成為核suboff暴露的主要物理場。圖12和圖13結(jié)果表明,在本文研究工況下,液體分層引起suboff自由面尾跡特征的變化主要集中在艇后3倍波長范圍之后,且主要表現(xiàn)為對suboff橫波尾跡的加劇。流體分層引起的suboff尾跡特征更為顯著,不利于suboff尾跡隱身性能。因此,在研究suboff尾跡隱身過程中,建議考慮液體分層引起的內(nèi)波尾跡對suboff自由面尾跡特征的影響。

    圖 12 2種工況下,自由面尾跡特征對比圖Fig. 12 Computed wave profiles for two cases

    圖 13 case B 內(nèi)波尾跡Fig. 13 Computed internal wave profiles for case B

    3 結(jié) 語

    本文以suboff模型為研究對象,基于非定常雷諾時(shí)均N-S方程(RANS),數(shù)值研究了均勻流體及強(qiáng)分層流體中運(yùn)動(dòng)suboff模型粘性流場及其尾跡特征。通過對2種工況下,suboff阻力、艇表壓力、尾跡流場速度分布及尾跡特征的分析與對比,直觀地給出了流體分層對suboff水動(dòng)力性能及其尾跡特征的影響規(guī)律。首先,以suboff模型在均勻流體中勻速前行工況進(jìn)行了數(shù)值試驗(yàn),計(jì)算結(jié)果與試驗(yàn)值吻合良好,說明本文采用的數(shù)值方法能夠正確地預(yù)報(bào)suboff水動(dòng)力性能。其次,通過網(wǎng)格收斂性驗(yàn)證,對計(jì)算結(jié)果的網(wǎng)格無關(guān)性進(jìn)行研究。最后,進(jìn)行suboff模型在強(qiáng)分層流體中勻速前行的數(shù)值試驗(yàn),并通過對比均勻流體及強(qiáng)分層流體工況下suboff水動(dòng)力性能及其尾跡特征,研究了內(nèi)波對suboff流場及其尾跡特征的影響。就本文研究工況而言,計(jì)算結(jié)果表明:

    1)2種工況下,suboff水面尾跡都呈現(xiàn)顯著的開爾文波系;

    2)2種工況下,suboff水動(dòng)力性能差異較小,即流體分層引起的內(nèi)波對近suboff周圍流場的影響微弱;

    3)2種工況下,suboff遠(yuǎn)場尾跡特征差異較大,即流體分層引起的內(nèi)波對suboff遠(yuǎn)后方流場影響較為顯著;

    4)流體分層引起的suboff自由面尾跡改變主要表現(xiàn)對suboff橫波尾跡的加強(qiáng),包括波幅的增加與傳播距離的增加。密度分層引起內(nèi)波對潛體尾跡的加強(qiáng),使得潛體尾跡在SAR圖像中呈現(xiàn)出明暗相間的條紋,而這種影響將導(dǎo)致潛體的暴露機(jī)率大幅增加;

    5)流體分層引起的內(nèi)波尾跡也呈現(xiàn)類似開爾文波系特征,但其縱波波長可達(dá)幾十倍船長,明顯大于橫波波長。即內(nèi)波尾跡具有持續(xù)時(shí)間長,傳播距離遠(yuǎn)等特點(diǎn),因此內(nèi)波尾跡有望和聲波一同成為能在水中遠(yuǎn)距離傳播信息的物理場,進(jìn)而成為核suboff暴露的主要物理場。

    基于上述結(jié)果,作者認(rèn)為流體分層對潛體尾跡特征具有加強(qiáng)作用,不利于其隱身。因此,為準(zhǔn)確研究潛體尾跡特征,保證潛體具備較高的尾跡隱身能力,建議考慮流體分層對潛體尾跡的影響,以加強(qiáng)潛體的安全性。此外,內(nèi)波尾跡特征具備超遠(yuǎn)距離傳播的能力,內(nèi)波尾跡探潛必將成為一種探潛新技術(shù)。相關(guān)資料表明,國外正在開展水下直接探測潛體內(nèi)波尾跡的研究。積極深入研究內(nèi)波尾跡產(chǎn)生機(jī)理及其抑制措施對潛體隱身性能提高具有重要意義。

    猜你喜歡
    內(nèi)波尾跡計(jì)算結(jié)果
    孤立內(nèi)波對過渡海域聲場干涉結(jié)構(gòu)的影響分析
    一種基于Radon 變換和尾跡模型的尾跡檢測算法
    內(nèi)波與死水,連潛艇都怕的海浪
    基于MODIS 遙感影像的安達(dá)曼海內(nèi)波特征參數(shù)分布及生成周期研究
    不等高軟橫跨橫向承力索計(jì)算及計(jì)算結(jié)果判斷研究
    甘肅科技(2020年20期)2020-04-13 00:30:40
    基于EEMD-Hilbert譜的渦街流量計(jì)尾跡振蕩特性
    SAR極化特征圖像與σ0圖像的海洋內(nèi)波可視性對比
    超壓測試方法對炸藥TNT當(dāng)量計(jì)算結(jié)果的影響
    噪聲對介質(zhì)損耗角正切計(jì)算結(jié)果的影響
    基于FABEMD和Goldstein濾波器的SAR艦船尾跡圖像增強(qiáng)方法
    国产av又大| 精品少妇一区二区三区视频日本电影| 亚洲中文av在线| 精品少妇一区二区三区视频日本电影| 99热6这里只有精品| av有码第一页| 欧美在线黄色| 在线十欧美十亚洲十日本专区| 18禁黄网站禁片午夜丰满| 人人妻,人人澡人人爽秒播| 亚洲国产中文字幕在线视频| 免费看美女性在线毛片视频| 国产男靠女视频免费网站| 亚洲在线自拍视频| 亚洲黑人精品在线| 国产成+人综合+亚洲专区| 中国美女看黄片| 精品福利观看| 在线永久观看黄色视频| 国产在线精品亚洲第一网站| 午夜福利在线观看吧| 国产91精品成人一区二区三区| 变态另类成人亚洲欧美熟女| 欧美另类亚洲清纯唯美| 国产亚洲精品久久久久5区| 亚洲自拍偷在线| 黑人操中国人逼视频| 亚洲一区高清亚洲精品| 国产精品av久久久久免费| 亚洲全国av大片| 波多野结衣巨乳人妻| 国产亚洲av嫩草精品影院| 两个人免费观看高清视频| 黄色成人免费大全| svipshipincom国产片| 日本 欧美在线| 婷婷六月久久综合丁香| 午夜两性在线视频| 高潮久久久久久久久久久不卡| 麻豆成人午夜福利视频| 久久久久久久久久黄片| 一级a爱视频在线免费观看| 亚洲国产中文字幕在线视频| 每晚都被弄得嗷嗷叫到高潮| 夜夜躁狠狠躁天天躁| 欧美黑人巨大hd| 免费在线观看成人毛片| 欧美激情极品国产一区二区三区| 国产精品影院久久| 亚洲最大成人中文| 精品第一国产精品| 女警被强在线播放| 90打野战视频偷拍视频| 日韩 欧美 亚洲 中文字幕| 久久久久久免费高清国产稀缺| 中文字幕av电影在线播放| 非洲黑人性xxxx精品又粗又长| 少妇裸体淫交视频免费看高清 | 国产亚洲欧美在线一区二区| 国产野战对白在线观看| 久久久久国产精品人妻aⅴ院| 国产久久久一区二区三区| 午夜福利在线观看吧| 国产三级黄色录像| 亚洲成人精品中文字幕电影| 手机成人av网站| 日本黄色视频三级网站网址| 久久久久久久久免费视频了| 久久九九热精品免费| 亚洲九九香蕉| 亚洲久久久国产精品| 国产野战对白在线观看| 男人舔女人的私密视频| 在线十欧美十亚洲十日本专区| 91大片在线观看| 免费在线观看亚洲国产| 国产精品久久视频播放| 美女高潮到喷水免费观看| 波多野结衣av一区二区av| 国产精品一区二区精品视频观看| xxx96com| 欧美丝袜亚洲另类 | 999久久久精品免费观看国产| 日韩欧美国产在线观看| 成人手机av| 久久久精品欧美日韩精品| 欧美黑人精品巨大| 国产视频一区二区在线看| 国产激情偷乱视频一区二区| 丝袜美腿诱惑在线| 色综合亚洲欧美另类图片| 每晚都被弄得嗷嗷叫到高潮| 变态另类成人亚洲欧美熟女| 一本大道久久a久久精品| 久久久久久九九精品二区国产 | 免费在线观看日本一区| 99热这里只有精品一区 | 国产人伦9x9x在线观看| 桃色一区二区三区在线观看| 一进一出抽搐动态| 欧美乱妇无乱码| 一边摸一边抽搐一进一小说| 日韩欧美免费精品| 男女做爰动态图高潮gif福利片| 久久久国产欧美日韩av| 丰满的人妻完整版| 黑丝袜美女国产一区| 午夜精品在线福利| 首页视频小说图片口味搜索| 精品午夜福利视频在线观看一区| tocl精华| 国产成人影院久久av| 亚洲黑人精品在线| 国产成人精品无人区| 亚洲av日韩精品久久久久久密| 久久久久久久久久黄片| 一进一出好大好爽视频| 一a级毛片在线观看| 2021天堂中文幕一二区在线观 | 国产黄色小视频在线观看| 免费高清在线观看日韩| 精品国产乱码久久久久久男人| av欧美777| 老汉色∧v一级毛片| 欧美不卡视频在线免费观看 | 一本久久中文字幕| 亚洲精品美女久久久久99蜜臀| 男人舔奶头视频| 少妇 在线观看| 国产亚洲精品av在线| 精品一区二区三区av网在线观看| 国产精品亚洲av一区麻豆| 欧美成狂野欧美在线观看| 男女之事视频高清在线观看| 免费在线观看完整版高清| 18禁观看日本| 老司机在亚洲福利影院| 久久 成人 亚洲| svipshipincom国产片| 最新在线观看一区二区三区| 1024视频免费在线观看| 亚洲人成77777在线视频| 色综合亚洲欧美另类图片| 久久精品人妻少妇| 日韩精品中文字幕看吧| 最新美女视频免费是黄的| av片东京热男人的天堂| 国产亚洲av嫩草精品影院| 亚洲无线在线观看| 亚洲精品国产精品久久久不卡| 中文字幕另类日韩欧美亚洲嫩草| 人成视频在线观看免费观看| 极品教师在线免费播放| 国产亚洲欧美在线一区二区| 亚洲男人天堂网一区| 亚洲 国产 在线| 免费在线观看影片大全网站| 亚洲va日本ⅴa欧美va伊人久久| 午夜日韩欧美国产| 我的亚洲天堂| 91老司机精品| 免费观看精品视频网站| 亚洲人成网站高清观看| 欧美大码av| 午夜久久久久精精品| 日本在线视频免费播放| 精品久久久久久成人av| 特大巨黑吊av在线直播 | 在线观看66精品国产| 天天一区二区日本电影三级| 超碰成人久久| 国产精品亚洲美女久久久| 亚洲国产日韩欧美精品在线观看 | 99久久久亚洲精品蜜臀av| 国产视频一区二区在线看| 国产视频一区二区在线看| 少妇 在线观看| 国产亚洲精品久久久久久毛片| 麻豆成人av在线观看| 久久精品91蜜桃| 夜夜看夜夜爽夜夜摸| 国产成人欧美| 丝袜美腿诱惑在线| 精品久久久久久久毛片微露脸| 日韩视频一区二区在线观看| 日韩视频一区二区在线观看| 亚洲九九香蕉| 18禁美女被吸乳视频| 国产精品自产拍在线观看55亚洲| 亚洲性夜色夜夜综合| 一夜夜www| 啪啪无遮挡十八禁网站| 国产成人精品无人区| 1024视频免费在线观看| 青草久久国产| 18禁裸乳无遮挡免费网站照片 | 国产成人一区二区三区免费视频网站| 国产极品粉嫩免费观看在线| 国内精品久久久久久久电影| aaaaa片日本免费| 中出人妻视频一区二区| 国产精品 欧美亚洲| 亚洲天堂国产精品一区在线| 夜夜躁狠狠躁天天躁| 成人精品一区二区免费| 亚洲av片天天在线观看| 亚洲一区二区三区不卡视频| 色婷婷久久久亚洲欧美| 欧美成人午夜精品| 久久婷婷成人综合色麻豆| 伊人久久大香线蕉亚洲五| 国产又爽黄色视频| 午夜免费观看网址| 成人免费观看视频高清| 日本a在线网址| 在线观看午夜福利视频| 中文字幕人妻熟女乱码| 两人在一起打扑克的视频| 欧美成人免费av一区二区三区| 久久久久精品国产欧美久久久| 亚洲人成网站高清观看| 91在线观看av| 欧美精品亚洲一区二区| 久久精品91蜜桃| 嫩草影视91久久| 国产又色又爽无遮挡免费看| 长腿黑丝高跟| 国产精品av久久久久免费| av超薄肉色丝袜交足视频| 香蕉国产在线看| 日韩一卡2卡3卡4卡2021年| 成人18禁在线播放| 成人欧美大片| 国产精品国产高清国产av| 在线观看66精品国产| 精品电影一区二区在线| 国产成人系列免费观看| 亚洲精品色激情综合| 极品教师在线免费播放| 中亚洲国语对白在线视频| 可以在线观看的亚洲视频| 人妻久久中文字幕网| 黑人欧美特级aaaaaa片| 黑人操中国人逼视频| 在线观看66精品国产| 亚洲全国av大片| 男人舔女人下体高潮全视频| 亚洲熟妇中文字幕五十中出| 国产精品av久久久久免费| 久久精品91无色码中文字幕| 99在线视频只有这里精品首页| 亚洲天堂国产精品一区在线| 欧美精品亚洲一区二区| 国产爱豆传媒在线观看 | 国产视频一区二区在线看| 久久婷婷成人综合色麻豆| 波多野结衣av一区二区av| 丰满的人妻完整版| 黄色视频不卡| 日本a在线网址| 青草久久国产| 久久久国产成人精品二区| 国产久久久一区二区三区| 免费在线观看日本一区| 搡老熟女国产l中国老女人| 欧美日韩黄片免| 在线视频色国产色| 宅男免费午夜| 久久久国产精品麻豆| 日本免费a在线| 久久性视频一级片| e午夜精品久久久久久久| 日韩欧美在线二视频| 精品久久久久久久人妻蜜臀av| 亚洲自偷自拍图片 自拍| 亚洲国产高清在线一区二区三 | 亚洲精品av麻豆狂野| 国产成人系列免费观看| 最新在线观看一区二区三区| 国产精品亚洲一级av第二区| a级毛片在线看网站| а√天堂www在线а√下载| 久久久久久久精品吃奶| 99在线视频只有这里精品首页| 老司机福利观看| 女人被狂操c到高潮| 久久久久国产一级毛片高清牌| 日本一本二区三区精品| 欧美久久黑人一区二区| 一本综合久久免费| 国产亚洲精品久久久久久毛片| 此物有八面人人有两片| 亚洲片人在线观看| 两个人免费观看高清视频| 久久久久久久午夜电影| 嫁个100分男人电影在线观看| 三级毛片av免费| 精品久久久久久久久久久久久 | 香蕉久久夜色| 亚洲熟女毛片儿| 午夜福利在线在线| 久久99热这里只有精品18| 亚洲av成人av| 国产成人精品无人区| 国产国语露脸激情在线看| 成熟少妇高潮喷水视频| 中文字幕av电影在线播放| 欧美精品亚洲一区二区| 男人舔奶头视频| 国产欧美日韩一区二区三| 欧美黑人欧美精品刺激| 欧美成人性av电影在线观看| 久久久国产成人免费| 亚洲男人的天堂狠狠| 99精品在免费线老司机午夜| 一边摸一边做爽爽视频免费| 亚洲自拍偷在线| 一二三四社区在线视频社区8| 色播亚洲综合网| 国内精品久久久久久久电影| 国产成人av激情在线播放| 国产精品1区2区在线观看.| 18禁观看日本| 深夜精品福利| 老鸭窝网址在线观看| 久久精品国产99精品国产亚洲性色| 人人澡人人妻人| 老司机在亚洲福利影院| 波多野结衣av一区二区av| 欧美性猛交黑人性爽| 后天国语完整版免费观看| 中文在线观看免费www的网站 | 国产又黄又爽又无遮挡在线| 国产伦在线观看视频一区| 日本免费一区二区三区高清不卡| 久久久久久久久久黄片| 久久天躁狠狠躁夜夜2o2o| 精品欧美一区二区三区在线| 免费一级毛片在线播放高清视频| 欧美日韩黄片免| 免费在线观看完整版高清| 99riav亚洲国产免费| www.精华液| 成年版毛片免费区| 中文字幕人成人乱码亚洲影| 免费看美女性在线毛片视频| 美女 人体艺术 gogo| 看免费av毛片| 999久久久精品免费观看国产| 欧美午夜高清在线| 母亲3免费完整高清在线观看| 欧美国产日韩亚洲一区| 一本大道久久a久久精品| 男女午夜视频在线观看| 亚洲熟妇熟女久久| 好男人在线观看高清免费视频 | 中文字幕人成人乱码亚洲影| 搡老妇女老女人老熟妇| 久久伊人香网站| 亚洲精品中文字幕在线视频| 久久天堂一区二区三区四区| 啪啪无遮挡十八禁网站| 夜夜夜夜夜久久久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲真实伦在线观看| 亚洲人成77777在线视频| bbb黄色大片| 99国产精品99久久久久| 欧美又色又爽又黄视频| 亚洲无线在线观看| 非洲黑人性xxxx精品又粗又长| 国产1区2区3区精品| 亚洲全国av大片| 19禁男女啪啪无遮挡网站| 青草久久国产| 久久国产精品人妻蜜桃| 亚洲五月色婷婷综合| 少妇被粗大的猛进出69影院| 丝袜人妻中文字幕| 高清在线国产一区| 777久久人妻少妇嫩草av网站| 男女做爰动态图高潮gif福利片| 男女视频在线观看网站免费 | 一本久久中文字幕| 制服人妻中文乱码| 亚洲 欧美 日韩 在线 免费| 亚洲精品久久成人aⅴ小说| 国产精品 国内视频| 精品日产1卡2卡| 制服人妻中文乱码| 熟女电影av网| 久久久久久人人人人人| 少妇被粗大的猛进出69影院| 色播亚洲综合网| 久久久久精品国产欧美久久久| 亚洲五月色婷婷综合| 88av欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 怎么达到女性高潮| 日韩精品青青久久久久久| 色综合站精品国产| 亚洲人成伊人成综合网2020| 亚洲国产日韩欧美精品在线观看 | 欧美成人性av电影在线观看| 欧美黄色淫秽网站| 超碰成人久久| 久久精品夜夜夜夜夜久久蜜豆 | 嫩草影视91久久| 国产精品一区二区三区四区久久 | 亚洲男人的天堂狠狠| 亚洲精品一卡2卡三卡4卡5卡| 国产91精品成人一区二区三区| 亚洲成av人片免费观看| 成人一区二区视频在线观看| 男人舔女人的私密视频| 丝袜美腿诱惑在线| 色综合欧美亚洲国产小说| 曰老女人黄片| 日韩欧美一区二区三区在线观看| 青草久久国产| 波多野结衣巨乳人妻| 18禁黄网站禁片免费观看直播| 国产成人精品无人区| 国产av一区在线观看免费| 侵犯人妻中文字幕一二三四区| 中文字幕人妻丝袜一区二区| 黄片大片在线免费观看| 国产亚洲欧美精品永久| 亚洲一区中文字幕在线| 国产精品影院久久| 一夜夜www| 国产精品永久免费网站| 精品国产超薄肉色丝袜足j| 久久久久久大精品| 91麻豆av在线| 国产精品免费视频内射| 欧美黄色片欧美黄色片| 精品无人区乱码1区二区| 欧美性长视频在线观看| 男人舔女人下体高潮全视频| 又大又爽又粗| 亚洲 欧美一区二区三区| 99国产极品粉嫩在线观看| av视频在线观看入口| 亚洲欧洲精品一区二区精品久久久| 亚洲成人久久性| 国产精品av久久久久免费| 亚洲欧美精品综合久久99| 久久精品91蜜桃| 久久久国产精品麻豆| 特大巨黑吊av在线直播 | 极品教师在线免费播放| 欧美日韩亚洲国产一区二区在线观看| 中文亚洲av片在线观看爽| 亚洲av日韩精品久久久久久密| 热99re8久久精品国产| 亚洲一区中文字幕在线| www.999成人在线观看| 人人妻人人澡人人看| 激情在线观看视频在线高清| 国产精品久久久久久精品电影 | 国内精品久久久久精免费| av天堂在线播放| 国产黄片美女视频| 午夜久久久久精精品| 亚洲三区欧美一区| 听说在线观看完整版免费高清| 精品久久久久久成人av| 久久久久久免费高清国产稀缺| 长腿黑丝高跟| 亚洲欧美精品综合久久99| 人人妻,人人澡人人爽秒播| 欧美成人一区二区免费高清观看 | 丰满的人妻完整版| 真人一进一出gif抽搐免费| 国产亚洲精品综合一区在线观看 | 精品一区二区三区视频在线观看免费| 久久久久久国产a免费观看| 成年女人毛片免费观看观看9| 亚洲专区国产一区二区| 一二三四社区在线视频社区8| 免费在线观看黄色视频的| netflix在线观看网站| 久久婷婷人人爽人人干人人爱| 给我免费播放毛片高清在线观看| 韩国av一区二区三区四区| 亚洲精品久久成人aⅴ小说| 男人操女人黄网站| 日韩国内少妇激情av| 黄片小视频在线播放| 日本撒尿小便嘘嘘汇集6| 欧美激情久久久久久爽电影| 国产一区二区在线av高清观看| 校园春色视频在线观看| 日本免费一区二区三区高清不卡| 欧美黑人巨大hd| 精品卡一卡二卡四卡免费| 国产精品,欧美在线| 少妇的丰满在线观看| 不卡一级毛片| 亚洲av熟女| 人人澡人人妻人| 女人高潮潮喷娇喘18禁视频| 欧美成狂野欧美在线观看| 久久久久国产精品人妻aⅴ院| 亚洲精品国产精品久久久不卡| svipshipincom国产片| 成人国产一区最新在线观看| 俺也久久电影网| netflix在线观看网站| 色综合亚洲欧美另类图片| av超薄肉色丝袜交足视频| 国产精品野战在线观看| 精品国产美女av久久久久小说| 亚洲一区高清亚洲精品| 天天一区二区日本电影三级| 男人舔女人下体高潮全视频| 欧美人与性动交α欧美精品济南到| 亚洲国产精品合色在线| 在线免费观看的www视频| 亚洲五月婷婷丁香| 中文字幕精品免费在线观看视频| 在线观看一区二区三区| 欧美成人一区二区免费高清观看 | 两个人免费观看高清视频| 桃色一区二区三区在线观看| 亚洲 欧美 日韩 在线 免费| 亚洲精品中文字幕一二三四区| 国产精品电影一区二区三区| 国产高清激情床上av| 久久久久久亚洲精品国产蜜桃av| 丝袜在线中文字幕| 母亲3免费完整高清在线观看| 色播在线永久视频| 亚洲精品美女久久久久99蜜臀| 欧美中文综合在线视频| 久久久国产欧美日韩av| 国产真实乱freesex| 亚洲最大成人中文| 国产精品九九99| 国产高清videossex| 国产不卡一卡二| 精品一区二区三区四区五区乱码| 亚洲中文字幕一区二区三区有码在线看 | 日本成人三级电影网站| av有码第一页| 夜夜躁狠狠躁天天躁| 亚洲第一av免费看| 一本一本综合久久| 亚洲精品一卡2卡三卡4卡5卡| 午夜亚洲福利在线播放| 一级毛片精品| 精品乱码久久久久久99久播| 91麻豆av在线| 欧美又色又爽又黄视频| 国产色视频综合| 色婷婷久久久亚洲欧美| 亚洲av电影在线进入| 一级毛片女人18水好多| 亚洲,欧美精品.| 欧洲精品卡2卡3卡4卡5卡区| 最好的美女福利视频网| 99在线视频只有这里精品首页| 在线观看免费视频日本深夜| 亚洲欧美一区二区三区黑人| 99国产精品99久久久久| 哪里可以看免费的av片| 亚洲人成网站高清观看| 婷婷六月久久综合丁香| 国产区一区二久久| 精品电影一区二区在线| 香蕉国产在线看| 一区二区三区激情视频| 亚洲精品久久国产高清桃花| 国产精品国产高清国产av| 国产亚洲精品久久久久5区| 久久久久久久久中文| 国产黄a三级三级三级人| 最新美女视频免费是黄的| 两个人免费观看高清视频| 99久久精品国产亚洲精品| 国产爱豆传媒在线观看 | 久久中文字幕一级| 中文字幕最新亚洲高清| 少妇被粗大的猛进出69影院| 精品午夜福利视频在线观看一区| 国产亚洲精品第一综合不卡| xxxwww97欧美| 12—13女人毛片做爰片一| 欧美成狂野欧美在线观看| 久热爱精品视频在线9| 听说在线观看完整版免费高清| 亚洲国产精品久久男人天堂| 亚洲国产精品999在线| 一个人观看的视频www高清免费观看 | 欧美国产日韩亚洲一区| 欧美成人午夜精品| 一区二区三区精品91| 亚洲欧美精品综合一区二区三区| 免费女性裸体啪啪无遮挡网站| 99久久99久久久精品蜜桃| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久人妻蜜臀av| 波多野结衣高清作品| 亚洲欧美日韩高清在线视频| 亚洲成av片中文字幕在线观看| 最新美女视频免费是黄的| 亚洲成人国产一区在线观看| xxx96com| 一二三四社区在线视频社区8| 中文字幕久久专区| 丰满的人妻完整版| 观看免费一级毛片| 亚洲久久久国产精品| 啦啦啦韩国在线观看视频|