• <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)計分析
    新生兒低流量頭罩吸氧的安全性研究
    熟女电影av网| 人妻丰满熟妇av一区二区三区| 午夜福利免费观看在线| 亚洲av美国av| 97人妻精品一区二区三区麻豆| 好男人电影高清在线观看| 久久久久久人人人人人| 婷婷精品国产亚洲av| 亚洲国产精品合色在线| 不卡av一区二区三区| 午夜日韩欧美国产| 看片在线看免费视频| 欧美一级a爱片免费观看看 | 搡老熟女国产l中国老女人| 两个人看的免费小视频| 正在播放国产对白刺激| 成人国语在线视频| 夜夜夜夜夜久久久久| e午夜精品久久久久久久| 99久久精品国产亚洲精品| 久久久久性生活片| 久久久精品国产亚洲av高清涩受| 国产精品美女特级片免费视频播放器 | 91麻豆av在线| 国产欧美日韩一区二区精品| 一区二区三区国产精品乱码| 怎么达到女性高潮| 亚洲第一电影网av| 成人国语在线视频| 狂野欧美激情性xxxx| 成人亚洲精品av一区二区| 国产成人啪精品午夜网站| 久久久久精品国产欧美久久久| 国产av一区在线观看免费| 正在播放国产对白刺激| 久久国产精品人妻蜜桃| 欧美一级毛片孕妇| 欧美极品一区二区三区四区| 蜜桃久久精品国产亚洲av| 中出人妻视频一区二区| 99热这里只有精品一区 | 亚洲熟妇熟女久久| 长腿黑丝高跟| 国产三级黄色录像| av有码第一页| 久久中文字幕人妻熟女| 欧美成人性av电影在线观看| 一个人观看的视频www高清免费观看 | 久久久久免费精品人妻一区二区| 亚洲一区二区三区不卡视频| 免费在线观看视频国产中文字幕亚洲| 国内精品一区二区在线观看| 日本熟妇午夜| 在线国产一区二区在线| 99久久99久久久精品蜜桃| 在线观看舔阴道视频| 9191精品国产免费久久| 在线国产一区二区在线| 久久久久免费精品人妻一区二区| 又黄又粗又硬又大视频| 亚洲午夜精品一区,二区,三区| 欧美日韩黄片免| 国产三级中文精品| 久久亚洲精品不卡| 91av网站免费观看| 好男人在线观看高清免费视频| 国产麻豆成人av免费视频| 热99re8久久精品国产| 国产成人一区二区三区免费视频网站| 免费在线观看亚洲国产| 女人被狂操c到高潮| 国内久久婷婷六月综合欲色啪| svipshipincom国产片| 国产av一区在线观看免费| 精品国产乱子伦一区二区三区| 波多野结衣巨乳人妻| 黑人操中国人逼视频| 亚洲欧美日韩东京热| 国产精品98久久久久久宅男小说| 夜夜看夜夜爽夜夜摸| 每晚都被弄得嗷嗷叫到高潮| 人人妻,人人澡人人爽秒播| av天堂在线播放| 亚洲自拍偷在线| 国产一区二区在线观看日韩 | 免费看a级黄色片| 最好的美女福利视频网| 此物有八面人人有两片| 成年版毛片免费区| 久久久久久久久中文| 色尼玛亚洲综合影院| 色尼玛亚洲综合影院| 日本一本二区三区精品| 亚洲第一电影网av| 亚洲在线自拍视频| 国产亚洲精品综合一区在线观看 | 久久香蕉国产精品| 亚洲av电影在线进入| 极品教师在线免费播放| 18禁美女被吸乳视频| av超薄肉色丝袜交足视频| 一本久久中文字幕| 欧美乱妇无乱码| 国产成人aa在线观看| 亚洲第一欧美日韩一区二区三区| 老司机福利观看| 午夜久久久久精精品| 国产探花在线观看一区二区| 欧美日韩国产亚洲二区| 成年版毛片免费区| 午夜视频精品福利| 欧美日本亚洲视频在线播放| 欧美av亚洲av综合av国产av| 亚洲片人在线观看| 在线永久观看黄色视频| 久久精品综合一区二区三区| 性色av乱码一区二区三区2| 亚洲国产欧美一区二区综合| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文字幕一区二区三区有码在线看 | 看片在线看免费视频| 波多野结衣巨乳人妻| 国产亚洲精品一区二区www| 女同久久另类99精品国产91| www.自偷自拍.com| 99re在线观看精品视频| 亚洲,欧美精品.| 国产精品日韩av在线免费观看| 欧美日韩亚洲综合一区二区三区_| 母亲3免费完整高清在线观看| 高清毛片免费观看视频网站| 欧美色欧美亚洲另类二区| 成人国语在线视频| 亚洲国产精品久久男人天堂| 久久热在线av| netflix在线观看网站| 国语自产精品视频在线第100页| 亚洲中文字幕一区二区三区有码在线看 | 国产探花在线观看一区二区| 色综合婷婷激情| 久久国产精品影院| 人妻夜夜爽99麻豆av| 亚洲成人精品中文字幕电影| 两个人免费观看高清视频| 亚洲熟妇熟女久久| 我要搜黄色片| 免费看十八禁软件| 一本一本综合久久| 99热6这里只有精品| www.熟女人妻精品国产| 两性午夜刺激爽爽歪歪视频在线观看 | 中文字幕高清在线视频| 国内精品一区二区在线观看| 亚洲成人久久性| 久久精品国产亚洲av香蕉五月| 亚洲中文av在线| 大型av网站在线播放| 99久久国产精品久久久| 国产人伦9x9x在线观看| 久久精品国产综合久久久| 色老头精品视频在线观看| 三级国产精品欧美在线观看 | 欧美av亚洲av综合av国产av| 中文字幕人妻丝袜一区二区| 叶爱在线成人免费视频播放| 精品久久久久久久末码| 午夜激情福利司机影院| 国产av一区在线观看免费| 日韩国内少妇激情av| 每晚都被弄得嗷嗷叫到高潮| 久久婷婷人人爽人人干人人爱| 舔av片在线| 久9热在线精品视频| √禁漫天堂资源中文www| 成人永久免费在线观看视频| 99国产极品粉嫩在线观看| 男女床上黄色一级片免费看| 亚洲av熟女| 日本撒尿小便嘘嘘汇集6| 中文在线观看免费www的网站 | 欧美性猛交╳xxx乱大交人| 国产av又大| 国产高清videossex| 成人永久免费在线观看视频| 一进一出抽搐gif免费好疼| 草草在线视频免费看| 叶爱在线成人免费视频播放| 欧美色视频一区免费| 伦理电影免费视频| 熟妇人妻久久中文字幕3abv| 成人av一区二区三区在线看| 91麻豆精品激情在线观看国产| 久久精品国产亚洲av高清一级| 一进一出抽搐动态| 免费观看精品视频网站| 全区人妻精品视频| 国产一区在线观看成人免费| 色综合站精品国产| 禁无遮挡网站| svipshipincom国产片| 女人爽到高潮嗷嗷叫在线视频| 一进一出抽搐gif免费好疼| 久久久久久久久久黄片| 欧美在线一区亚洲| 亚洲一区高清亚洲精品| 免费在线观看日本一区| 成人特级黄色片久久久久久久| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久毛片微露脸| 久久天躁狠狠躁夜夜2o2o| 国产v大片淫在线免费观看| 色综合婷婷激情| 老熟妇仑乱视频hdxx| 亚洲av成人av| 日本一二三区视频观看| 亚洲男人的天堂狠狠| 两性午夜刺激爽爽歪歪视频在线观看 | 日本 av在线| 久久精品国产亚洲av高清一级| 国产视频内射| 丁香六月欧美| 日本黄大片高清| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av高清一级| 性欧美人与动物交配| 在线国产一区二区在线| 成人av一区二区三区在线看| 久久精品影院6| 最近在线观看免费完整版| 最近最新中文字幕大全电影3| 欧美日韩瑟瑟在线播放| 久久久久亚洲av毛片大全| 成人国产一区最新在线观看| 日韩欧美三级三区| 午夜福利视频1000在线观看| 制服人妻中文乱码| 一二三四在线观看免费中文在| 男人舔女人下体高潮全视频| 国产精品av久久久久免费| 欧美人与性动交α欧美精品济南到| 国产高清videossex| 女同久久另类99精品国产91| 精品人妻1区二区| 亚洲狠狠婷婷综合久久图片| 青草久久国产| 一本精品99久久精品77| 成年版毛片免费区| 亚洲成人中文字幕在线播放| 99国产综合亚洲精品| tocl精华| 国产精品一区二区免费欧美| 女同久久另类99精品国产91| 国产一区在线观看成人免费| 欧美成人性av电影在线观看| 夜夜夜夜夜久久久久| 成人精品一区二区免费| 欧美日本视频| 麻豆久久精品国产亚洲av| 国产精品,欧美在线| 在线观看午夜福利视频| 国产av一区在线观看免费| 国产精品一区二区精品视频观看| 国产av在哪里看| 国产亚洲欧美在线一区二区| 成年人黄色毛片网站| 法律面前人人平等表现在哪些方面| 日韩精品中文字幕看吧| 日韩精品免费视频一区二区三区| 岛国在线免费视频观看| 国产激情久久老熟女| 午夜日韩欧美国产| 国产精品久久久久久亚洲av鲁大| 亚洲av电影在线进入| 97超级碰碰碰精品色视频在线观看| 真人做人爱边吃奶动态| 国产精品久久久久久亚洲av鲁大| 俄罗斯特黄特色一大片| 国产区一区二久久| 欧美精品啪啪一区二区三区| 欧美黑人精品巨大| a级毛片在线看网站| 欧洲精品卡2卡3卡4卡5卡区| 成人国产综合亚洲| 他把我摸到了高潮在线观看| 午夜福利视频1000在线观看| 高潮久久久久久久久久久不卡| 国产高清视频在线播放一区| 国产av不卡久久| 久久久国产成人免费| 成人欧美大片| 亚洲真实伦在线观看| 国产亚洲精品久久久久5区| 少妇的丰满在线观看| 99国产精品一区二区三区| 国内揄拍国产精品人妻在线| 最近最新免费中文字幕在线| 欧美久久黑人一区二区| 成人特级黄色片久久久久久久| 国产91精品成人一区二区三区| 18禁裸乳无遮挡免费网站照片| 亚洲在线自拍视频| 波多野结衣高清无吗| 国产黄色小视频在线观看| 校园春色视频在线观看| 国产亚洲av高清不卡| 国产一区二区激情短视频| 非洲黑人性xxxx精品又粗又长| 好看av亚洲va欧美ⅴa在| 五月伊人婷婷丁香| 男女视频在线观看网站免费 | 国产精品av视频在线免费观看| 69av精品久久久久久| 最近最新中文字幕大全电影3| 午夜视频精品福利| 听说在线观看完整版免费高清| 久久久久国产一级毛片高清牌| 婷婷精品国产亚洲av| 日韩欧美三级三区| 三级毛片av免费| 特级一级黄色大片| 男女下面进入的视频免费午夜| 久久性视频一级片| 成人一区二区视频在线观看| 成年版毛片免费区| 在线视频色国产色| 免费人成视频x8x8入口观看| 亚洲第一电影网av| 精品国产亚洲在线| 欧美又色又爽又黄视频| 淫秽高清视频在线观看| 1024手机看黄色片| 国产久久久一区二区三区| 久久久久久免费高清国产稀缺| 天堂√8在线中文| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av| 国产高清视频在线观看网站| 老司机午夜十八禁免费视频| 精品福利观看| 长腿黑丝高跟| 亚洲成a人片在线一区二区| 少妇熟女aⅴ在线视频| 久久精品国产清高在天天线| 好男人在线观看高清免费视频| 久久久久久久久中文| 搡老岳熟女国产| 人成视频在线观看免费观看| 一本精品99久久精品77| 国产午夜福利久久久久久| 久久久久久久午夜电影| 精品久久蜜臀av无| 精品不卡国产一区二区三区| 亚洲18禁久久av| 午夜激情福利司机影院| 国产精品一区二区三区四区免费观看 | 午夜福利免费观看在线| 很黄的视频免费| 国产精品一区二区三区四区免费观看 | 国内精品久久久久精免费| 日韩中文字幕欧美一区二区| 美女免费视频网站| 久久久久九九精品影院| 久久午夜综合久久蜜桃| 国产三级中文精品| 亚洲无线在线观看| 亚洲五月婷婷丁香| e午夜精品久久久久久久| 免费观看精品视频网站| 国产亚洲精品久久久久久毛片| 国产三级在线视频| 天堂动漫精品| cao死你这个sao货| 热99re8久久精品国产| 日本一本二区三区精品| 亚洲一区中文字幕在线| 日韩精品青青久久久久久| 午夜福利欧美成人| 中文字幕熟女人妻在线| 最近在线观看免费完整版| www.熟女人妻精品国产| 国产亚洲精品av在线| 亚洲成人久久性| 日韩欧美一区二区三区在线观看| 精品国产乱码久久久久久男人| 操出白浆在线播放| 国内精品一区二区在线观看| 亚洲精品美女久久久久99蜜臀| 久久久久性生活片| 在线a可以看的网站| 精品久久久久久久毛片微露脸| 日韩欧美在线二视频| 亚洲全国av大片| 国产精品 国内视频| 免费一级毛片在线播放高清视频| 国产一区二区在线观看日韩 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲中文字幕日韩| 91成年电影在线观看| 国产片内射在线| 人人妻,人人澡人人爽秒播| 亚洲精华国产精华精| 免费一级毛片在线播放高清视频| 国产一区在线观看成人免费| 两性夫妻黄色片| 成人高潮视频无遮挡免费网站| 深夜精品福利| 国产亚洲欧美98| 欧美三级亚洲精品| 精品免费久久久久久久清纯| 国产成人av激情在线播放| 色在线成人网| 俺也久久电影网| 国产亚洲精品一区二区www| 精品乱码久久久久久99久播| 精品一区二区三区四区五区乱码| 白带黄色成豆腐渣| 久久久精品欧美日韩精品| 午夜久久久久精精品| 国产亚洲精品第一综合不卡| 1024视频免费在线观看| 99久久99久久久精品蜜桃| 国产精品一及| 国产精品亚洲av一区麻豆| 12—13女人毛片做爰片一| 啪啪无遮挡十八禁网站| 一级毛片高清免费大全| 一a级毛片在线观看| 一进一出抽搐动态| 精品少妇一区二区三区视频日本电影| 舔av片在线| 夜夜躁狠狠躁天天躁| av国产免费在线观看| 久9热在线精品视频| 国产乱人伦免费视频| a级毛片a级免费在线| 高清在线国产一区| 国产精品久久电影中文字幕| 国内精品久久久久精免费| 草草在线视频免费看| 久久这里只有精品19| 日韩中文字幕欧美一区二区| 99国产综合亚洲精品| 最新美女视频免费是黄的| 亚洲成人精品中文字幕电影| 午夜精品在线福利| 欧美性长视频在线观看| 久久久久久久午夜电影| 欧美大码av| 男人舔女人下体高潮全视频| aaaaa片日本免费| 一边摸一边做爽爽视频免费| 欧美日韩乱码在线| 99国产精品99久久久久| 亚洲美女视频黄频| 日本一区二区免费在线视频| 麻豆成人av在线观看| 亚洲性夜色夜夜综合| 久久久久久久久免费视频了| 99re在线观看精品视频| 亚洲人成伊人成综合网2020| 精品少妇一区二区三区视频日本电影| 久久中文字幕一级| 可以在线观看毛片的网站| 国产蜜桃级精品一区二区三区| 少妇的丰满在线观看| 大型黄色视频在线免费观看| 日本一二三区视频观看| 欧美日韩瑟瑟在线播放| 麻豆国产97在线/欧美 | 欧美日韩国产亚洲二区| 亚洲精品美女久久久久99蜜臀| 久久久久国产一级毛片高清牌| 亚洲真实伦在线观看| 亚洲国产精品成人综合色| 脱女人内裤的视频| 男女下面进入的视频免费午夜| 国产日本99.免费观看| 动漫黄色视频在线观看| 国产黄色小视频在线观看| 色综合婷婷激情| 亚洲精品色激情综合| 国产精品 欧美亚洲| 国产在线观看jvid| xxxwww97欧美| 天天躁夜夜躁狠狠躁躁| 制服诱惑二区| 中文亚洲av片在线观看爽| 日日干狠狠操夜夜爽| 一区二区三区国产精品乱码| 欧美黄色片欧美黄色片| 久久久久精品国产欧美久久久| 亚洲欧美日韩无卡精品| 一进一出抽搐动态| 欧美成人午夜精品| cao死你这个sao货| 国产av麻豆久久久久久久| 国产蜜桃级精品一区二区三区| 老司机午夜十八禁免费视频| 国产99白浆流出| 在线观看美女被高潮喷水网站 | 18禁裸乳无遮挡免费网站照片| 国产伦在线观看视频一区| 国产欧美日韩精品亚洲av| 欧美大码av| 日本 av在线| 桃色一区二区三区在线观看| 国产爱豆传媒在线观看 | 无限看片的www在线观看| 国产精品日韩av在线免费观看| 久久久久久久久免费视频了| 国产精品久久久人人做人人爽| 欧美丝袜亚洲另类 | 国产成人影院久久av| 欧美激情久久久久久爽电影| 听说在线观看完整版免费高清| 免费在线观看日本一区| 五月伊人婷婷丁香| aaaaa片日本免费| 伊人久久大香线蕉亚洲五| 国产亚洲av嫩草精品影院| 三级男女做爰猛烈吃奶摸视频| 欧美zozozo另类| 亚洲自偷自拍图片 自拍| 中文字幕熟女人妻在线| 国产精品电影一区二区三区| 18禁国产床啪视频网站| 精品久久久久久成人av| 精华霜和精华液先用哪个| 非洲黑人性xxxx精品又粗又长| 一区二区三区高清视频在线| 99riav亚洲国产免费| 757午夜福利合集在线观看| 午夜亚洲福利在线播放| 亚洲国产欧美网| 国产午夜精品久久久久久| 欧美久久黑人一区二区| 久久99热这里只有精品18| 久久久国产成人精品二区| 人人妻人人澡欧美一区二区| 黄色片一级片一级黄色片| 麻豆久久精品国产亚洲av| 后天国语完整版免费观看| 动漫黄色视频在线观看| 亚洲五月婷婷丁香| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲五月天丁香| 成人高潮视频无遮挡免费网站| 叶爱在线成人免费视频播放| 国产伦在线观看视频一区| 最近视频中文字幕2019在线8| 99精品在免费线老司机午夜| 一二三四社区在线视频社区8| 精品久久久久久久久久免费视频| 亚洲激情在线av| 国产精品 欧美亚洲| 国产成人系列免费观看| 最新在线观看一区二区三区| 亚洲真实伦在线观看| 在线看三级毛片| 欧美日本亚洲视频在线播放| 在线看三级毛片| 欧美日本亚洲视频在线播放| 亚洲成a人片在线一区二区| 久久人人精品亚洲av| 成人国语在线视频| 村上凉子中文字幕在线| 首页视频小说图片口味搜索| 国产成人系列免费观看| 亚洲av电影不卡..在线观看| 国产一区在线观看成人免费| 人妻久久中文字幕网| 亚洲成人中文字幕在线播放| 国产成人欧美在线观看| 国产成人系列免费观看| 欧美成人免费av一区二区三区| 亚洲精品久久成人aⅴ小说| 成人18禁在线播放| 一进一出好大好爽视频| 久久午夜综合久久蜜桃| 五月玫瑰六月丁香| 亚洲av成人精品一区久久| 老汉色av国产亚洲站长工具| 国产激情欧美一区二区| 91成年电影在线观看| 午夜日韩欧美国产| 黄频高清免费视频| 国内精品一区二区在线观看| 亚洲黑人精品在线| 国产精品亚洲一级av第二区| 91国产中文字幕| 人妻夜夜爽99麻豆av| 一级作爱视频免费观看| 天天一区二区日本电影三级| 黄频高清免费视频| 曰老女人黄片| 十八禁人妻一区二区| 亚洲成人久久性| 久久这里只有精品19| 在线视频色国产色| 精品日产1卡2卡| 99久久精品热视频| 淫妇啪啪啪对白视频| 19禁男女啪啪无遮挡网站| 国产成人精品久久二区二区免费| 中文字幕熟女人妻在线| 老司机深夜福利视频在线观看| а√天堂www在线а√下载| av免费在线观看网站| 禁无遮挡网站| 好男人电影高清在线观看| 久久久久性生活片| 精华霜和精华液先用哪个| 黄色丝袜av网址大全| 亚洲精品美女久久久久99蜜臀|