• <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
    在线观看免费午夜福利视频| 俺也久久电影网| 999精品在线视频| 一本精品99久久精品77| 在线a可以看的网站| 一边摸一边抽搐一进一小说| 国产精品野战在线观看| 欧美最黄视频在线播放免费| 国产亚洲欧美98| 免费人成视频x8x8入口观看| 久久精品国产清高在天天线| 国产一级毛片七仙女欲春2| netflix在线观看网站| 黄色片一级片一级黄色片| 亚洲国产欧洲综合997久久,| 国产av在哪里看| 久久国产精品影院| 小蜜桃在线观看免费完整版高清| 99热6这里只有精品| 在线观看美女被高潮喷水网站 | 国内少妇人妻偷人精品xxx网站 | 亚洲精品色激情综合| 日韩欧美精品v在线| 欧美不卡视频在线免费观看| av黄色大香蕉| av女优亚洲男人天堂 | 高清在线国产一区| 午夜免费激情av| 国产精品av久久久久免费| 他把我摸到了高潮在线观看| 我的老师免费观看完整版| 国产精品国产高清国产av| 国产成人精品久久二区二区免费| 亚洲精品乱码久久久v下载方式 | 免费观看的影片在线观看| 免费观看人在逋| 他把我摸到了高潮在线观看| 999久久久国产精品视频| 亚洲成av人片在线播放无| 亚洲第一电影网av| 久9热在线精品视频| 岛国在线免费视频观看| 91av网站免费观看| 国产伦精品一区二区三区视频9 | 国产三级黄色录像| 成人高潮视频无遮挡免费网站| 男女视频在线观看网站免费| 成人特级黄色片久久久久久久| 精品熟女少妇八av免费久了| 午夜福利在线观看免费完整高清在 | 久久中文字幕一级| 亚洲 欧美一区二区三区| 99国产极品粉嫩在线观看| 变态另类丝袜制服| 久久精品国产综合久久久| 国产综合懂色| 亚洲av片天天在线观看| 久久久久国产精品人妻aⅴ院| 不卡一级毛片| 精品人妻1区二区| 最新中文字幕久久久久 | 在线观看日韩欧美| 亚洲国产精品久久男人天堂| 99久国产av精品| 成人三级做爰电影| 91av网站免费观看| 欧美日韩乱码在线| 国产一区二区激情短视频| 欧美黑人欧美精品刺激| 国产成年人精品一区二区| 亚洲18禁久久av| 国产真人三级小视频在线观看| 熟妇人妻久久中文字幕3abv| 精品一区二区三区av网在线观看| 一进一出抽搐动态| 狠狠狠狠99中文字幕| 18禁国产床啪视频网站| 身体一侧抽搐| 老鸭窝网址在线观看| 国产精品av视频在线免费观看| 欧美乱色亚洲激情| 欧美高清成人免费视频www| 国产乱人伦免费视频| 999精品在线视频| 国产野战对白在线观看| 亚洲成人精品中文字幕电影| 国产乱人视频| 搡老岳熟女国产| 亚洲片人在线观看| 午夜福利欧美成人| 国产主播在线观看一区二区| 日本 av在线| 久久精品aⅴ一区二区三区四区| 久久精品夜夜夜夜夜久久蜜豆| 18禁国产床啪视频网站| 久久中文看片网| 一级a爱片免费观看的视频| 午夜a级毛片| 精品国产三级普通话版| 最新中文字幕久久久久 | 国产精品综合久久久久久久免费| 激情在线观看视频在线高清| 国产成人精品久久二区二区免费| 国产又黄又爽又无遮挡在线| 日韩大尺度精品在线看网址| 淫妇啪啪啪对白视频| 日韩大尺度精品在线看网址| 国产精品电影一区二区三区| 精品国产乱子伦一区二区三区| 国产 一区 欧美 日韩| 国产aⅴ精品一区二区三区波| 国产精品 欧美亚洲| 中文在线观看免费www的网站| 色精品久久人妻99蜜桃| 亚洲国产精品成人综合色| 精华霜和精华液先用哪个| 成人国产综合亚洲| 99视频精品全部免费 在线 | 久久久水蜜桃国产精品网| 国产一级毛片七仙女欲春2| 久久伊人香网站| 蜜桃久久精品国产亚洲av| 欧美中文综合在线视频| 久久伊人香网站| 久久久久精品国产欧美久久久| 天堂网av新在线| 亚洲av成人一区二区三| 老司机在亚洲福利影院| 两性午夜刺激爽爽歪歪视频在线观看| 色综合婷婷激情| 亚洲av片天天在线观看| 好看av亚洲va欧美ⅴa在| 国产乱人伦免费视频| 国产亚洲av高清不卡| 亚洲 欧美 日韩 在线 免费| 成年女人看的毛片在线观看| 免费看a级黄色片| 欧美黑人欧美精品刺激| 日本免费a在线| 国产精品av久久久久免费| 久久久国产精品麻豆| 亚洲国产欧洲综合997久久,| 久久久久久久久久黄片| 久久精品国产亚洲av香蕉五月| 人人妻,人人澡人人爽秒播| 999精品在线视频| 欧美+亚洲+日韩+国产| 男女午夜视频在线观看| 亚洲精品一区av在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 伊人久久大香线蕉亚洲五| 不卡一级毛片| 国产精品久久久av美女十八| 女生性感内裤真人,穿戴方法视频| 欧美黑人巨大hd| 夜夜爽天天搞| 欧美zozozo另类| 日日干狠狠操夜夜爽| 国产成年人精品一区二区| 18禁国产床啪视频网站| 久久这里只有精品中国| 国产精品av视频在线免费观看| 美女大奶头视频| 不卡一级毛片| 又粗又爽又猛毛片免费看| 国产男靠女视频免费网站| 亚洲国产精品sss在线观看| 国产精品av视频在线免费观看| 黄片小视频在线播放| 亚洲av日韩精品久久久久久密| 老鸭窝网址在线观看| 别揉我奶头~嗯~啊~动态视频| 99久久精品国产亚洲精品| svipshipincom国产片| 9191精品国产免费久久| 久久久久国产精品人妻aⅴ院| 黄色 视频免费看| 动漫黄色视频在线观看| 久久亚洲精品不卡| 97超视频在线观看视频| 俄罗斯特黄特色一大片| 免费人成视频x8x8入口观看| 国产精品美女特级片免费视频播放器 | 一本久久中文字幕| 操出白浆在线播放| 日本黄色片子视频| 级片在线观看| 成人精品一区二区免费| 亚洲一区二区三区不卡视频| 久久精品亚洲精品国产色婷小说| 亚洲最大成人中文| 国内精品一区二区在线观看| 久久久久国产精品人妻aⅴ院| 一个人看的www免费观看视频| 午夜两性在线视频| 国产午夜福利久久久久久| 午夜两性在线视频| 麻豆久久精品国产亚洲av| 天天躁狠狠躁夜夜躁狠狠躁| 一进一出好大好爽视频| 一二三四社区在线视频社区8| 国产一区二区三区在线臀色熟女| 国产高清有码在线观看视频| tocl精华| 国产麻豆成人av免费视频| 国内精品久久久久久久电影| e午夜精品久久久久久久| 中出人妻视频一区二区| 精品久久久久久久久久免费视频| 中文字幕久久专区| 日本一本二区三区精品| 日本在线视频免费播放| 免费在线观看日本一区| 少妇熟女aⅴ在线视频| 999精品在线视频| 最新中文字幕久久久久 | 国产伦人伦偷精品视频| 免费人成视频x8x8入口观看| 男人和女人高潮做爰伦理| 精品久久久久久成人av| 亚洲欧美日韩高清在线视频| 白带黄色成豆腐渣| 午夜福利欧美成人| 午夜亚洲福利在线播放| 国产成人影院久久av| 性色avwww在线观看| 国产美女午夜福利| 色av中文字幕| 最近最新中文字幕大全电影3| 亚洲午夜理论影院| 午夜日韩欧美国产| av福利片在线观看| 久久欧美精品欧美久久欧美| 欧美黑人欧美精品刺激| av视频在线观看入口| 亚洲aⅴ乱码一区二区在线播放| 日韩高清综合在线| 国产欧美日韩一区二区精品| 老司机午夜福利在线观看视频| 丁香六月欧美| 三级男女做爰猛烈吃奶摸视频| 色在线成人网| 成年女人永久免费观看视频| 99久久无色码亚洲精品果冻| 精品一区二区三区四区五区乱码| 日本三级黄在线观看| 国产免费男女视频| 村上凉子中文字幕在线| 国产精品一区二区三区四区久久| 国产精品 国内视频| www国产在线视频色| 久久久久免费精品人妻一区二区| 亚洲成人久久性| 国产精品一及| 国产精品一区二区三区四区久久| 久久久久久久午夜电影| 国产成人精品久久二区二区91| 1024手机看黄色片| 婷婷丁香在线五月| 一级毛片女人18水好多| 熟妇人妻久久中文字幕3abv| 日本黄大片高清| 日韩成人在线观看一区二区三区| 色吧在线观看| 91av网站免费观看| 真人一进一出gif抽搐免费| 18禁黄网站禁片午夜丰满| 亚洲精品乱码久久久v下载方式 | 一卡2卡三卡四卡精品乱码亚洲| 日日夜夜操网爽| 国产精品一区二区精品视频观看| 伊人久久大香线蕉亚洲五| 草草在线视频免费看| 亚洲性夜色夜夜综合| 嫩草影院入口| 别揉我奶头~嗯~啊~动态视频| 亚洲国产色片| 丰满人妻熟妇乱又伦精品不卡| 岛国在线观看网站| 亚洲电影在线观看av| 亚洲色图av天堂| 又黄又粗又硬又大视频| 日韩欧美在线乱码| 日韩人妻高清精品专区| 久久久成人免费电影| 成在线人永久免费视频| 午夜激情欧美在线| 久久午夜综合久久蜜桃| 成年女人看的毛片在线观看| 亚洲自拍偷在线| 97超视频在线观看视频| 成人午夜高清在线视频| 无人区码免费观看不卡| 国产高清三级在线| 岛国在线观看网站| 神马国产精品三级电影在线观看| 国产精品免费一区二区三区在线| 一个人免费在线观看电影 | 国产精品乱码一区二三区的特点| ponron亚洲| 在线永久观看黄色视频| 叶爱在线成人免费视频播放| 最近在线观看免费完整版| 高清在线国产一区| 不卡av一区二区三区| 人人妻,人人澡人人爽秒播| 日本精品一区二区三区蜜桃| 色精品久久人妻99蜜桃| 精品日产1卡2卡| 国产精品乱码一区二三区的特点| 最好的美女福利视频网| 琪琪午夜伦伦电影理论片6080| 久久性视频一级片| 一本精品99久久精品77| 老司机在亚洲福利影院| 99re在线观看精品视频| 一夜夜www| 久久久国产成人精品二区| 9191精品国产免费久久| 日本免费一区二区三区高清不卡| 丰满的人妻完整版| 午夜免费激情av| 精品国产乱子伦一区二区三区| 亚洲专区中文字幕在线| 色综合站精品国产| 欧美色视频一区免费| 国产精品电影一区二区三区| 女人高潮潮喷娇喘18禁视频| 亚洲无线在线观看| 久99久视频精品免费| 真人做人爱边吃奶动态| 伊人久久大香线蕉亚洲五| 国产成人福利小说| 午夜精品在线福利| www.熟女人妻精品国产| 色老头精品视频在线观看| 亚洲欧美日韩无卡精品| 久久久久久国产a免费观看| 成年人黄色毛片网站| 在线播放国产精品三级| 国内精品一区二区在线观看| 免费看a级黄色片| 综合色av麻豆| bbb黄色大片| 天天一区二区日本电影三级| 国产不卡一卡二| 亚洲欧美日韩卡通动漫| 一卡2卡三卡四卡精品乱码亚洲| 精品福利观看| 国产高清视频在线观看网站| 久久久色成人| 可以在线观看毛片的网站| 女人高潮潮喷娇喘18禁视频| 国产精品精品国产色婷婷| 日本黄色片子视频| 亚洲一区二区三区不卡视频| 亚洲乱码一区二区免费版| 岛国在线免费视频观看| 成人18禁在线播放| av天堂中文字幕网| 99精品欧美一区二区三区四区| 亚洲人与动物交配视频| 国产成+人综合+亚洲专区| 欧美在线黄色| 国产精品久久久久久人妻精品电影| 欧美一级毛片孕妇| 亚洲自偷自拍图片 自拍| 午夜福利欧美成人| 久久性视频一级片| 亚洲 欧美一区二区三区| 国产av在哪里看| 国产成人精品久久二区二区91| 国产高清videossex| 亚洲熟女毛片儿| 视频区欧美日本亚洲| 特大巨黑吊av在线直播| 黄色成人免费大全| 亚洲国产欧洲综合997久久,| 女警被强在线播放| 99国产精品99久久久久| 五月伊人婷婷丁香| 日韩高清综合在线| 国产av麻豆久久久久久久| 亚洲自拍偷在线| 母亲3免费完整高清在线观看| 久久国产精品人妻蜜桃| 97碰自拍视频| 亚洲国产精品成人综合色| 日日摸夜夜添夜夜添小说| 激情在线观看视频在线高清| 午夜影院日韩av| 国产亚洲av嫩草精品影院| 午夜成年电影在线免费观看| 一级黄色大片毛片| 搡老岳熟女国产| 免费看美女性在线毛片视频| 亚洲第一电影网av| 757午夜福利合集在线观看| 亚洲18禁久久av| 十八禁人妻一区二区| 人人妻人人看人人澡| 少妇人妻一区二区三区视频| 18禁裸乳无遮挡免费网站照片| 免费一级毛片在线播放高清视频| 亚洲精品在线观看二区| 性色av乱码一区二区三区2| 九色国产91popny在线| 亚洲av成人不卡在线观看播放网| 嫩草影院精品99| 高清在线国产一区| 一进一出抽搐gif免费好疼| 国产精品女同一区二区软件 | 丰满人妻熟妇乱又伦精品不卡| 老汉色∧v一级毛片| 亚洲av日韩精品久久久久久密| 欧美日韩一级在线毛片| 91九色精品人成在线观看| 长腿黑丝高跟| 精品久久久久久成人av| www日本在线高清视频| www日本黄色视频网| 国产亚洲精品久久久com| 高清在线国产一区| 成人特级黄色片久久久久久久| 淫秽高清视频在线观看| 色尼玛亚洲综合影院| 免费在线观看影片大全网站| 国产精品av久久久久免费| 亚洲国产精品久久男人天堂| 国产一级毛片七仙女欲春2| 亚洲狠狠婷婷综合久久图片| 中出人妻视频一区二区| 熟女电影av网| 亚洲无线在线观看| 99热这里只有精品一区 | 欧美在线黄色| www.www免费av| 国产精品电影一区二区三区| 亚洲欧美日韩无卡精品| 麻豆国产av国片精品| 亚洲自偷自拍图片 自拍| 日韩欧美免费精品| 日韩av在线大香蕉| 亚洲 欧美一区二区三区| 亚洲欧美日韩无卡精品| 国产成人影院久久av| 少妇的逼水好多| 午夜福利18| 1024手机看黄色片| 老熟妇乱子伦视频在线观看| 国产av不卡久久| 成熟少妇高潮喷水视频| av在线蜜桃| 男女之事视频高清在线观看| 高潮久久久久久久久久久不卡| 搡老岳熟女国产| 一级作爱视频免费观看| 网址你懂的国产日韩在线| 女生性感内裤真人,穿戴方法视频| e午夜精品久久久久久久| 色老头精品视频在线观看| 88av欧美| 成人鲁丝片一二三区免费| 熟女电影av网| 91在线精品国自产拍蜜月 | 99热这里只有精品一区 | 亚洲av免费在线观看| 久久久久国产一级毛片高清牌| 12—13女人毛片做爰片一| 日韩av在线大香蕉| 看片在线看免费视频| 免费av毛片视频| 精品久久蜜臀av无| 悠悠久久av| 中亚洲国语对白在线视频| 国产欧美日韩精品一区二区| 九九久久精品国产亚洲av麻豆 | 国产又黄又爽又无遮挡在线| 国产野战对白在线观看| 国产精品自产拍在线观看55亚洲| 欧美中文综合在线视频| 老汉色∧v一级毛片| 久久精品综合一区二区三区| 欧美极品一区二区三区四区| 可以在线观看毛片的网站| 麻豆久久精品国产亚洲av| 日本精品一区二区三区蜜桃| 一夜夜www| 亚洲 欧美一区二区三区| 亚洲真实伦在线观看| 欧美成人性av电影在线观看| 久久久久精品国产欧美久久久| 12—13女人毛片做爰片一| 夜夜爽天天搞| 丁香六月欧美| 日韩欧美一区二区三区在线观看| 日韩有码中文字幕| 女同久久另类99精品国产91| 在线观看美女被高潮喷水网站 | 成人国产综合亚洲| 九色成人免费人妻av| 欧美国产日韩亚洲一区| 久久久国产欧美日韩av| 亚洲国产色片| 全区人妻精品视频| 日韩欧美免费精品| 精品国产乱子伦一区二区三区| 男人舔女人下体高潮全视频| x7x7x7水蜜桃| 国产99白浆流出| 这个男人来自地球电影免费观看| 在线十欧美十亚洲十日本专区| 国产黄a三级三级三级人| 中文字幕最新亚洲高清| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利视频1000在线观看| 757午夜福利合集在线观看| 色综合站精品国产| 日本一本二区三区精品| 国内揄拍国产精品人妻在线| 美女高潮的动态| 日韩人妻高清精品专区| 久久性视频一级片| 国产黄片美女视频| 香蕉丝袜av| 一级毛片高清免费大全| 国产精品免费一区二区三区在线| 欧美一级毛片孕妇| 婷婷精品国产亚洲av在线| 欧美乱妇无乱码| 一二三四在线观看免费中文在| 91av网站免费观看| 日韩欧美在线乱码| 亚洲18禁久久av| 午夜久久久久精精品| 国产伦精品一区二区三区视频9 | 国产精品爽爽va在线观看网站| 久久热在线av| 国产精品久久久av美女十八| 熟女少妇亚洲综合色aaa.| 久久精品aⅴ一区二区三区四区| 美女扒开内裤让男人捅视频| 亚洲乱码一区二区免费版| 一级作爱视频免费观看| 身体一侧抽搐| 高清毛片免费观看视频网站| 亚洲专区国产一区二区| 欧美zozozo另类| 欧美日韩国产亚洲二区| 免费看美女性在线毛片视频| 91在线观看av| 国产精品爽爽va在线观看网站| 看黄色毛片网站| 国产亚洲欧美在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 搡老岳熟女国产| 国产视频内射| 久久精品人妻少妇| 国产欧美日韩一区二区精品| 一区二区三区激情视频| 桃红色精品国产亚洲av| 亚洲中文字幕一区二区三区有码在线看 | 听说在线观看完整版免费高清| 成人精品一区二区免费| 日韩三级视频一区二区三区| 最新美女视频免费是黄的| 久久精品夜夜夜夜夜久久蜜豆| 女警被强在线播放| 国产精品综合久久久久久久免费| 禁无遮挡网站| 国产精品久久久久久久电影 | cao死你这个sao货| 亚洲午夜精品一区,二区,三区| 51午夜福利影视在线观看| 国产亚洲精品av在线| 国产三级黄色录像| 一级毛片精品| 亚洲精品粉嫩美女一区| 中文字幕久久专区| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 色吧在线观看| 欧美大码av| 午夜免费观看网址| 丰满人妻熟妇乱又伦精品不卡| 国产精品香港三级国产av潘金莲| 亚洲第一欧美日韩一区二区三区| 国产乱人视频| 美女 人体艺术 gogo| 国内少妇人妻偷人精品xxx网站 | 日本熟妇午夜| 岛国在线免费视频观看| 国产成+人综合+亚洲专区| 久久中文看片网| 国产69精品久久久久777片 | 精品乱码久久久久久99久播| x7x7x7水蜜桃| 小说图片视频综合网站| 久久久精品欧美日韩精品| 天堂影院成人在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲精品久久国产高清桃花| 欧美高清成人免费视频www| 国产蜜桃级精品一区二区三区| 黄色片一级片一级黄色片| 欧美日韩综合久久久久久 | 午夜福利视频1000在线观看| 国产伦在线观看视频一区| 香蕉av资源在线| 狂野欧美激情性xxxx| 国产高潮美女av| 久久久久国内视频| 欧美乱色亚洲激情| 国产亚洲欧美在线一区二区| 亚洲av成人一区二区三|