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

    球錐形頭罩高超聲速繞流流場數(shù)值模擬

    2013-11-04 03:05:26許東李欣徐力王亞輝
    飛行力學 2013年6期
    關(guān)鍵詞:頭罩迎角湍流

    許東, 李欣, 徐力, 王亞輝

    (1.北京航空航天大學 儀器科學與光電工程學院, 北京 100191; 2.北京航天自動控制研究所 宇航與智能控制技術(shù)國家級重點實驗室, 北京 100854)

    球錐形頭罩高超聲速繞流流場數(shù)值模擬

    許東1, 李欣1, 徐力2, 王亞輝2

    (1.北京航空航天大學 儀器科學與光電工程學院, 北京 100191; 2.北京航天自動控制研究所 宇航與智能控制技術(shù)國家級重點實驗室, 北京 100854)

    研究了不同飛行參數(shù)下球錐形光學頭罩的高超聲速繞流流場的數(shù)值模擬問題。首先,分別利用k-ε湍流模型和SA湍流模型求解N-S方程,獲得了球錐形光學頭罩繞流流場的模擬結(jié)果。然后,將模擬結(jié)果與風洞測試得到的實驗數(shù)據(jù)進行比較分析。結(jié)果表明,在高馬赫數(shù)、大迎角情況下采用k-ε湍流模型可以得到相對更準確的計算結(jié)果。最后,對流場的溫度、壓強、密度分布特性進行了分析,給出了球錐形頭罩繞流流場的壓強、溫度和密度隨飛行參數(shù)的變化規(guī)律。

    球錐形頭罩; 流場; 數(shù)值模擬; 湍流模型

    0 引言

    隨著現(xiàn)代化戰(zhàn)爭對導彈飛行速度和突防能力要求的不斷提高,高超聲速導彈已成為當前的研究熱點。導彈頭罩與周圍氣體劇烈作用會形成復雜的流場,進而產(chǎn)生氣動光學效應(yīng)[1],影響光學制導系統(tǒng)對目標的探測、識別和跟蹤的能力。球錐形光學頭罩具有較好的空氣動力學特性,可以減小導彈在高速飛行時受到的空氣阻力和摩擦熱,現(xiàn)已廣泛應(yīng)用于空空導彈(如“霹靂”7、“楊樹”R-27T、“怪蛇”3等)。進行球錐形頭罩的繞流流場數(shù)值模擬,不僅對研究導彈的氣動特性具有重要意義,而且還是實現(xiàn)氣動光學效應(yīng)校正、提高導彈制導精度的關(guān)鍵技術(shù)之一。

    對流場參數(shù)計算的準確性在很大程度上取決于所選擇的湍流模型。目前對高速流場的數(shù)值模擬主要采用的是k-ε湍流模型和SA(Spalart-Allmaras)湍流模型,并且研究對象主要集中于飛行器整體、翼型、平板形側(cè)窗等[2-7]。如劉景源[4]、于子文[5]利用k-ε湍流模型分別對平板繞流、前飛旋翼流場進行了數(shù)值模擬,周大高[6]、段卓毅[7]利用SA湍流模型分別對翼型S825繞流、某機翼副翼舵面流場進行了數(shù)值模擬。這些研究都得到了與實驗數(shù)據(jù)吻合較好的計算結(jié)果,但是這兩種湍流模型對于球錐形頭罩高超聲速繞流流場是否適用并不確定。因此,本文將首先對這兩種湍流模型在計算球錐形頭罩高超聲速繞流流場數(shù)值模擬中的有效性進行討論;然后通過對模擬結(jié)果的分析,研究不同馬赫數(shù)和不同迎角條件下繞流流場的溫度、壓強、密度等參數(shù)的分布特性,為分析球錐形頭罩的氣動特性和氣動光學效應(yīng)提供數(shù)據(jù)依據(jù)。

    1 計算方法

    1.1 控制方程

    流體運動的基本方程是建立在質(zhì)量守恒定律、牛頓運動定律和能量守恒定律的基礎(chǔ)上的。習慣上將連續(xù)方程、動量方程和能量方程組成的流體運動方程稱為N-S方程組。在三維笛卡兒直角坐標系中,可壓縮粘性流體N-S方程組的守恒積分形式為:

    (1)

    式中,W為流場守恒矢量;V為任意控制體體積;?V為控制體邊界;H為矢量通量,可分解為無粘對流矢量通量Hc和粘性矢量通量Hv之差的形式;n為?V表面的單位外法矢量;S為面積。W,Hc,Hv的具體表達式可參見文獻[1]。

    為封閉方程組,還需引入4個熱力學關(guān)系式,分別為完全氣體狀態(tài)方程:

    p=ρRT

    (2)

    單位質(zhì)量氣體的總能量:

    e=p/[(γ-1)ρ]+ujuj/2

    (3)

    由Sutherland公式計算得到的動力粘性系數(shù)為:

    μ=C1T3/2/(T+C2)

    (4)

    式中,C1=1.458×10-6kg/m·s·K1/2;C2=110.4 K。

    對于各向同性流體,其導熱系數(shù)為:

    k=μCP/Pr=μγR/[(γ-1)Pr]

    (5)

    式中,γ為氣體比熱系數(shù),對于完全氣體通常γ=1.4;R為氣體常數(shù);Pr為普朗特數(shù),對于完全氣體通常Pr=0.72。

    1.2 湍流模型

    k-ε二方程湍流模型是由Launder等人在1972年提出的,因其穩(wěn)定性、經(jīng)濟性和計算精度較高,現(xiàn)已成為湍流模型中應(yīng)用最廣泛、最為人熟知的。SA湍流模型是由Spalart和Allmaras在1992年提出的,是相對簡單的一方程模型,對流動分離區(qū)附近的計算效果較好,常用于模擬翼形、壁面邊界等流動,現(xiàn)已廣泛應(yīng)用于航空領(lǐng)域。

    在k-ε二方程湍流模型中,通過求解湍流動能k方程和湍動耗散率ε方程,得到k和ε的解,然后再用k和ε的值計算湍流粘度μt,最終通過Boussinesq假設(shè)得到雷諾應(yīng)力的解。其相應(yīng)的方程和函數(shù)表達式為[8-9]:

    (6)

    (7)

    (8)

    式中,模型系數(shù)Cμ=0.09,Cε1=1.44,Cε2=1.92,σk=1.0,σε=1.3。

    (9)

    式中,右邊四項分別為湍流粘性的生成項、耗散項、自定義源項和解體項。

    (10)

    (11)

    式(9)~式(11)中的各項具體表達式和模型系數(shù)可參見文獻[10-11]。

    1.3 網(wǎng)格劃分

    圖1 計算網(wǎng)格示意圖Fig.1 Scheme of computational grids

    1.4 邊界條件

    上游邊界:采用無窮遠處自由可壓來流作為壓力遠場邊界條件,給定來流的靜壓、靜溫和馬赫數(shù);出口邊界:采用壓力出口邊界條件,對于達到高超聲速的流動,參數(shù)都由流場內(nèi)部通過插值外推得到;物面邊界:采用無滑移壁面邊界條件。

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

    2.1 湍流模型有效性分析

    以文獻[12]中風洞測試的實驗條件為輸入?yún)?shù),設(shè)來流的初始壓強為1.5 MPa,初始溫度為370 K,導彈飛行馬赫數(shù)為5和8,迎角α為0°,5°,10°,15°,對圖1中所示光學頭罩的繞流流場進行數(shù)值模擬。計算過程中分別采用k-ε湍流模型和SA湍流模型對N-S方程進行求解,并將計算得到的頭罩表面壓力系數(shù)與風洞實驗測得的頭罩表面壓力系數(shù)進行了比較。

    圖2和圖3中給出了z=0對稱截面上頭罩表面壓力系數(shù)沿軸向距離的分布情況,其中Cp0表示風洞實驗測得的壓力系數(shù),Cp1表示采用k-ε湍流模型時計算得到的壓力系數(shù),Cp2表示采用SA湍流模型時計算得到的壓力系數(shù)。從圖中可以看出,當飛行馬赫數(shù)為5、迎角為0°~15°時,兩種模型的計算結(jié)果與實驗數(shù)據(jù)吻合程度都比較好;當飛行馬赫數(shù)為8、迎角為0°和5°時,兩種模型的吻合程度也都比較好;但當馬赫數(shù)為8、迎角為10°和15°時,k-ε湍流模型的計算結(jié)果比SA湍流模型的計算結(jié)果更為準確。這是由于SA模型運輸方程中的常系數(shù)是基于簡單流動分析得到的,在這些簡單流動中湍流運輸特性一般處于平衡狀態(tài)。但在高馬赫數(shù)、大迎角情況下湍流的非平衡運輸現(xiàn)象明顯,湍流粘性的生成和耗散之間的平衡關(guān)系發(fā)生了改變,導致計算結(jié)果出現(xiàn)了偏差??梢娫陲w行條件較為苛刻的情況下,k-ε湍流模型比SA湍流模型更為適合。

    圖2 Ma=5時壓力系數(shù)分布Fig.2 Distribution of the pressure coefficient when Ma=5

    圖3 Ma=8時壓力系數(shù)分布Fig.3 Distribution of the pressure coefficient when Ma=8

    2.2 流場參數(shù)的分布情況

    以飛行高度為5 km、馬赫數(shù)分別為5和8、迎角分別為0°和15°的飛行參數(shù)為例,利用k-ε湍流模型對光學頭罩的繞流流場進行了數(shù)值模擬。圖4~圖6給出了z=0對稱截面上頭罩周圍流場的壓強、溫度、密度的等值分布云圖。從圖中可以看出,頭罩頂部區(qū)域的壓強、溫度、密度都最大,在Ma=5時壓強高達1.6 MPa,溫度高達1 500 K,密度高達4 kg/m3;在Ma=8時則分別高達4 MPa,3 500 K,4.5 kg/m3。這是由于來流與光學頭罩相遇時受到強烈壓縮,壓強、密度大為增加;空氣本身又具有粘性,與頭罩表面相接觸時受到阻滯,使得氣流的速度大為降低,在頭罩表面附近形成邊界層,在邊界層內(nèi)來流的動能幾乎全部轉(zhuǎn)化為內(nèi)能,使光學頭罩周圍的氣流溫度迅速升高。而頭罩后部的流場充分發(fā)展為湍流,形成弓形脫體激波,分子運動耗散掉湍動能,使得壓強、溫度、密度都有所降低。

    同時從圖中還可以看出,隨著迎角的增加,駐點位置逐漸移至迎風面,頭罩弓形激波在迎風面更加貼近物面,脫體距離減少,同一位置對應(yīng)的壓強、溫度、密度都有所增大;在背風面激波則遠離物面,脫體距離增大,同一位置對應(yīng)的壓強、溫度、密度有所減小。另外,隨著馬赫數(shù)的增大,同一位置對應(yīng)的壓強、溫度、密度都有所增大,激波脫體距離有所減小。

    圖4 壓強分布云圖Fig.4 Contour map of pressure

    圖5 溫度分布云圖Fig.5 Contour map of temperature

    圖6 密度分布云圖Fig.6 Contour map of density

    3 結(jié)束語

    本文的分析表明,利用數(shù)值計算的方法模擬光學頭罩高超聲速繞流流場是有效的,且在高馬赫數(shù)、大迎角情況下采用k-ε湍流模型計算得到的結(jié)果比采用SA湍流模型計算得到的結(jié)果更為準確。通過分析不同飛行速度、不同迎角情況下流場壓強、溫度、密度等參數(shù)的分布情況可以看出,頭罩頂部駐點位置的壓強、溫度、密度最大。由此可知,此處的氣動加熱效應(yīng)和熱輻射效應(yīng)最為嚴重,對光傳輸及目標成像的影響最為明顯,需重點研究氣動光學效應(yīng)的校正方法。此外,在大迎角飛行時,迎風面的壓強、溫度、密度較大,也需考慮氣動光學效應(yīng)的影響。

    總的來說,數(shù)值模擬方法與風洞實驗相比更具有靈活性和經(jīng)濟性,針對導彈不同的飛行參數(shù),只需改變計算的初始條件、邊界條件即可得到相應(yīng)的數(shù)值解算結(jié)果,從而為定量研究球錐形頭罩的氣動特性和氣動光學效應(yīng)提供數(shù)據(jù)支持。

    [1] 殷興良.氣動光學原理[M].北京:中國宇航出版社,2003:55-62.

    [2] Tilmann C P,Buter T A,Bowersox R D W. Characterization of the flowfield near a wrap-around fin at mach number 2.8 [R].AIAA-97-0522,1997.

    [3] 陳澄,費錦東.側(cè)窗頭罩高速層流流場光學傳輸效應(yīng)數(shù)值模擬[J].紅外與激光工程,2005,34(5):548-552.

    [4] 劉景源,李椿萱.基于Reynolds平均的高超聲速二方程湍流模型[J].航空學報,2008,29(1):28-33.

    [5] 于子文,曹義華.前飛旋翼三維湍流場的數(shù)值模擬[J].北京航空航天大學學報,2006,32(7):751-755.

    [6] 周大高,柳陽威.改進SA模型對翼型分離流動的數(shù)值模擬[J].北京航空航天大學學報,2012,38(10):1384-1388.

    [7] 段卓毅,王小震,陳迎春.三維舵面繞流的N-S方程計算研究[J].飛行力學,2007,25(3):67-70.

    [8] Jones W P,Launder B E.The prediction of laminarization with a two-equation model of turbulence[J].Journal of Heat and Mass Transfer,1972,15(2):301-314.

    [9] Jones W P,Launder B E.The calculation of low-reynolds number phenomena with a two-equation model of turbulence[J].Journal of Heat and Mass Transfer,1973,16(6):1119-1130.

    [10] Spalart P R,Allmaras S R.A one-equation turbulence transport model for aerodynamic flows[R].AIAA-92-0439,1992.

    [11] Paciorri R,Dieudonne W,Degrez G,et al.Validation of the spalart-allmaras turbulence model for application in hypersonic flows[R].AIAA-97-2323,1997.

    [12] 李素循.典型外形高超聲速流動特性[M].北京:國防工業(yè)出版社,2007:80-109.

    Numericalsimulationofhypersonicflowfieldsaroundsphere-cone-shapeddome

    XU Dong1, LI Xin1, XU Li2, WANG Ya-hui2

    (1.School of Instrument Science and Opto-electronics Engineering, BUAA, Beijing 100191, China; 2.National Key Laboratory of Science and Technology on Aerospace Intelligence Control, Beijing Aerospace Automatic Control Institute, Beijing 100854, China)

    The numerical simulation of hypersonic flow fields around sphere-cone-shaped dome under different flight conditions was studied. Firstly, by separately applyingk-εturbulence model and SA (Spalart-Allmaras) turbulence model to solve the N-S equations, the flow fields around sphere-cone-shaped optical dome were simulated. Secondly, the simulation results were compared with related experimental data, and the comparison shows that, under the condition of high Mach number and attack-angle more accurate calculation results can be obtained when usingk-εturbulence model. Lastly, the distribution of pressure, temperature and density of the flow fields around sphere-cone-shaped dome was analyzed, and the variation of the three parameters with flight condition was also presented.

    sphere-cone-shaped dome; flow fields; numerical simulation; turbulence model

    V211.3

    A

    1002-0853(2013)06-0491-05

    2013-03-29;

    2013-09-01; < class="emphasis_bold">網(wǎng)絡(luò)出版時間

    時間:2013-10-22 14:15

    國家自然科學基金資助(61378077)

    許東(1973-),男,江蘇徐州人,副教授,博士,主要研究方向為氣動光學效應(yīng)及目標探測識別;

    李欣(1988-),女,山東濟南人,碩士研究生,主要研究方向為氣動光學及流場仿真。

    (編輯:姚妙慧)

    猜你喜歡
    頭罩迎角湍流
    連續(xù)變迎角試驗數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    基于非結(jié)構(gòu)嵌套網(wǎng)格的低空大動壓頭罩分離數(shù)值模擬
    重氣瞬時泄漏擴散的湍流模型驗證
    USB接口式防花粉頭罩
    小康(2014年4期)2015-01-29 01:12:28
    超聲速氟化鎂紅外頭罩熱強度失效判別準則研究
    上海航天(2014年3期)2014-12-31 11:48:24
    失速保護系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    弱分層湍流輸運特性的統(tǒng)計分析
    新生兒低流量頭罩吸氧的安全性研究
    久久精品夜色国产| 偷拍熟女少妇极品色| 最近最新中文字幕免费大全7| 午夜福利在线在线| 国产高潮美女av| 99re6热这里在线精品视频| 国产精品国产av在线观看| 亚洲丝袜综合中文字幕| 亚洲在线观看片| 日韩亚洲欧美综合| 黄色怎么调成土黄色| 日韩av免费高清视频| 啦啦啦在线观看免费高清www| 亚洲av免费在线观看| 亚洲伊人久久精品综合| 最新中文字幕久久久久| 18禁裸乳无遮挡动漫免费视频 | 日韩一区二区三区影片| 伊人久久国产一区二区| 国产午夜精品一二区理论片| 欧美激情久久久久久爽电影| 免费大片黄手机在线观看| 特级一级黄色大片| 麻豆成人av视频| 国产精品99久久久久久久久| 人人妻人人爽人人添夜夜欢视频 | 亚洲成色77777| 国产精品99久久久久久久久| 天天一区二区日本电影三级| 97超碰精品成人国产| 欧美老熟妇乱子伦牲交| 欧美老熟妇乱子伦牲交| 91精品一卡2卡3卡4卡| 国产精品国产av在线观看| 99久久精品一区二区三区| 精品视频人人做人人爽| 看黄色毛片网站| 国产成人一区二区在线| 久久亚洲国产成人精品v| 人人妻人人爽人人添夜夜欢视频 | 日韩一本色道免费dvd| 欧美精品一区二区大全| 国产女主播在线喷水免费视频网站| 免费观看av网站的网址| 亚洲成人av在线免费| 成人免费观看视频高清| 少妇高潮的动态图| 亚洲av二区三区四区| 在线a可以看的网站| 插逼视频在线观看| av网站免费在线观看视频| 久久精品人妻少妇| 久久久久久九九精品二区国产| 三级男女做爰猛烈吃奶摸视频| 三级男女做爰猛烈吃奶摸视频| 成人亚洲精品一区在线观看 | 亚洲av不卡在线观看| 建设人人有责人人尽责人人享有的 | 久久久久久久久久久免费av| 国产毛片a区久久久久| 日韩国内少妇激情av| 看免费成人av毛片| 精品酒店卫生间| 国产欧美亚洲国产| 不卡视频在线观看欧美| 午夜视频国产福利| 国产成人一区二区在线| 亚洲在久久综合| 久久女婷五月综合色啪小说 | 国产精品av视频在线免费观看| 亚洲国产精品专区欧美| 国产精品人妻久久久影院| 亚洲精品色激情综合| 有码 亚洲区| 在线观看美女被高潮喷水网站| 91久久精品国产一区二区成人| 好男人视频免费观看在线| av卡一久久| 中文乱码字字幕精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 久久99热6这里只有精品| 亚洲在久久综合| 免费av观看视频| 日韩,欧美,国产一区二区三区| 夫妻性生交免费视频一级片| 国产精品久久久久久久电影| 成人毛片60女人毛片免费| 交换朋友夫妻互换小说| 亚洲真实伦在线观看| 亚洲熟女精品中文字幕| 美女国产视频在线观看| 国产男人的电影天堂91| 一级a做视频免费观看| 深爱激情五月婷婷| 色播亚洲综合网| 国模一区二区三区四区视频| 熟女人妻精品中文字幕| 如何舔出高潮| 亚洲精品乱码久久久v下载方式| 国产大屁股一区二区在线视频| 搞女人的毛片| 成年女人看的毛片在线观看| 日韩一本色道免费dvd| av国产免费在线观看| 国产成人免费观看mmmm| 久久久久久国产a免费观看| 欧美性感艳星| 午夜亚洲福利在线播放| 禁无遮挡网站| 性色av一级| 别揉我奶头 嗯啊视频| 欧美潮喷喷水| 尾随美女入室| 久久久久久久久久久免费av| 极品少妇高潮喷水抽搐| 久久久亚洲精品成人影院| 中国三级夫妇交换| 熟女人妻精品中文字幕| 亚洲欧美精品自产自拍| 纵有疾风起免费观看全集完整版| 久久99蜜桃精品久久| 精品一区二区免费观看| 国产亚洲91精品色在线| 视频中文字幕在线观看| 午夜福利视频1000在线观看| 51国产日韩欧美| 亚洲人成网站在线播| 夜夜爽夜夜爽视频| 成人高潮视频无遮挡免费网站| 成人毛片a级毛片在线播放| 国产精品秋霞免费鲁丝片| 国产成人福利小说| 国产成人a∨麻豆精品| 亚洲欧美中文字幕日韩二区| 一级av片app| 在线观看人妻少妇| 亚洲人成网站在线观看播放| 人人妻人人澡人人爽人人夜夜| 成人毛片60女人毛片免费| 精品国产乱码久久久久久小说| 2018国产大陆天天弄谢| 黄色怎么调成土黄色| 简卡轻食公司| 熟妇人妻不卡中文字幕| 日韩电影二区| 内射极品少妇av片p| 肉色欧美久久久久久久蜜桃 | 国产美女午夜福利| 欧美高清性xxxxhd video| 赤兔流量卡办理| 欧美最新免费一区二区三区| 日韩av在线免费看完整版不卡| 国产成人91sexporn| 有码 亚洲区| 久久韩国三级中文字幕| 男人和女人高潮做爰伦理| 美女高潮的动态| a级一级毛片免费在线观看| 99久久九九国产精品国产免费| 午夜日本视频在线| 亚洲国产av新网站| 午夜免费鲁丝| 亚洲国产精品成人久久小说| 成人亚洲精品av一区二区| 亚洲精华国产精华液的使用体验| 欧美日韩综合久久久久久| 欧美亚洲 丝袜 人妻 在线| 18禁裸乳无遮挡动漫免费视频 | 久久久久久久亚洲中文字幕| 国产毛片在线视频| 亚洲欧美一区二区三区黑人 | 在线观看美女被高潮喷水网站| 亚洲最大成人手机在线| 欧美 日韩 精品 国产| 国产探花在线观看一区二区| 国产精品一区二区三区四区免费观看| 日本-黄色视频高清免费观看| 国产真实伦视频高清在线观看| 美女国产视频在线观看| 久久人人爽人人爽人人片va| 日韩成人av中文字幕在线观看| 亚洲欧美清纯卡通| 成人鲁丝片一二三区免费| 亚洲天堂国产精品一区在线| 又黄又爽又刺激的免费视频.| 91精品一卡2卡3卡4卡| 老司机影院毛片| 交换朋友夫妻互换小说| 26uuu在线亚洲综合色| av在线天堂中文字幕| 亚洲综合精品二区| 成人二区视频| 看黄色毛片网站| 在线天堂最新版资源| 免费观看a级毛片全部| 大片电影免费在线观看免费| 黄色怎么调成土黄色| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av天美| 日韩一区二区视频免费看| 九草在线视频观看| 我的老师免费观看完整版| 哪个播放器可以免费观看大片| 亚洲国产欧美人成| 视频中文字幕在线观看| 成年人午夜在线观看视频| 亚洲av在线观看美女高潮| 日韩成人伦理影院| 亚洲精品aⅴ在线观看| 国产 一区精品| 国产黄a三级三级三级人| 18禁在线无遮挡免费观看视频| 晚上一个人看的免费电影| 18禁动态无遮挡网站| 欧美老熟妇乱子伦牲交| 夜夜看夜夜爽夜夜摸| 欧美日韩精品成人综合77777| 超碰av人人做人人爽久久| 人妻少妇偷人精品九色| 乱码一卡2卡4卡精品| 91狼人影院| 久久久久久久久久人人人人人人| 日韩一本色道免费dvd| 日日撸夜夜添| 天美传媒精品一区二区| 久久精品夜色国产| 国产色爽女视频免费观看| 成人高潮视频无遮挡免费网站| 亚洲精品视频女| 嫩草影院精品99| 亚洲精品成人av观看孕妇| 天堂网av新在线| 在线a可以看的网站| 久久久久精品久久久久真实原创| 爱豆传媒免费全集在线观看| 嫩草影院入口| 在线 av 中文字幕| 老师上课跳d突然被开到最大视频| 天天一区二区日本电影三级| 久久国产乱子免费精品| 男女边吃奶边做爰视频| 亚洲欧美中文字幕日韩二区| 亚洲美女搞黄在线观看| 久久久久国产网址| 成人黄色视频免费在线看| 99九九线精品视频在线观看视频| 亚洲欧美中文字幕日韩二区| 国产91av在线免费观看| 欧美日本视频| 久久久午夜欧美精品| 国产精品爽爽va在线观看网站| 亚洲婷婷狠狠爱综合网| 又爽又黄无遮挡网站| 国产伦精品一区二区三区四那| 精品熟女少妇av免费看| 菩萨蛮人人尽说江南好唐韦庄| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品古装| 麻豆成人av视频| 国产精品成人在线| 黄色怎么调成土黄色| 亚洲国产成人一精品久久久| 不卡视频在线观看欧美| 日本午夜av视频| 国产免费视频播放在线视频| 国产v大片淫在线免费观看| 国产成人freesex在线| 久久精品国产自在天天线| 欧美精品人与动牲交sv欧美| 18禁裸乳无遮挡动漫免费视频 | 在线播放无遮挡| 男女边摸边吃奶| 久久影院123| 又爽又黄无遮挡网站| 国产精品99久久久久久久久| 亚洲自拍偷在线| 久久久久久久精品精品| 欧美3d第一页| 欧美+日韩+精品| 亚洲精品日韩在线中文字幕| 精品一区在线观看国产| 国产乱来视频区| 色综合色国产| xxx大片免费视频| 一本一本综合久久| 男女那种视频在线观看| 亚洲av日韩在线播放| 精品久久久久久久久av| 国产精品久久久久久精品电影小说 | 亚洲av一区综合| 亚洲精品成人av观看孕妇| 人妻一区二区av| 看非洲黑人一级黄片| 日日摸夜夜添夜夜爱| 美女视频免费永久观看网站| 五月开心婷婷网| 亚洲欧洲国产日韩| 一级毛片我不卡| 精品久久久噜噜| 国产免费福利视频在线观看| 麻豆国产97在线/欧美| 秋霞在线观看毛片| 人妻一区二区av| 一区二区av电影网| 18禁裸乳无遮挡免费网站照片| 精品一区二区三卡| 亚洲精品日韩在线中文字幕| 观看美女的网站| 亚洲久久久久久中文字幕| 国产精品久久久久久久电影| 建设人人有责人人尽责人人享有的 | 久久精品熟女亚洲av麻豆精品| 国产午夜福利久久久久久| 97超视频在线观看视频| 国语对白做爰xxxⅹ性视频网站| 高清毛片免费看| 亚洲成人av在线免费| 日韩三级伦理在线观看| 99久国产av精品国产电影| 亚洲av成人精品一区久久| av在线蜜桃| 国国产精品蜜臀av免费| 三级国产精品片| 国产毛片a区久久久久| 日韩不卡一区二区三区视频在线| 亚洲va在线va天堂va国产| 又黄又爽又刺激的免费视频.| 天天躁夜夜躁狠狠久久av| 日韩成人av中文字幕在线观看| 特大巨黑吊av在线直播| 国产高清不卡午夜福利| 各种免费的搞黄视频| 亚洲精品日本国产第一区| 日韩在线高清观看一区二区三区| 欧美日韩在线观看h| 一二三四中文在线观看免费高清| 国产综合懂色| 三级国产精品欧美在线观看| 亚洲最大成人av| 香蕉精品网在线| 美女脱内裤让男人舔精品视频| 我要看日韩黄色一级片| 亚洲精品456在线播放app| 国产毛片在线视频| 亚洲精品成人久久久久久| 久久97久久精品| 中文字幕亚洲精品专区| 大陆偷拍与自拍| 777米奇影视久久| 亚洲人与动物交配视频| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 久久久色成人| 成年女人在线观看亚洲视频 | 最近2019中文字幕mv第一页| 搞女人的毛片| 免费电影在线观看免费观看| 人妻少妇偷人精品九色| 久久久久久久久久成人| 七月丁香在线播放| 国产精品久久久久久精品电影小说 | 久久久久网色| 国产又色又爽无遮挡免| 亚洲精品日韩av片在线观看| 欧美成人一区二区免费高清观看| 大又大粗又爽又黄少妇毛片口| 午夜福利在线在线| 国产成人免费无遮挡视频| 亚洲在线观看片| 亚洲av日韩在线播放| 97在线视频观看| 午夜精品国产一区二区电影 | 3wmmmm亚洲av在线观看| 一区二区三区乱码不卡18| 少妇的逼水好多| eeuss影院久久| 免费观看的影片在线观看| 亚洲av中文字字幕乱码综合| av播播在线观看一区| 97热精品久久久久久| 亚洲丝袜综合中文字幕| av国产免费在线观看| 欧美zozozo另类| 欧美97在线视频| 亚洲自拍偷在线| 国产免费视频播放在线视频| 国产免费又黄又爽又色| 在线观看av片永久免费下载| 精品少妇久久久久久888优播| av黄色大香蕉| 男女啪啪激烈高潮av片| 嫩草影院新地址| 国产综合懂色| 自拍偷自拍亚洲精品老妇| 久久久久九九精品影院| 看黄色毛片网站| 亚洲怡红院男人天堂| 国产精品熟女久久久久浪| 午夜视频国产福利| 看非洲黑人一级黄片| 五月天丁香电影| 国产探花极品一区二区| 国产精品人妻久久久影院| 啦啦啦啦在线视频资源| 肉色欧美久久久久久久蜜桃 | 久久精品国产鲁丝片午夜精品| 国产女主播在线喷水免费视频网站| 色播亚洲综合网| 免费人成在线观看视频色| 日韩中字成人| 国产人妻一区二区三区在| 特级一级黄色大片| 高清av免费在线| 成人鲁丝片一二三区免费| 插逼视频在线观看| 精品少妇久久久久久888优播| 欧美一级a爱片免费观看看| 欧美日韩视频精品一区| 日本欧美国产在线视频| 麻豆成人av视频| 在线观看三级黄色| 晚上一个人看的免费电影| 日韩av不卡免费在线播放| 精品久久久精品久久久| 国产毛片在线视频| 一级片'在线观看视频| 欧美极品一区二区三区四区| 交换朋友夫妻互换小说| 婷婷色综合www| 黄色配什么色好看| 精品久久久久久久久av| 又黄又爽又刺激的免费视频.| 久久99热这里只有精品18| 国产黄色视频一区二区在线观看| 日本一本二区三区精品| 中文字幕av成人在线电影| 久久午夜福利片| 在线播放无遮挡| 在线观看三级黄色| 内地一区二区视频在线| 欧美精品人与动牲交sv欧美| 男男h啪啪无遮挡| 麻豆成人av视频| 大片电影免费在线观看免费| 香蕉精品网在线| 亚洲激情五月婷婷啪啪| 婷婷色综合www| 看十八女毛片水多多多| 久久影院123| 国产免费又黄又爽又色| 久久精品国产自在天天线| 国产精品99久久99久久久不卡 | 超碰97精品在线观看| 两个人的视频大全免费| 97热精品久久久久久| 精品少妇黑人巨大在线播放| 精品午夜福利在线看| 九九爱精品视频在线观看| 18在线观看网站| 免费在线观看视频国产中文字幕亚洲 | 少妇猛男粗大的猛烈进出视频| 丝袜脚勾引网站| 人人妻人人澡人人爽人人夜夜| 精品一区在线观看国产| 天天影视国产精品| 亚洲欧洲国产日韩| 亚洲精品在线美女| 日韩一本色道免费dvd| 夜夜骑夜夜射夜夜干| 不卡视频在线观看欧美| √禁漫天堂资源中文www| 极品少妇高潮喷水抽搐| 黄片无遮挡物在线观看| 亚洲国产精品成人久久小说| 夫妻性生交免费视频一级片| 久久综合国产亚洲精品| 久久久精品免费免费高清| 七月丁香在线播放| 999久久久国产精品视频| 青草久久国产| 水蜜桃什么品种好| 免费观看人在逋| 啦啦啦在线免费观看视频4| 黄片小视频在线播放| 国产精品亚洲av一区麻豆 | 亚洲精品,欧美精品| 中文乱码字字幕精品一区二区三区| 欧美最新免费一区二区三区| 国产成人av激情在线播放| 考比视频在线观看| 大码成人一级视频| 亚洲精品一区蜜桃| 日本午夜av视频| 热re99久久精品国产66热6| 亚洲欧美一区二区三区黑人| 久久国产亚洲av麻豆专区| 捣出白浆h1v1| 9色porny在线观看| 国产av一区二区精品久久| www.精华液| 国产日韩欧美在线精品| 亚洲国产av影院在线观看| 婷婷色av中文字幕| 日韩一本色道免费dvd| 下体分泌物呈黄色| 男女之事视频高清在线观看 | 人人妻人人澡人人看| 老司机靠b影院| 久久久亚洲精品成人影院| 热99国产精品久久久久久7| 最近中文字幕2019免费版| 精品少妇久久久久久888优播| 国产av一区二区精品久久| 成年女人毛片免费观看观看9 | 天天躁夜夜躁狠狠久久av| 亚洲欧美清纯卡通| 中文字幕另类日韩欧美亚洲嫩草| 毛片一级片免费看久久久久| 性高湖久久久久久久久免费观看| 91精品国产国语对白视频| 亚洲国产日韩一区二区| 制服人妻中文乱码| 纯流量卡能插随身wifi吗| 一区在线观看完整版| 黄网站色视频无遮挡免费观看| 汤姆久久久久久久影院中文字幕| 男人爽女人下面视频在线观看| 人体艺术视频欧美日本| 99九九在线精品视频| 久久国产精品大桥未久av| 在线观看免费高清a一片| 成人免费观看视频高清| 卡戴珊不雅视频在线播放| 亚洲在久久综合| 看免费av毛片| 国产探花极品一区二区| 综合色丁香网| 99国产精品免费福利视频| 青春草亚洲视频在线观看| 久热这里只有精品99| 日本猛色少妇xxxxx猛交久久| 免费少妇av软件| xxx大片免费视频| 美女中出高潮动态图| 大片电影免费在线观看免费| 日韩视频在线欧美| 又粗又硬又长又爽又黄的视频| 国产又爽黄色视频| 国产一区亚洲一区在线观看| 日本91视频免费播放| 各种免费的搞黄视频| 一个人免费看片子| 午夜激情久久久久久久| 久热这里只有精品99| 国产精品 欧美亚洲| 久久精品aⅴ一区二区三区四区| √禁漫天堂资源中文www| h视频一区二区三区| 男女高潮啪啪啪动态图| 宅男免费午夜| 亚洲欧洲精品一区二区精品久久久 | 国产在线一区二区三区精| 99久久精品国产亚洲精品| 女人高潮潮喷娇喘18禁视频| 一级毛片 在线播放| 女人高潮潮喷娇喘18禁视频| 精品久久久久久电影网| 男女午夜视频在线观看| av片东京热男人的天堂| 搡老岳熟女国产| 国产精品免费大片| 国产精品国产av在线观看| 国产高清国产精品国产三级| 熟女av电影| 久久青草综合色| 五月开心婷婷网| 丝袜美足系列| 男女之事视频高清在线观看 | 成人影院久久| 欧美老熟妇乱子伦牲交| 黄片播放在线免费| www.av在线官网国产| 久久韩国三级中文字幕| 一级毛片 在线播放| 女人高潮潮喷娇喘18禁视频| 久久久久视频综合| 久久狼人影院| 国产片内射在线| 亚洲国产精品999| 亚洲色图综合在线观看| 成人亚洲精品一区在线观看| 日韩欧美精品免费久久| 男女之事视频高清在线观看 | 成人三级做爰电影| 日韩精品免费视频一区二区三区| www.熟女人妻精品国产| 一本一本久久a久久精品综合妖精| 精品少妇久久久久久888优播| 亚洲精品美女久久久久99蜜臀 | 毛片一级片免费看久久久久| 免费高清在线观看日韩| 久久精品亚洲熟妇少妇任你| 亚洲av欧美aⅴ国产| 日韩一区二区三区影片| av免费观看日本| 嫩草影视91久久| 七月丁香在线播放| 成年动漫av网址| 中文字幕最新亚洲高清| 亚洲精品美女久久久久99蜜臀 | 国产1区2区3区精品| av有码第一页| 啦啦啦视频在线资源免费观看| 亚洲久久久国产精品| 国产精品欧美亚洲77777|