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

    返回艙垂直自由入水砰擊過程的數(shù)值模擬

    2015-05-10 05:47:04李少偉屈秋林劉沛清周丹杰
    中國空間科學(xué)技術(shù) 2015年5期
    關(guān)鍵詞:返回艙水面峰值

    李少偉屈秋林劉沛清周丹杰

    (1北京機電工程研究所,北京100074) (2北京航空航天大學(xué)航空科學(xué)與工程學(xué)院,北京100191)

    返回艙垂直自由入水砰擊過程的數(shù)值模擬

    李少偉1屈秋林2劉沛清2周丹杰1

    (1北京機電工程研究所,北京100074) (2北京航空航天大學(xué)航空科學(xué)與工程學(xué)院,北京100191)

    文章數(shù)值研究返回艙垂直自由入水砰擊過程中水氣流場與返回艙運動之間的動力學(xué)與運動學(xué)耦合問題。水氣兩相流的流動方程選為非定常雷諾平均Navier-Stokes (URANS)方程和Realizable k-ε湍流模型,返回艙的運動方程選為剛體的一維平動方程,水氣交界面的追蹤采用流體體積函數(shù)(VOF)方法,返回艙與水面之間的相對運動采用動網(wǎng)格技術(shù)實現(xiàn)。在計算方法得到試驗驗證的基礎(chǔ)上,數(shù)值研究了返回艙的不同觸水速度、質(zhì)量對入水過程中砰擊力、加速度等參數(shù)的影響規(guī)律;給出了入水過程中加速度峰值與觸水速度及質(zhì)量之間的定量關(guān)系式;發(fā)現(xiàn)了返回艙底部的入水砰擊壓強峰值發(fā)生在入水初期,且壓強峰值始終位于噴濺射流的根部。

    水動力學(xué);水上降落;砰擊載荷;流體體積函數(shù)模型;返回艙

    1 引言

    入水砰擊問題廣泛存在于船舶、水下兵器、海洋結(jié)構(gòu)等研究領(lǐng)域,在航空航天領(lǐng)域也會經(jīng)常遇到,例如水上飛機的水面降落、航天返回艙的水上降落回收、民用飛機的水上迫降[1]等。結(jié)構(gòu)體入水砰擊過程是一個水體流場、結(jié)構(gòu)體運動和結(jié)構(gòu)體變形三者相互耦合的復(fù)雜過程,其中水體受到結(jié)構(gòu)體的沖擊力,水面會發(fā)生大的變形,形成飛濺,同時結(jié)構(gòu)體受到水體的砰擊力,其運動參數(shù)變化和結(jié)構(gòu)變形都很大,特別是在入水初期砰擊載荷極大,可能使結(jié)構(gòu)體發(fā)生嚴(yán)重破壞。

    入水砰擊問題的研究方法主要分為:試驗法,分析法與數(shù)值模擬方法。早期的研究以試驗和分析法為主。文獻(xiàn)[2]和文獻(xiàn)[3]提出了二維楔形和圓柱形剛體入水的分析法,文獻(xiàn)[4]修正了文獻(xiàn)[3]的沾濕因子,建立了semi-Wagner法,2004年文獻(xiàn)[5]又在semi-Wagner法的基礎(chǔ)上加入了浮力的影響。這些分析方法簡單快速,常被應(yīng)用于工程實踐中。但是試驗與分析法難以精確給出入水過程中流場和結(jié)構(gòu)體變形的細(xì)節(jié)變化,不利于開展深入的機理研究,而數(shù)值模擬方法剛好彌補這些不足,文獻(xiàn)[6]對發(fā)生在1929—2003年間的航空航天結(jié)構(gòu)入水問題研究進(jìn)行了總結(jié),指出數(shù)值模擬方法已逐漸成為入水砰擊問題研究的一種重要方法。

    返回艙入水砰擊問題的研究始于美國的Mercury、Gemini和Apollo計劃[6-7],初期主要采用試驗與經(jīng)驗公式的方法,目前多采用數(shù)值模擬方法。文獻(xiàn)[8]對類“聯(lián)盟號”返回艙入水問題進(jìn)行了試驗研究和理論分析,獲得了返回艙入水過程中的加速度和表面三個特征點壓力隨時間的變化過程;文獻(xiàn)[9]和文獻(xiàn)[10]均采用有限元方法數(shù)值模擬返回艙的入水問題,重點研究了入水過程中加速度隨時間的變化規(guī)律;文獻(xiàn)[11]采用有限元方法數(shù)值模擬了返回艙入水的流固耦合過程,主要關(guān)注了返回艙底部中心點處沖擊壓強的時間歷程。

    上述研究主要關(guān)注返回艙入水過程中加速度的變化規(guī)律,但是對觸水速度和質(zhì)量等參數(shù)對其的影響規(guī)律關(guān)注較少,也很少給出詳盡的流場結(jié)構(gòu)從而分析加速度的變化原因。本文針對此問題,采用有限體積法,應(yīng)用動網(wǎng)格技術(shù)和多相流流體體積函數(shù)(VOF)方法數(shù)值模擬類“聯(lián)盟號”返回艙垂直自由入水砰擊過程,研究返回艙的觸水速度、質(zhì)量對入水砰擊過程中載荷的影響規(guī)律以及相應(yīng)的流場結(jié)構(gòu)的變情況。

    2 物理模型及網(wǎng)格劃分

    由于各種返回艙的主體幾何外形相似,本文選用類聯(lián)盟號返回艙模型[8]作為代表,開展入水問題數(shù)值研究,如圖1(a)所示。返回艙底面為半徑R=2.5 m的球面,艙體柱面最大半徑B=1.25 m,模型關(guān)于其中心軸對稱。柱坐標(biāo)系固連于大地,y軸沿返回艙軸線鉛垂向上(見圖1(b))。

    由于返回艙垂直入水過程中流場具有軸對稱性,因此計算中選取二維軸對稱模型進(jìn)行計算。計算區(qū)域的劃分如圖1(b)所示,上方為空氣域,下方為水域。為了保證計算精度,采用結(jié)構(gòu)網(wǎng)格,并對返回艙表面及水面附近的網(wǎng)格進(jìn)行了加密,返回艙壁面第一層網(wǎng)格高度為5×10-5m,網(wǎng)格總數(shù)為140萬。圖2(a)為整體的網(wǎng)格劃分情況,圖2(b)為返回艙壁面附近的網(wǎng)格劃分情況。

    圖1 計算模型描述Fig.1 Computational model description

    圖2 計算網(wǎng)格分布Fig.2 Grid in computation zone

    初始時刻,返回艙以一定的速度和高度釋放后垂直自由下落,觸水之前受重力和空氣阻力的作用,觸水后則受到重力、空氣阻力和水體作用力的共同作用。為了對入水過程中返回艙底部的沖擊壓強進(jìn)行詳細(xì)研究,在其底部設(shè)置壓強監(jiān)測點P1~P10,如圖3所示。P1位于底部中心,每個監(jiān)測點水平間距為0.12 m。

    圖3 壓強監(jiān)測點的布置Fig.3 Position of pressure sensors

    3 計算方法

    本文假設(shè)返回艙為剛體,采用動力學(xué)與運動學(xué)耦合的數(shù)值方法研究返回艙垂直入水砰擊問題。采用數(shù)值方法求解多相流體的非定常雷諾平均Navier-Stokes(URANS)方程,同時采用多相流VOF模型以及Realizable k-ε湍流模型[12]對URANS方程進(jìn)行封閉;對于返回艙的運動問題,數(shù)值求解剛體的一維平動方程;采用動網(wǎng)格技術(shù),模擬返回艙和水面之間的相對運動。計算過程中空氣為可壓流體,水為不可壓流體,忽略水的表面張力,同時計算中尚不考慮返回艙的變形。

    3.1 流動控制方程

    流動控制方程選為軸對稱的可壓縮的RANS方程,其包括:連續(xù)方程、軸向動量方程、徑向動量方程、能量方程。

    連續(xù)性方程:

    軸向動量方程:

    徑向動量方程:

    能量方程:

    3.2 VOF模型

    采用多相流VOF模型捕捉自由水面,通過求解第二相的連續(xù)方程式(5)和式(6),以及修正的HRIC(High Resolution Interface Capturing)方法,確定自由水面位置。

    式中 αq、ρq、vq分別為q相流體體積分?jǐn)?shù)、密度和運動速度,其中:角標(biāo)q=0,1分別表示空氣相與水相。

    3.3 湍流模型

    計算中,湍流模型選取Realizable k-ε模型,其湍動能k及其耗散率ε輸運方程為

    式中 湍動黏性μt=Cμρk2/ε,Cμ為與剪切、旋轉(zhuǎn)速率、湍流信息相關(guān)的參數(shù);σk和σε分別為湍流強度k與湍流耗散率ε的普朗特數(shù)。

    3.4 返回艙動力學(xué)方程

    通過以上方程求解流場信息,可計算出返回艙所受的合力,然后根據(jù)返回艙當(dāng)前步計算的受力和當(dāng)前步的運動速度確定下一步的速度,以確定返回艙的運動。

    式中 Frm為返回艙所受軸向合力;Vrm為返回艙的軸向速度;y為返回艙的軸向位置坐標(biāo);n為當(dāng)前計算步數(shù);m為返回艙質(zhì)量。

    4 計算結(jié)果與分析

    本文對不同觸水速度V0(V0為接觸水面時的瞬時速度)、質(zhì)量下返回艙垂直入水問題進(jìn)行了研究,得到了多種工況下的返回艙載荷、運動規(guī)律以及相應(yīng)的砰擊流場結(jié)構(gòu)。

    圖4 不同觸水速度下返回艙的加速度峰值Fig.4 Acceleration peaks versus different water entry velocities

    4.1 數(shù)值方法驗證

    通過數(shù)值計算,得到了質(zhì)量m=3000 kg不同觸水速度下返回艙的加速度峰值amax,計算值、試驗值[8]及理論模型結(jié)果的比較如圖4所示。由圖可以看出:von Karman理論基于動量守恒,但未能考慮水面涌起,使得沾濕面積小于真實情況,從而所得結(jié)果小于試驗值;Wagner在von Karman理論基礎(chǔ)上考慮了涌起效應(yīng),但是沒有考慮入水砰擊產(chǎn)生的噴濺射流會帶走部分能量,所以其結(jié)果大于試驗值;Miloh對Wagner理論的沾濕因子進(jìn)行修正,所得結(jié)果小于Wagner結(jié)果而略大于試驗值;Laverty在Miloh理論的基礎(chǔ)上考慮了浮力的影響,所得結(jié)果略大于Miloh結(jié)果;本文的數(shù)值計算結(jié)果與試驗結(jié)果符合得很好且精度明顯高于理論結(jié)果,如: V0=12 m/s時,加速度峰值的試驗均值為279 m/s2,本文CFD結(jié)果為275 m/s2,誤差1.4%, Miloh理論值為312 m/s2,誤差11.4%。

    4.2 觸水速度的影響

    初始時刻,質(zhì)量為3000 kg的返回艙距離水面0.1 m,以不同初始速度下落,接觸水面時的瞬時速度為V0,觸水后速度V隨著時間的變化歷程見圖5。觸水后返回艙主要受到水體的砰擊力與自身重力的作用,由于水體的砰擊力較大,返回艙將發(fā)生減速,觸水速度越大,水體砰擊力越大,返回艙減速越快,待到返回艙底部完全入水以后(比如:V0=10 m/s時,t=30 ms;V0=12 m/s時, t=25 ms),砰擊力明顯減小,返回艙以較小的減速度向下運動。由于本文主要研究入水砰擊問題,而砰擊的大載荷發(fā)生在入水初期,所以為了節(jié)約計算資源,只對返回艙入水初期過程進(jìn)行計算。

    圖6為質(zhì)量為3 000 kg不同觸水速度下返回艙著水過程中的加速度a變化。在較大觸水速度下,返回艙觸水后加速度迅速增大,達(dá)到峰值后逐漸減小,最終趨于平緩。其中曲線變化過程中較高頻率的震蕩是由于水面波動引起的。在返回艙著水過程中,觸水速度越大,其加速度峰值越大,即載荷峰值也越大,這將對航天器以及航天員產(chǎn)生極大的危害。

    圖5 返回艙速度隨時間變化歷程Fig.5 Velocities of recovery module with different initial impact velocity versus time

    圖6 返回艙加速度隨時間變化歷程Fig.6 Accelerations of recovery module with different initial impact velocity versus time

    4.3 質(zhì)量的影響

    為了研究返回艙不同裝載對入水過程的影響,本文選取質(zhì)量分別為3000kg、4500 kg、6000 kg的返回艙在觸水速度為10 m/s的條件下進(jìn)行著水過程模擬。圖7和圖8分別為不同質(zhì)量返回艙垂向砰擊力F與加速度的時間歷程。不同質(zhì)量返回艙觸水后,其垂向砰擊力均迅速增加達(dá)到峰值后再衰弱,由于質(zhì)量的影響,大質(zhì)量的返回艙入水砰擊力比小質(zhì)量的返回艙大,而大質(zhì)量的返回艙其加速度較小,這是由于質(zhì)量增加引起的砰擊力增量的貢獻(xiàn)小于質(zhì)量增加的貢獻(xiàn),因而加速度較小。

    圖7 不同質(zhì)量返回艙垂向砰擊力的時間歷程Fig.7 Force of recovery module with different weights versus time

    圖8 不同質(zhì)量返回艙加速度的時間歷程Fig.8 Acceleration of recovery module with different weights versus time

    4.4 加速度峰值擬合公式

    為了便于工程應(yīng)用,現(xiàn)對類“聯(lián)盟號”返回艙在入水過程中的加速度峰值隨觸水速度和質(zhì)量的變化曲線進(jìn)行擬合。

    由文獻(xiàn)[2-5]可知,加速度峰值與結(jié)構(gòu)的觸水速度及質(zhì)量有如下關(guān)系:

    amax=bm-2/3,其中b為待定常系數(shù)。所以本文采用該公式擬合類“聯(lián)盟號”返回艙入水加速度峰值amax與返回艙質(zhì)量m及觸水速度V0之間的關(guān)系:

    圖9給出了不同觸水速度及質(zhì)量的情況下返回艙的加速度峰值CFD結(jié)果與擬合結(jié)果的對比圖。其中圖9(a)為質(zhì)量3000 kg條件下,加速度峰值隨速度的變化情況,圖9(b)為速度10 m/s條件下,加速度峰值隨質(zhì)量的變化情況。由圖可以看出,擬合公式與CFD計算值偏差在m=3 000 kg、V0=10 m/s時最大,約為6.5%,其余點誤差均小于6%,可以用于類“聯(lián)盟號”返回艙不同觸水速度及質(zhì)量的情況下的加速度峰值快速預(yù)測。

    圖9 不同觸水速度及質(zhì)量下返回艙的加速度峰值Fig.9 Acceleration peaks versus water entry velocities and weights

    4.5 砰擊流場結(jié)構(gòu)

    通過有限體積法計算返回艙入水的多相流場,能夠詳盡地給出返回艙入水過程中的水面變化以及流場內(nèi)部信息。

    圖10給出了V0=10 m/s時,質(zhì)量為3000 kg的返回艙入水過程中的監(jiān)測點壓強p隨時間的變化曲線及壓強峰值連線。返回艙開始觸水后,底部中心的1號監(jiān)測點壓強迅速增大,達(dá)到峰值后迅速衰減;隨著入水深度的增加,沿著徑向監(jiān)測點的壓強相繼出現(xiàn)峰值并迅速衰減,且其壓強峰值幅度也迅速減小。

    圖11為返回艙自由入水過程中,艙底不同位置的壓強峰值pmax曲線。可以看出在入水過程中,返回艙底部中心位置承受巨大的壓強(約10.5 MPa),而遠(yuǎn)離中心位置,壓強迅速減小,所以在返回艙設(shè)計過程中應(yīng)該注意加強底部中心位置結(jié)構(gòu),以滿足強度要求。

    圖10 返回艙監(jiān)測點壓強的時間變化歷程Fig.10 Pressure of the tensors versus time

    圖11 返回艙入水過程中艙底各位置壓強峰值Fig.11 Pressure peaks of recovery module

    圖12為返回艙入水過程中不同時刻的流體壓力云圖,其中白色線為自由水面,上方為空氣,下方為水體。在入水初期,返回艙底部的壓力峰值始終位于噴濺射流的根部位置(見圖12(a)~(d)),這與文獻(xiàn)[13]所述一致;隨著入水深度的繼續(xù)增加,返回艙底部完全浸沒水中,壓力最大值遠(yuǎn)離自由水面附近,轉(zhuǎn)移至返回艙底部正下方(見圖12(e)),且此時返回艙底部壓力已經(jīng)很小,沒有明顯的峰值。

    圖12 返回艙周圍的壓力場分布圖Fig.12 Pressure contour of fluid around recovery module

    圖13給出了V0=10m/s、質(zhì)量為3000kg、t=17ms時返回艙周圍流線及流體的速度云圖(圖中量值單位為Pa)。結(jié)合入水過程中的壓力云圖(見圖12(c))可以看出,在入水過程中,返回艙觸水后擠壓水體運動,由于下方靜止水體的限制,水體向返回艙四周水面運動,在返回艙與水面交界附近產(chǎn)生了較大的壓力差,從而形成高速的水體射流,在射流的剪切作用下,附近氣流形成了明顯的漩渦。

    圖13 t=17ms時返回艙周圍流線及速度云圖Fig.13 Streamlines and velocity magnitude contour of the fluid around recovery module at t=17ms

    5 結(jié)束語

    本文采用有限體積法對返回艙垂直自由入水砰擊問題進(jìn)行數(shù)值模擬,研究返回艙的觸水速度和質(zhì)量對入水過程中力學(xué)性能的影響規(guī)律,得到了以下結(jié)論:

    1)相同質(zhì)量的返回艙,觸水速度越大砰擊力和加速度越大;相同觸水速度的返回艙,質(zhì)量越小砰擊力越小,但加速度越大。

    2)本文擬合的加速度峰值計算公式能夠較好地用于類“聯(lián)盟號”返回艙在不同觸水速度及質(zhì)量下入水砰擊過程中加速度峰值的快速預(yù)測。

    3)在入水過程中,由于返回艙擠壓,水體流線向外圍偏轉(zhuǎn),并在返回艙與自由水面交界處形成高速射流,且返回艙底部的壓力峰值位于射流根部。

    [1] GUO B D,LIU P Q,QU Q L,et al.Effect of pitch angle on the initial stage of a transport aircraft ditching[J].Chinese Journal of Aeronautics,2013,26(1):17-26.

    [2] VON KARMAN T.The impact of seaplane floats during landing[R].Technical Note 321,NACA, Washington,D.C.,1929.

    [3] WAGNER V H.Phenomena associated with impacts and sliding on liquid surfaces[R].NACA,1947.

    [4] MILOH T.On the initial-stage slamming of a rigid sphere in a vertical water entry[J].Applied Ocean Research,1991,13(1):43-48.

    [5] LAVERTY S M.Experimental hydrodynamics of spherical projectiles impacting on a free surface using high speed imaging techniques[D].Boston:Massachusetts Institute of Technology,2004.

    [6] SEDDON C M,MOATAMEDI M.Review of water entry with applications to aerospace structures[J]. International Journal of Impact Engineering,2006,32(7):1045-1067.

    [7] BENSON H E.Water impact of the Apollo spacecraft[J].Journal of Spacecraft and Rockets,1966,3(8): 1282-1284.

    [8] 宣建明,繆弋,程軍,等.返回艙水上沖擊特性的試驗研究與理論計算[J].水動力學(xué)研究與進(jìn)展,2000, 15(3):276-286.

    XUAN JIANMING,MIU YI,CHENG JUN,et al.Test and theoretical research on water landing characteristics of space capsule[J].Journal of Hydrodynamics,2000,15(3):276-286.

    [9] JOHN T W,KAREN H LYLE.Simulating space capsule water landing with explicit finite element method[R]. AIAA-2007-1779,2007.

    [10] 張?zhí)搼?杜匯良,馬春生.返回艙著水的數(shù)值模擬及沖擊特性分析 [J].清華大學(xué)學(xué)報:自然科學(xué)版, 2010,50(8):1297-1301.

    ZHANG XUHUAI,DU HUILIANG,MA CHUNSHENG.Water impact simulations and analyses of space capsule response characteristics[J].Journal of Tsinghua University(Science and Technology), 2010,50(8):1297-1301.

    [11] 溫坤,李書.返回艙入水沖擊數(shù)值建模與響應(yīng)分析[J].飛機設(shè)計,2007,27(6):19-23.

    WEN KUN,LI SHU.The numerical modeling and response analysis of water impact of a recovery module[J].Aircraft Dessign,2007,27(6):19-23.

    [12] 王福軍.計算流體動力學(xué)分析-CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2011:272.

    WANG FUJUN.Analysis of computational fluid dynamics-CFD theory and application[M].Beijing: Tsinghua University Press,2011:272.

    [13] PIERRE M,DAVID L T,GUILLAUME O,et al.SPH high-performance computing simulations of rigid solids impacting the free-surface of water[J].Journal of Hydraulic Research,2010,48:126-134.

    Numerical Simulation on Vertical Water Entry Impact of a Recovery Module

    LI Shaowei1QU Qiulin2LIU Peiqing2ZHOU Danjie1
    (1 Beijing Institute of Mechanical and Electrical Engineering,Beijing 100074)
    (2 School of Aeronautic Science and Engineering,Beihang University,Beijing 100191)

    The complex multiphase and fluid-structure interaction problem of the vertical water entry of a recovery module was numerically researched.The two-dimensional axisymmetrical unsteady Reynolds averaged Navier-Stokes(URANS)equation and Realizable k-εturbulent equations were solved by the SIMPLE algorithm.The water impact on the recovery module in single-degree-of-freedom motion was simulated by dynamic mesh method.The air-water interface was tracked using the volume of fluid (VOF)model.The employed numerical methods were firstly validated by comparing simulation results with experimental data.The effects of slamming velocity and weight on slamming force were investigated.A relationship formula of the acceleration peak to variable velocities and weights was obtained.The maximum shock pressure on the recovery module occurred in the early stage when entrying into the water and the pressure peak is always located in the splash root close to the recovery module.

    Hydrodynamics;Water landing;Impact overload;Volume of fluid model; Recovery module

    10.3780/j.issn.1000-758X.2015.05.006

    (編輯:楊嬋、范真真)

    2014-10-21。收修改稿日期:2015-05-06

    李少偉 1987年生,2013年獲北京航空航天大學(xué)流體力學(xué)專業(yè)碩士學(xué)位,工程師。研究方向為飛行器氣動布局設(shè)計。

    猜你喜歡
    返回艙水面峰值
    “四單”聯(lián)動打造適齡兒童隊前教育峰值體驗
    少先隊活動(2022年9期)2022-11-23 06:55:52
    “ 神舟十三號”返回艙安全著陸
    軍事文摘(2022年12期)2022-07-13 03:11:16
    水黽是怎樣浮在水面的
    創(chuàng)造足以亂真的水面反光
    爭奪水面光伏
    能源(2016年3期)2016-12-01 05:11:02
    實踐十號返回艙回家
    太空探索(2016年5期)2016-07-12 15:17:55
    多用途飛船縮比返回艙成功著陸
    太空探索(2016年7期)2016-07-10 12:10:15
    寬占空比峰值電流型準(zhǔn)PWM/PFM混合控制
    基于峰值反饋的電流型PFM控制方法
    一塊水面
    滇池(2014年5期)2014-05-29 07:30:25
    国产精品综合久久久久久久免费 | 亚洲avbb在线观看| av有码第一页| 亚洲av成人av| 欧美老熟妇乱子伦牲交| 99国产极品粉嫩在线观看| 黑人操中国人逼视频| 久久伊人香网站| 桃色一区二区三区在线观看| 日韩av在线大香蕉| 日本a在线网址| 日本a在线网址| а√天堂www在线а√下载| 美女高潮喷水抽搐中文字幕| 久久久久久人人人人人| 色哟哟哟哟哟哟| 亚洲专区中文字幕在线| 日韩中文字幕欧美一区二区| 757午夜福利合集在线观看| 搡老熟女国产l中国老女人| 国产免费男女视频| 亚洲情色 制服丝袜| 男人的好看免费观看在线视频 | 亚洲 国产 在线| 妹子高潮喷水视频| 亚洲五月色婷婷综合| 不卡av一区二区三区| 人人澡人人妻人| 亚洲精品国产色婷婷电影| 久久这里只有精品19| 日本a在线网址| 亚洲一码二码三码区别大吗| 午夜精品在线福利| 中文字幕色久视频| 国产单亲对白刺激| 亚洲成国产人片在线观看| 国产精品一区二区免费欧美| 9色porny在线观看| svipshipincom国产片| 精品无人区乱码1区二区| av片东京热男人的天堂| 夜夜爽天天搞| 嫩草影院精品99| 久久精品国产亚洲av香蕉五月| 18禁观看日本| 欧美成人午夜精品| 精品国产国语对白av| netflix在线观看网站| 亚洲精品久久午夜乱码| 久久国产精品人妻蜜桃| 久久影院123| 欧美另类亚洲清纯唯美| 久久香蕉激情| 51午夜福利影视在线观看| 亚洲久久久国产精品| 成人三级黄色视频| 欧美日韩黄片免| 欧美国产精品va在线观看不卡| 91老司机精品| 精品久久久久久成人av| 色婷婷av一区二区三区视频| 久久久久国内视频| 国产精品 国内视频| 亚洲七黄色美女视频| 久久香蕉国产精品| 国产精品久久久人人做人人爽| 村上凉子中文字幕在线| 国产亚洲精品第一综合不卡| 午夜两性在线视频| 极品人妻少妇av视频| 三级毛片av免费| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品sss在线观看 | 欧美+亚洲+日韩+国产| 欧美不卡视频在线免费观看 | 大型av网站在线播放| 亚洲精品国产精品久久久不卡| 亚洲中文av在线| 黑丝袜美女国产一区| 国产三级在线视频| 国产精品1区2区在线观看.| 国产aⅴ精品一区二区三区波| 亚洲九九香蕉| 国产一卡二卡三卡精品| 国产精品综合久久久久久久免费 | 一区二区三区精品91| 午夜精品久久久久久毛片777| 亚洲欧美激情综合另类| 性少妇av在线| 国产精品久久久久久人妻精品电影| 狠狠狠狠99中文字幕| 亚洲国产精品一区二区三区在线| 国产午夜精品久久久久久| 久久精品国产亚洲av高清一级| 亚洲一码二码三码区别大吗| 欧美日韩乱码在线| 久久久国产精品麻豆| 香蕉丝袜av| 午夜福利在线免费观看网站| 亚洲av成人av| 亚洲熟妇熟女久久| 亚洲性夜色夜夜综合| 啪啪无遮挡十八禁网站| 18美女黄网站色大片免费观看| 免费女性裸体啪啪无遮挡网站| 啦啦啦免费观看视频1| 12—13女人毛片做爰片一| 99久久人妻综合| 欧美不卡视频在线免费观看 | 中文字幕人妻丝袜制服| 国产精品久久久人人做人人爽| 国产av一区二区精品久久| 欧美黑人欧美精品刺激| 国内毛片毛片毛片毛片毛片| 男女下面插进去视频免费观看| 一级,二级,三级黄色视频| 视频在线观看一区二区三区| 亚洲欧美一区二区三区黑人| 午夜亚洲福利在线播放| 亚洲一区二区三区欧美精品| 满18在线观看网站| 国产有黄有色有爽视频| 黄色视频,在线免费观看| 18禁美女被吸乳视频| 国产精品av久久久久免费| 日本免费a在线| 中文字幕人妻丝袜制服| 极品教师在线免费播放| 欧美成人性av电影在线观看| 久久精品亚洲av国产电影网| 国产麻豆69| 精品午夜福利视频在线观看一区| 色综合欧美亚洲国产小说| 桃色一区二区三区在线观看| 99久久精品国产亚洲精品| 亚洲免费av在线视频| 色综合欧美亚洲国产小说| 夫妻午夜视频| 黄色毛片三级朝国网站| 欧美日韩瑟瑟在线播放| 搡老熟女国产l中国老女人| 美女福利国产在线| 一区二区三区国产精品乱码| 国产精品野战在线观看 | 亚洲美女黄片视频| 久久国产精品男人的天堂亚洲| 精品国产一区二区三区四区第35| 国产欧美日韩一区二区三| 亚洲av美国av| 91成年电影在线观看| 长腿黑丝高跟| 亚洲av成人av| 国产一区二区激情短视频| 日韩三级视频一区二区三区| 十八禁网站免费在线| 天堂√8在线中文| 19禁男女啪啪无遮挡网站| 久久香蕉国产精品| 精品国产乱码久久久久久男人| 免费高清在线观看日韩| 少妇裸体淫交视频免费看高清 | 亚洲伊人色综图| 国产男靠女视频免费网站| 国产高清视频在线播放一区| 成人三级黄色视频| 免费高清视频大片| 亚洲av第一区精品v没综合| 日韩av在线大香蕉| 18禁裸乳无遮挡免费网站照片 | 亚洲国产欧美网| 91在线观看av| 色播在线永久视频| 欧美乱码精品一区二区三区| 丝袜在线中文字幕| 免费在线观看影片大全网站| 久久狼人影院| videosex国产| 亚洲五月天丁香| 69精品国产乱码久久久| 亚洲人成伊人成综合网2020| 99久久综合精品五月天人人| 精品国产乱子伦一区二区三区| 手机成人av网站| 成人av一区二区三区在线看| 精品一品国产午夜福利视频| 国产一区二区在线av高清观看| 国产亚洲精品久久久久5区| 69av精品久久久久久| 大香蕉久久成人网| 黄色片一级片一级黄色片| 欧美丝袜亚洲另类 | 首页视频小说图片口味搜索| 国产高清videossex| 久久午夜综合久久蜜桃| 欧美乱码精品一区二区三区| 精品国产亚洲在线| 18禁观看日本| 村上凉子中文字幕在线| avwww免费| 一二三四社区在线视频社区8| 黄色视频,在线免费观看| av网站免费在线观看视频| 99riav亚洲国产免费| 国产国语露脸激情在线看| 色婷婷av一区二区三区视频| 亚洲欧美精品综合久久99| 午夜免费鲁丝| 午夜91福利影院| a级毛片黄视频| 丰满饥渴人妻一区二区三| 大型黄色视频在线免费观看| 村上凉子中文字幕在线| 亚洲精品久久成人aⅴ小说| 在线免费观看的www视频| 亚洲性夜色夜夜综合| 麻豆av在线久日| 亚洲成a人片在线一区二区| 两个人看的免费小视频| 色婷婷久久久亚洲欧美| 亚洲人成网站在线播放欧美日韩| 国产精品免费视频内射| 变态另类成人亚洲欧美熟女 | 女人爽到高潮嗷嗷叫在线视频| 午夜激情av网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲欧美在线一区二区| 一边摸一边抽搐一进一小说| x7x7x7水蜜桃| 人人妻人人添人人爽欧美一区卜| 亚洲精品一卡2卡三卡4卡5卡| 乱人伦中国视频| 男女床上黄色一级片免费看| 国产伦人伦偷精品视频| 午夜精品久久久久久毛片777| 久久草成人影院| 亚洲欧美日韩高清在线视频| 欧美不卡视频在线免费观看 | 亚洲全国av大片| 国产精品影院久久| 久久国产精品男人的天堂亚洲| 欧美黄色片欧美黄色片| 亚洲av熟女| 亚洲一区高清亚洲精品| 亚洲色图综合在线观看| 久热爱精品视频在线9| 成人三级黄色视频| 美女高潮到喷水免费观看| 麻豆av在线久日| 国产单亲对白刺激| 久久精品国产综合久久久| 亚洲中文字幕日韩| 国产精品影院久久| 欧美日本中文国产一区发布| 久久久久九九精品影院| 高清在线国产一区| 激情在线观看视频在线高清| 亚洲全国av大片| 久久精品国产亚洲av香蕉五月| 国产一区二区三区在线臀色熟女 | 搡老熟女国产l中国老女人| 亚洲专区中文字幕在线| 亚洲av美国av| 国产成人精品无人区| 午夜精品在线福利| 国产精品野战在线观看 | 日本a在线网址| 高清av免费在线| 国产免费男女视频| 国产精品香港三级国产av潘金莲| cao死你这个sao货| 夜夜躁狠狠躁天天躁| 精品一区二区三区视频在线观看免费 | 亚洲精品国产色婷婷电影| 成人永久免费在线观看视频| 亚洲国产欧美日韩在线播放| 婷婷六月久久综合丁香| 久久中文字幕一级| 啦啦啦 在线观看视频| 日韩欧美三级三区| 久久精品影院6| 国产国语露脸激情在线看| 黄色怎么调成土黄色| 午夜视频精品福利| 国产欧美日韩精品亚洲av| 午夜精品在线福利| 国产aⅴ精品一区二区三区波| 曰老女人黄片| 久久久久亚洲av毛片大全| 亚洲,欧美精品.| 精品免费久久久久久久清纯| 99热国产这里只有精品6| 午夜成年电影在线免费观看| 欧美 亚洲 国产 日韩一| 国产深夜福利视频在线观看| 91精品三级在线观看| 午夜精品在线福利| av天堂久久9| 亚洲欧美精品综合一区二区三区| avwww免费| av在线天堂中文字幕 | 美女 人体艺术 gogo| 丰满迷人的少妇在线观看| 妹子高潮喷水视频| 美女高潮喷水抽搐中文字幕| 国产av又大| 人妻丰满熟妇av一区二区三区| av免费在线观看网站| 亚洲国产毛片av蜜桃av| 亚洲国产精品一区二区三区在线| 中文字幕色久视频| 欧美激情久久久久久爽电影 | 动漫黄色视频在线观看| 国产亚洲精品第一综合不卡| 欧美激情极品国产一区二区三区| 国产欧美日韩一区二区精品| 久久青草综合色| 悠悠久久av| 国产精品九九99| 中文欧美无线码| 久久久久久亚洲精品国产蜜桃av| 国产精品野战在线观看 | 高清毛片免费观看视频网站 | 欧美激情久久久久久爽电影 | 久久狼人影院| 麻豆av在线久日| 午夜久久久在线观看| 曰老女人黄片| av欧美777| 国产精品 国内视频| 精品午夜福利视频在线观看一区| 日韩免费高清中文字幕av| 久久亚洲精品不卡| 最近最新中文字幕大全免费视频| www国产在线视频色| 欧美日韩av久久| 在线永久观看黄色视频| 在线观看免费日韩欧美大片| 91老司机精品| 亚洲 欧美 日韩 在线 免费| 亚洲精品美女久久久久99蜜臀| 大香蕉久久成人网| 我的亚洲天堂| 成人国产一区最新在线观看| 欧美日韩亚洲综合一区二区三区_| 久久久国产成人精品二区 | 国产精品久久久人人做人人爽| 亚洲av第一区精品v没综合| 精品国产一区二区久久| 真人做人爱边吃奶动态| 久久人人97超碰香蕉20202| 久久久久久久午夜电影 | 久久婷婷成人综合色麻豆| 亚洲久久久国产精品| 一区二区三区精品91| 国产成人av激情在线播放| 麻豆成人av在线观看| 一二三四社区在线视频社区8| 99在线视频只有这里精品首页| av天堂久久9| 亚洲欧美激情在线| 1024香蕉在线观看| www.www免费av| 制服人妻中文乱码| 中文亚洲av片在线观看爽| 国产一卡二卡三卡精品| www.www免费av| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩另类电影网站| 一进一出好大好爽视频| 亚洲精华国产精华精| 国产精品永久免费网站| 少妇 在线观看| 1024香蕉在线观看| 在线观看免费日韩欧美大片| 人人澡人人妻人| 不卡一级毛片| 久久久水蜜桃国产精品网| 精品熟女少妇八av免费久了| 丁香六月欧美| 日韩欧美国产一区二区入口| 精品久久久久久久毛片微露脸| 国产麻豆69| 黑人操中国人逼视频| 一级a爱片免费观看的视频| 日本a在线网址| 最新美女视频免费是黄的| 纯流量卡能插随身wifi吗| 一夜夜www| www.熟女人妻精品国产| 成人av一区二区三区在线看| 男女床上黄色一级片免费看| 成人黄色视频免费在线看| 麻豆一二三区av精品| 亚洲国产精品999在线| 天堂中文最新版在线下载| 免费av中文字幕在线| 麻豆久久精品国产亚洲av | 久久国产精品人妻蜜桃| 欧美日韩亚洲高清精品| 午夜福利免费观看在线| 亚洲国产中文字幕在线视频| 老汉色∧v一级毛片| 熟女少妇亚洲综合色aaa.| 一级a爱视频在线免费观看| 欧美激情 高清一区二区三区| 制服诱惑二区| 啦啦啦在线免费观看视频4| 99国产极品粉嫩在线观看| 黄频高清免费视频| av天堂久久9| 久久欧美精品欧美久久欧美| 精品福利永久在线观看| 99精品在免费线老司机午夜| 80岁老熟妇乱子伦牲交| x7x7x7水蜜桃| 精品福利观看| 亚洲七黄色美女视频| 亚洲av成人一区二区三| 亚洲熟妇中文字幕五十中出 | 国产主播在线观看一区二区| 精品熟女少妇八av免费久了| 黑人操中国人逼视频| 天堂√8在线中文| www.熟女人妻精品国产| 乱人伦中国视频| 色精品久久人妻99蜜桃| 亚洲欧美一区二区三区久久| 亚洲精品一二三| 欧美黄色淫秽网站| 久久精品91蜜桃| 一二三四社区在线视频社区8| 国产精品秋霞免费鲁丝片| 亚洲av五月六月丁香网| 久久久水蜜桃国产精品网| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看影片大全网站| 性色av乱码一区二区三区2| 级片在线观看| av天堂在线播放| 18禁裸乳无遮挡免费网站照片 | 日韩中文字幕欧美一区二区| 叶爱在线成人免费视频播放| 最新美女视频免费是黄的| 精品人妻1区二区| 麻豆av在线久日| 韩国av一区二区三区四区| 久久中文看片网| 亚洲成国产人片在线观看| 天堂俺去俺来也www色官网| 久久人妻熟女aⅴ| 精品国产乱子伦一区二区三区| 国产野战对白在线观看| 国产精品 国内视频| 国产在线精品亚洲第一网站| 亚洲美女黄片视频| 久久午夜亚洲精品久久| 色老头精品视频在线观看| 国产精品爽爽va在线观看网站 | 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲avbb在线观看| 最新在线观看一区二区三区| 亚洲专区国产一区二区| 国产三级黄色录像| 久久人人爽av亚洲精品天堂| 国产不卡一卡二| 夜夜夜夜夜久久久久| 亚洲精品在线美女| 人成视频在线观看免费观看| 午夜日韩欧美国产| 久久天堂一区二区三区四区| 老司机午夜十八禁免费视频| 亚洲精品成人av观看孕妇| 亚洲九九香蕉| 亚洲精品久久午夜乱码| 亚洲精品久久成人aⅴ小说| 国产精品亚洲av一区麻豆| 老司机靠b影院| 免费看十八禁软件| 在线十欧美十亚洲十日本专区| 欧美日本中文国产一区发布| 法律面前人人平等表现在哪些方面| www.自偷自拍.com| 久久久久国内视频| 高清欧美精品videossex| 国产黄a三级三级三级人| 最近最新中文字幕大全电影3 | a级片在线免费高清观看视频| 三级毛片av免费| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品美女久久av网站| 国产亚洲欧美精品永久| 国产亚洲av高清不卡| 18禁美女被吸乳视频| 久久 成人 亚洲| 亚洲精品一二三| 国产亚洲欧美精品永久| 精品人妻1区二区| 如日韩欧美国产精品一区二区三区| a级毛片黄视频| 国产av精品麻豆| 99精品在免费线老司机午夜| 80岁老熟妇乱子伦牲交| 性少妇av在线| 亚洲精品中文字幕在线视频| 一夜夜www| 黑人操中国人逼视频| 国产一区二区在线av高清观看| 黑人猛操日本美女一级片| 美女午夜性视频免费| 97超级碰碰碰精品色视频在线观看| 久久香蕉精品热| 久久99一区二区三区| 看片在线看免费视频| 国产精品一区二区三区四区久久 | 欧美人与性动交α欧美软件| 国产三级黄色录像| 亚洲欧美日韩高清在线视频| 久久中文看片网| 高清在线国产一区| 国产主播在线观看一区二区| 中文字幕色久视频| 嫩草影视91久久| 波多野结衣一区麻豆| 国产精品一区二区精品视频观看| 级片在线观看| 国产野战对白在线观看| 久99久视频精品免费| 久久久久久免费高清国产稀缺| 最近最新免费中文字幕在线| 欧美激情久久久久久爽电影 | 欧美黑人欧美精品刺激| 久久伊人香网站| 国产不卡一卡二| 欧美中文日本在线观看视频| 乱人伦中国视频| 国产欧美日韩综合在线一区二区| 亚洲人成网站在线播放欧美日韩| 免费搜索国产男女视频| 亚洲va日本ⅴa欧美va伊人久久| 黑丝袜美女国产一区| 黄色毛片三级朝国网站| 亚洲熟妇中文字幕五十中出 | 99国产极品粉嫩在线观看| 免费在线观看黄色视频的| 动漫黄色视频在线观看| 丝袜人妻中文字幕| 在线十欧美十亚洲十日本专区| 99精国产麻豆久久婷婷| 91精品国产国语对白视频| 天堂动漫精品| 国产日韩一区二区三区精品不卡| 亚洲片人在线观看| www国产在线视频色| 午夜免费成人在线视频| 久9热在线精品视频| 日本黄色视频三级网站网址| 中国美女看黄片| 欧美日韩精品网址| 久久99一区二区三区| 中文字幕精品免费在线观看视频| 久久久久久大精品| 看片在线看免费视频| 天堂√8在线中文| 亚洲在线自拍视频| 在线免费观看的www视频| 伦理电影免费视频| 80岁老熟妇乱子伦牲交| 亚洲精品国产一区二区精华液| 亚洲国产精品999在线| 最新在线观看一区二区三区| 国产精品影院久久| 人人妻人人澡人人看| 欧美亚洲日本最大视频资源| 国产精华一区二区三区| 国产精品亚洲一级av第二区| 电影成人av| 欧美精品亚洲一区二区| 免费搜索国产男女视频| a级毛片黄视频| 日韩av在线大香蕉| 老熟妇乱子伦视频在线观看| 久久99一区二区三区| 美女国产高潮福利片在线看| 久久99一区二区三区| 老熟妇乱子伦视频在线观看| 亚洲国产欧美一区二区综合| 一级黄色大片毛片| 黄片大片在线免费观看| 不卡一级毛片| 久99久视频精品免费| 99国产极品粉嫩在线观看| 免费不卡黄色视频| 涩涩av久久男人的天堂| 两人在一起打扑克的视频| 神马国产精品三级电影在线观看 | 欧美黄色淫秽网站| 在线观看一区二区三区激情| 高清av免费在线| 日韩欧美一区二区三区在线观看| www.熟女人妻精品国产| 亚洲精品一二三| 久久人人爽av亚洲精品天堂| 男女床上黄色一级片免费看| 国产激情欧美一区二区| 一个人免费在线观看的高清视频| 美国免费a级毛片| 国产精品久久久久久人妻精品电影| 午夜两性在线视频| 真人一进一出gif抽搐免费| 9热在线视频观看99| 老鸭窝网址在线观看| 成人av一区二区三区在线看| 99精品久久久久人妻精品| av天堂在线播放| 国产精品久久电影中文字幕|